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UNE  METHODE  DE  CONSTRUCTION  D’ELEMENTS  FINIS  RATIONNELS 

DE  CLASSE  DANS  IR^ 


K.  Nicolas  AGBEVE 


{UniversHe  de  Nante^ 
3/135,  Place  de  Venise 
59650  Villeneuve  d’Ascq 


Soient  9.  tin  ouvert  borne  de  et  5  =  {o,-,  z  =  1,  ■  •  - ,  n}  un  ensemble  de 
points  isoles  de  0.  Notons  T  una triangulation  basee  sur  5,  c’est  a  dire  dont  les 
sommets  des  simplexes  sont  dans  S. 

Pour  construire  une  fonction  de  classe  avec  des  donnees  d’interpolation 
sur  les  points  de  S,  on  pent  utiliser  des  methodes  globales,  (  Duchon,  Bouhamidi,...), 
des  methodes  radiales,  (Dyn  Sz  Levin,  Powel,...)  et  aussi  des  fonctions  splines 
polynomiales. 

Si  les  deux  premieres  methodes  ont  le  desavantage  d’etre  difficiles  a  mettre 
au  point,  la  demiere  elle,  engendre  des  couts  eleves  a  I’utilisation. 

Le  but  vise  ici,  est  la  construction  des  elements  finis  de  meme  regularite 
que  celle  que  I’on  pent  avoir  avec  les  methodes  de  splines  polynomiales  mais  en 
utilisant  des  fonctions  rationnelles  pour  reduire  le  nombre  de  degres  de  liberte. 
La  methode  consiste  a  chercher  un  element  fini  sous  la  forme  d’un  interpolant  V 
qui  s’ecrit  :  V  =  C  +  %  ou  C  est  une  fonctionnelle  polynomiale  pure  et  U  une 
fonctionnelle  rationnelle  pure  c’est  a  dire  irreductible.  L’interpolant  71  apparaitra 
alors  comme  tm  correcteur  de  I’interpolant  C  ,  done  doit  conserver  les  proprietes 
de  C. 


RATE  OF  CONVERGENCE  OF  INTERPOLATING  MEANS 
WITH  EXPONENTIAL-TYPE  WEIGHTS 

G.  ALLASIA  -R.  BESENGHI 

Dipartimento  di  Matematica,  Universita  di  Torino,  Via  Carlo  Alberto  10, 1-10123  Torino,  Italy. 


Given  a  set  5„  C  [a,  6]  C  IR  of  irregularly-spaced  points,  we  consider  the  interpolation  of  a  function 
/  G  C’’[a,  6]  by  the  weighted  arithmetic  mean 

n 

/i  =  /r(i;r,S„)  =  ^  Tk(x;f,r)wk(.x;SJ, 

k=l 

that  is  a  positive  linear  functional.  In  order  to  increase  the  degree  of  approximation  of  /x  to  /,  we 
use  in  the  operator  /i  ,  instead  of  the  value  /(x*),  (k  =  1,2,..,  n),  of  the  interpolated  function,  the 
truncated  Taylor  expansion  of  /  up  to  the  derivative  of  order  r  about  the  point  Xk 

Tk  =  Tk(x-, /,  r)  =  ^  (x  -  XkY. 


The  weights  multiplying  the  Taylor  expansions  are 


Wk  =  njkCx,  Sn)  = 


|x-xt|  ^exp(-^|x- xtl) 
ELi  1=^  -  Xh\-^exp(-l3\x  -  Xftl)  ’ 


with  q  and  p  nonnegative  real  numbers  .  They  satisfy  the  properties 


Wkix;Sn)>0,  y^u;t(x;g„)  =  1, 

k=l 

depend  on  the  data  points  and  are  approximately  exponentially  decreasing  with  the  distance.  These 
weights  are  interesting,  both  theoretically  and  practically,  and  include  noteworthy  particular  cases. 
At  first  we  prove  that  the  derivatives  D’wk,  (s  =  0, 1 ,  ...r),  of  the  weights  vanish  at  the  interpolation 
nodes.  As  a  conseguence,  we  obtain  that  the  derivatives  of  the  interpolating  operator  and  the 
corresponding  derivatives  D’f  of  the  interpolated  function  are  equal  at  the  nodes. 

We  then  examine  the  degree  of  approximation  of  /ir  to  /,  giving  upper  bounds  for  the  difference 
—  D^f\  by  means  of  the  modulus  of  continuity  of  D^f. 
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Abstract: 

We  present  a  new  vector  quasi-interpolant  based  on  a  vector  spline  func¬ 
tion  introduced  by  the  authors  (  J.  Approx.  Th.,  Vol.  67,  No.  1,  Oct.  91). 
Following  a  recent  approach  used  in  quasi-interpolation  (C.  Rabut,  These 
d’Etat,  90),  the  elementary  matrix  approximant  is  obtained  by  convolution 
of  a  discretized  version  of  the  partial  differential  operator,  which  caracterize 
the  vector  spline  function,  and  the  exact  fundamental  solution  associated  to 
the  operator.  Some  general  regularity  properties  of  the  elementary  matrix 
are  given  along  with  results  on  convergence  of  the  quasi-interpolant.  Specific 
features  related  to  the  vector  structure  ot  the  problem  are  showed,  rinaliy 
numerical  exemples  of  the  elementary  matrix  are  given. 
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Authors:  I  J  Anderson,  M  G  Cox,  S  Harbour  and  J  C  Mason 

Speaker:  J  C  Mason 

Address:  Applied  and  Computational  Mathematics  Group, 

RMCS  (Cranfield),  Shrivenham,  Swindon  SN6  SLA,  England 


ABSTRACT: 

Structured  data,  such  as  data  whose  ordinates  lie  on  or  near  a  set  of  parallel  lines  or  a 
family  of  curves,  should  be  fitted  in  such  a  way  as  to  exploit  the  data  structure  -  even 
when  the  ordinates  are  irregularly  spaced  and  in  different  positions  on  each  line/ curve. 
Existing  algorithms  of  this  type,  which  effectively  separate  the  variables,  include  the  1963 
algorithm  of  C  W  Clenshaw  and  J  G  Hayes  for  Chebyshev  polynomial  approximation  of 
lines  of  data,  the  1990  algorithm  of  J  C  Mason  and  R  P  Bennell  which  effectively  extends 
this  to  a  procedure  for  curves  of  data,  and  the  1993  algorithm  (to  appear  in  "Algorithms 
for  Approximation  3"  -  Proceedings  of  1992  NATO  Symposium,  Oxford)  of  I J  Anderson, 
M  G  Cox,  and  J  C  Mason  which  iteratively  interpolates  a  B-spline  through  data  near  to 
lines.  These  algorithms  are  now  being  extended  in  two  principal  ways:  first,  B-splines 
are  being  fitted  in  a  least  squares  sense  to  (overspecified)  data  on  or  near  lines/ curves, 
and  second,  parallel  processing  procedures  are  being  investigated  and  implemented  for 
exploiting  the  separability  of  the  data/algorithms.  It  is  relatively  straightforward  to 
envisage  the  principles  of  parallel  algorithms  for  these  problems,  but  it  is  more  important 
and  difficult  to  carry  this  out  in  practice  and  to  assess  the  communications  and 
organisational  overheads  for  specific  parallel  computer  platforms. 
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Smoothing  Noisy  Data  with  Coiflets 
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ABSTRACT 

This  paper  is  concerned  with  an  orthogonal  wavelet  series  re¬ 
gression  estimator  of  an  unknown  smooth  regression  function  ob¬ 
served  with  noise  on  a  bounded  interval.  The  method  is  based  on 
applying  results  of  the  recently  developed  theory  of  wavelets  and 
uses  the  specific  asymptotic  interpolating  properties  of  the  wavelet 
approximation  generated  by  a  particular  wavelet  basis,  Daubechie  s 
coiflets.  Conditions  are  given  for  the  estimator  to  attain  optimal 
convergence  rates  in  the  integrated  mean  square  sense  as  the  sam¬ 
ple  size  increases  to  infinity.  Moreover,  the  estimator  is  shown  to 
be  pointwise  consistent  and  asymptotically  normal.  The  numerical 
implementation  of  the  estimation  procedure  relies  on  the  discrete 
wavelet  transform  and  the  algorithm  for  smoothing  a  noisy  sample 
of  size  n  requires  order  0(11)  operations.  The  general  theory  is  il¬ 
lustrated  with  simulated  and  real  examples  and  a  comparison  with 
other  nonpar ametric  smoothers  is  made. 

Keywords  and  phrases:  Nonparametric  regression,  curve  smoothing,  wavelets, 
multiresolution  analysis,  splines. 


QUASI-INTERPOLANTS,  BOX-SPLINES  AND  (QUASI-)WAVELETS  [P(D)]-manifolds. 
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Let  C  =  (Cl.  -.Cn)  €  IR”  and  PijiO,  1  <  ^  <  m,  I  <  j  <  q,  he  a.  real  polynomial  iu  the  variables  Ci,  •■•,Cn 
with  constant  coefficients.  We  set  P(C)  =  (PijiC))- 

In  [1],  with  convenient  assumptions,  we  have  proved  the  existence  of  a  spline  manifold  a  relating  to 
P{D)  (where  D  =  (  -z — , ...,  —  ))  and  that  cr  in  characterized  by  a  hilbertian  (semi-)kernel  £  relating  to 

\C/Xn  J 

*pp. 

In  this  lecture,  we  define  and  study  quasi-interpolants  (resp.  box-splines)  manifolds  and  (quasi-)wavelets 
manifolds,  when  P  in  homogeneous.  They  are  obtained  by  approximating  in  a  convenient  way  the  formula 
*PPP  =  6  {6  is  the  Dirac  functional)  on  a  sequence  of  regular  meshes  of  IR". 

We  give  some  examples. 


[1]  Spline  manifolds,  Curves  and  surfaces,  Chamonix  1990,  pp.  1-4,  Academic  Press  N.Y.  1990 


Affine  Classification  of  Bezier  Cubics 


M.  L.  Baart 
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Abstract: 

Cubics  occupy  an  intermediate  position  in  the  hierarchy  of  algebraic  curves.  They  allow 
considerably  more  freedom  than  conics,  but  limit  the  types  of  complications  that  can 
occur  for  higher  order  curves.  We  review  some  relevant  algebraic  properties  of  cubics 
in  general,  and  polynomial  cubics  in  particular.  Of  these  properties,  those  that  remain 
invariant  under  affine  transformations  are  of  particular  interest  in  applications  in 
computer  graphics.  These  affine-invariant  properties  can  be  used  to  find  a  practical 
classification  of  planar  polynomial  cubics.  As  an  application,  a  complete  partition  of  the 
parameter-plane  of  such  Bezier  cubics  are  presented.  In  conclusion  we  consider  additional 
classes  introduced  by  allowing  rational  parametrization. 

References: 
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ABSTRACT 

The  problem  of  reconstructing  a  function  from  its  line  integrals  has  arisen  indepen¬ 
dently  in  many  scientific  fields,  astrophysics,  geophysics,  and  particularly  in  medicine. 
In  practice  there  is  an  unknown  compact-supported  function  /  defined  in  B?  ,  repre¬ 
senting  relative  tissue  density,  a  finite  number  of  measurements  are  given  approximating 
integrals  of  /  along  lines  and  we  have  to  reconstruct  /  from  these  data. 

Two  approaches  have  been  used  to  solve  this  problem,  the  first  one  deals  with 
methods  aiming  to  approximate  the  inverse  of  the  transformation  defined  by  the  line 
integrals  or  Radon  Transform;  these  type  of  methods  are  colled  transform  methods 
and  one  of  them,  convolution-backprojection,  has  been  universally  adopted  for  medical 
purposes.  The  sencond  approach  consists  of  discretizing  the  model  at  an  early  stage 
of  the  process,  leading  to  a  finite  dimensional  algebraic  system  of  linear  equations  for 
which  a  solution  has  to  be  choosen  using  some  iterative  methods  .  Both  these  methods 
involve  large  amount  of  data  when  dealing  with  practical  applications. 

In  this  work  we  have  studied  and  experimented  the  use  of  orthogonal  and  biortog- 
onal  wavelet  and  wavelet  packet  bases  with  the  aim  to  construct  efficient  parallel  algo¬ 
rithms  for  the  solution  of  these  reconstruction  problems.  Numerical  results  obtained  on 
the  distributed  memory  system  iPSC/860  hypercube  axe  presented. 
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ABSTRACT 

Wavelet  decomposition  has  recently  been  used  to  convert  dense  matrices  of  a  cer¬ 
tain  type  into  higly  sparse  matrices,  with  consequent  reduction  of  the  computational 
cost  applying  them  to  arbitrary  vectors  [1],  [2].  Moreover,  efEcient  nonstandard  matrix 
multiplication  algorithms  have  been  obtained  by  doing  a  wavelet  packet  dec  omposition 
of  the  matrices  and  setting  all  coefficients  to  zero  whose  absolute  values  are  less  than 
a  fixed  tollerance  [3].  In  this  paper  we  take  advantage  of  the  typical  block-structure 
presented  by  the  wavelet  packet  decomposed  matrices  in  order  to  construct  parallel 
numerical  algorithms  whose  hight  efficiency  is  mostly  due  to  the  strong  diagonal  dom¬ 
inance  of  the  trasformed  matrices  and  the  same  structure  of  the  different  blocks.  An 
extensive  experimentation  taken  on  a  distributed  memory  multiprocessor  iPSC/860  is 
presented, 
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Given  a  point  on  a  B-spline  curve,  one  sometimes  wishes  to  find  which  control  vertices 
most  affect  that  point.  The  question  arises:  if  we  find  the  control  vertex  with  maximum 
influence,  will  vertices  with  indices  near  the  index  of  the  vertex  of  maximum  influence 
aricici  the  curve  more  than  vertices  with  indices  farther  away?  More  specifically,  given  a 
sequence  of  B-splines  {Nj},  and  a  domain  value  t,  if  k  is  an  index  such  that  Nkit)  >  Nj{t) 
for  all  j,  is  it  true  that  Nj(i)  <  iVj-)-i(i)  for  all  j  <  k,  and  Nj{t)  >  for  all  j  >  k? 

This  talk  will  show  that  this  property  does  hold  for  many  types  of  B-splines,  but  not  for 
all. 


Object  Oriented  Spline  Software 
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Abstract 

Object  oriented  programming  provides  software  with  a  facility  for  imitating  the  math¬ 
ematician’s  tool  of  theoretical  abstraction.  The  communality  of  a  collection  of  related 
entities  can  be  described  in  terms  of  a  base  class,  which  may  be  added  to  an  object  ori¬ 
ented  language  as  a  new  data  type.  Each  concrete  entity  of  the  collection  can  be  described 
as  a  derived  class  by  specifying  whatever  features  distinguish  it  from  the  base,  and  the 
entity  can  be  added  to  the  language  as  a  further  data  type.  Algorithms  written  in  terms 
of  the  base  type  can  be  applied  to  any  of  the  derived  types. 

This  talk  will  introduce  the  concepts  of  object  oriented  programming  and  illustrate  the  use 
of  its  facilities  for  abstraction  by  presenting  class  designs  for  spline  bases,  spline  functions 
(curves  and  surfaces),  and  spline  knot  insertion.  Some  examples  of  the  use  of  these  classes, 
of  their  further  extension,  or  of  alternative  approaches  to  class  design  will  be  given  in 
associated  talks  scheduled  for  the  Workshop  on  Software  Infrastructure  for  Computer  Aided 
Design. 
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Abstract 

This  workshop  will  cover  some  issues  on  the  design,  implementation,  and  use  of  software  for 
spline  curves,  surfaces,  and  finite  elements.  Techniques  for  object  oriented  programming 
will  be  a  primary  focus  of  the  talks.  Bartels  will  cover  the  basic  concepts  of  object  oriented 
programming,  specifically  its  capacity  to  support  the  abstraction  of  algorithms  and  data 
types,  and  will  offer  basic  examples  in  the  design  of  classes  for  spline  bases,  refinement, 
and  curves  and  surfaces.  Greiner  and  Slusallek  will  show  how  the  abstraction  of  surface 
software  facilitates  the  implementation  of  surface  blending.  Hickey  and  Seidel  will  cover  the 
use  of  object  oriented  classes  and  symbolic  algebra  to  build  a  general  system  for  connection 
matrix  splines.  Vermeulen  will  cover  the  use  of  splines  and  object  oriented  techniques  to 
tackle  the  geometry  and  finite  element  analysis  of  physical  structures. 
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An  Algorithm  for  Bivariate 
Simultaneous  Polynomial  Approximation 


G.  Baszenski 
Dortmund,  Germany 


We  consider  the  problem  of  approximating  a  smooth  function  f(x,y) 
on  the  square  [-1,1]  x  [-1,1]  by  a  polynomial  p(x,y)  such  that  derivatives 
of  f(x,y)  are  approximated  by  derivatives  of  p(x,y).  By  blending 
interpolation  mediods  we  construct  a  preapproximation  based  on  the 
boundary  values  of  f(x,y).  The  remainder  fimction  is  then  interpolated 
on  a  Chebyshev  node  grid.  Test  computation  results  show  that  this 
procedure  yields  an  approximant  with  high  convergence  orders  for  the 
function  values  and  derivatives. 
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On  polynomial  functions  defining  the 
geometrical  continuity  between  two  (SBR). 


J.L.  Bauchat 
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This  talk  develops  works  about  the  and  geometrical  continuity  be¬ 
tween  rational  surfaces  put  into  {SBR)  form.([l]  ,  [2]  ,  [3])  . 

The  study  of  G^  and  G^  geometrical  continuity  between  two  surfaces  which 
are  Iin  -  projections  of  polynomial  surfaces  ,  leads  to  necessary  and  sufficient 
conditions  pertaining  to  the  parametrizations  C{u,v)  and  TZ{u,v)  of  the 
surfaces  “above” .They  are  expressed  by  polynomial  functions  {/i,/2,  /a,  A}, 
for  the  G^  continuity  ,  and  {gi , 92 , 93, 94}  ,  for  the  G"^  continuity,  giving 
relations  between  the  last  rows  of  the  massic  net  controlling  Si  and  the  first 
of  the  massic  net  of  Sr- 

An  example  deals  with  the  G^  geometrical  connection  between  two  quadrics: 
the  obtained  functions  show  how  the  massic  vectors  of  Si  and  Sr  are  inter¬ 
related. 

Then  we  analyse  the  constaints  to  be  satisfied  by  the  coefficients  of  the 
polynomial  functions  {/i ,  /z,  /a,  /i}  in  order  to  guarantee  the  G^  geometrical 
continuity  between  a  given  {SBR)  surface  Si  and  a  new  {SBR)  surface  to 
be  construct  .  The  results  are  illustrated  by  an  example. 

[1]  Fiorot ,  J.C. ,  P.  Jeannin  ,  Courbes  et  Surfaces  rationnelles  .  Applications 
a  la  C.A.O.  ,  Coll.  R.M.A.  12  ,  Masson  ,  Paris  1989. 

[2]  Bauchat  ,  J.L.  ,  Etude  du  raccordement  geometrique  G^  de  surfaces  rsi- 
tionnelles  mises  sous  forme  {SBR),  These  Universite  des  Sciences  et  Tech¬ 
nologies  de  LiUe  ,  Laboratoire  A.N.O.  ,  mars  1992. 

[3]  Bauchat, J.L.  ,  J.C.  Fiorot  ,  P.  Jeannin  ,  G^  and  continuity  of  {SBR) 
surfaces  in  Curves  and  Surfaces  ,  P.J.  Laurent  ,  A.  le  Mehaute  and  L.L. 
Schumaker  (eds)  ,  Academic  Press  ,  New  York  1991. 

‘Laboratoire  de  conception  et  de  realisation  mecaniques  assistees  par  ordinateur 


32. 


CONSTRUCTION  OF  ITERATION  FUNCTIONS 
FOR  THE  SIMULTANEOUS  COMPUTATION  OF  THE  SOLUTIONS 
OF  SYSTEMS  OF  ALGEBRAIC  EQUATIONS 

Anne-Mercedes  BELLIDO 
L.A.N.  U.P.S.  31062  Toulouse  Cedex 
E-Mail  :  bellido@cix.cict.fr 

The  material  of  this  talk  is  part  of  a  thesis  advised  by  J.C.  Yakoubsohn  ([2]).  Our  purpose 
was  to  generahze  to  the  multivariate  case  the  well-known  Weierstrass’  iteration  function  (iOi), 
for  the  simultaneous  computation  of  the  roots  of  a  monic  polynomial  P  of  degree  d  in  the  case 
of  simple  roots  ;  this  function  is  attractive  not  only  for  its  local  quadratic  convergence  but  above 
all  because  numerical  experimentation  incited  a  few  authors  to  advance  the  conjecture  that  it 
converges  with  almost  every  given  starting  point  (i4]). 

The  unified  approach  of  the  methods  of  simultaneous  computation  of  roots  of  generalized 
polynomials  developed  by  A.  Frommer  ([3|)  particularly  well  fitted  with  our  aim. 

Our  contribution  to  this  approach  in  the  univariate  case  (il])  essentially  rehed  on  three  points  ; 

-  the  first  one  is  the  use  of  spaces  V  containing  the  functions  of  the  studied  system 
whose  dimension  D  may  be  greater  than  d  -r  1-  by  focussing  our  attention  on  the  part  played 
by  the  kernel  of  a  linear  functionnal  A. 

-  the  second  one  is  the  introduction  of  interpolation  basis  of  this  kernel  (depending 
on  the  iteration  points  )  to  construct  in  practice  the  iteration  function,  in  the  case  in  which, 
from  a  theoretical  point  of  view,  basis  are  supposed  to  be  independent  of  these  points. 

-  the  third  one  is  the  extension  of  the  theoretical  construction  of  iteration  functions 
with  basis  depending  on  the  iteration  points. 

Here,  we  show  how  the  two  first  points  are  the  keys  to  construct  iteration  functions  for  the 
simultaneous  computation  of  the  solutions  of  a  system  of  algebraic  equations  (5)  /'^  =  0 
A  =  1, . . . ,  n,  in  the  case  of  zero-dimensional  sets  of  solutions. 

We  obtain  in  particular,  to  compute  the  d  solutions  of  (S),  a  «  generalized  Weierstrass’  iteration 
function  »  from  into  C"'^  in  the  form  X  =  {Xi, . . .  ,Xd)  with  : 


1  -1 


r,(r)  =  . .  .  ,  =  r(‘)  -  ^j(0  (T(t))J  •  F  (r^')) 

^'here  is  the  transposed  vector  of  J^‘^(T(r))  the  Jacobian 


matrix  evaluated  at  the  point  of  a  system  T(t)  of  n  equations  and  n  unknowns,  depending 
on  T  and  identical  to  (S)  at  the  solution  point. 

We  sum  up  the  first  experimentations  and  show  how  to  improve  convergence  results  in  order 
to  obtain  a  conjecture  of  global  convergence  in  the  real  case. 

We  end  applying  one  of  these  algorithms  to  the  solution  of  a  problem  of  robotic. 
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3D  Generalized  Voronoi  Diagram 
E.  BERTIN  and  J.M.  CHASSERY 

Laboratoire  TIMC-IMAG,  URA  D  1618,  Universite  Joseph  FOURIER,GRENOBLE. 

Abstract:  This  paper  describes  a  study  of  the  3D  Generalized  Voronoi  Diagram  for  polyhedra 
(3D  GVD)  as  a  set  of  quadric  surfaces  and  planes.  The  separator  between  two  distinct  objects 
Obji  and  Obj2  is  defined  as  the  set  of  points  equidistant  to  Obji  and  Obj2.  The  3D  GVD  is  the 
union  of  such  separators.  Two  main  research  axis  will  be  presented.  The  first  axis  is  a  theoritical 
study  of  the  equations  of  separators  between  3D  polyhedral  objects.  This  study  will  be  illustrated 
with  a  first  theorem.  The  second  axis  involves  presentation  of  an  algorithm  computing  an 
approximated  solution  of  the  3D  GVD  based  on  uses  of  3D  Voronoi  Diagram  for  Ponctual  seeds. 
The  convergence  of  this  approach  is  provided  in  a  second  theorem. 

Theorem  1:  Let  S  be  a  set  of  polyhedral  objects,  the  separators  between  the  objects  of  S  are 
composed  of  following  quadric  surfaces:  Paraboloid  of  revolution:  x^  +  y^  +  z  =  1,  Hyperbolic 
paraboloid:  x^  -  y2  =  z.  Parabolic  surface:  x^  =  z.  Quadric  which  equation  is  x^  +  y2  -  z^  =  0  and 
r.une:  a.x  +  b.y  +  c.z  +  d  =  0. 

A  second  result  concerns  computation  of  the  3D  GVD  of  a  set  of  polyhedric  volumes.  Our 
approach  is  based  on  discretization  of  the  polyhedral  surfaces  and  utilisation  of  the  3D  Voronoi 
Diagram  for  Points.  The  result  obtained  with  this  algorithm  is  an  approximation  of  the  3D  GVD 
of  the  set  of  objects  and  the  proof  of  this  algorithm  is  based  on  the  following  theorem. 

Theorem  2:  The  discretisation  of  the  3D  GVD  by  the  3D  Voronoi  Diagram  for  Points  converges 
toward  the  3D  GVD  when  the  step  of  discretisation  converges  towards  zero. 

Reference:  Hu  H.T.,  Diagramme  de  Voronoi'  Generalise  Pour  un  ensemble  de  Polygones, 


3D  Voronoi  Diagram  for  one  point  and  a  square. 


3D  Voronoi  Diagram  for  the  border  of  a  square  and  one  segment. 
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Spline  under  tension,  interpolants  of  scattered  data 


A.  Bouhamidi  and  A.Le  Mehaute 
Abstract 

We  recall  briefly  some  properties  of  -spline  and  we  show  how  these  functions  enable  us  to  charac¬ 
terize  thin  plate  spline  under  tension.  Some  examples  are  presented,  showing  the  effect  of  the  parameter  of 
tension. 

The  concept  of  spline  curve  under  tension  was  introduced  by  Schweikert  and  later  D.  Franke  investigates 
a  scheme  for  the  representation  of  surfaces  which  incorporates  the  same  idea  of  tension..  The  method  of  Franke 
presents  an  elegant  aspect  but,  unlike  Duchon’s  thin  plates  and  other  radial  basis  functions  methods, it  is 
lacking  in  some  mathematical  theory  and  presents  also  a  problem  of  a  numerical  matter,  because  it  involves 
Bessel  functions  that  need  to  be  approximate  for  each  entry  of  a  full  matrix.  Using  another  differential 
operator,  and  the  associated  reproducing  kernel,  we  are  able  to  construct  surface  spline  under  tension  as  a 
particular  case  of  -spline,  a  generalization  of  thin  plates  splines.  We  recall  briefly  some  properties  of 
these  functions  and  show  how  it  works  on  some  examples  The  data,  with  scattered  data  points  are  taken 
from  Franke  and  have  respectively  23  points,  33  points  and  130  points. 

All  of  the  examples  illustrate  the  behavior  of  the  surfaces  for  various  values  of  the  tension  parameter. 
We  note  that  as  the  value  tension  is  increased,  the  overshoot  becomes  less,  and  finally  disappears. 

The  same  idea,  restricting  to  the  univariate  case,  give  a  theoretical  framework  (Hilbert  spaces  and 
reproducing  kernels)  for  spline  curves  under  tension. 

A.  Bouhamidi  and  A.  Le  Mehaute 
Departement  de  Mathematiques  et  d’lnformatique 
Faculte  des  Sciences,  Universite  de  Nantes 
2  rue  de  la  Houssiniere 
44072  NANTES  CedexOS  -  FRANCE 
email :  alm@  ci  cb  .fr 
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SINGULARITY  DETECTION  FROM  NOISED  DATA  WITH  WAVELETS 


M.  T.  Bozzini^  F.  De  Tisi^  M.  Rossini^ 

We  have  studied  the  problem  of  detecting  the  locations  of  irregularity  points  of  a 
function  y  =  f{x)  on  the  base  of  a  data  set  corrupted  by  noise: 

yi  —  fi Ci  z 

Cj-  r.v.  idd  with  E{ei)  =  0  =  a 

The  problem  is  particulary  complex  because  it  is  difficult  to  distinguish  discontinuities  of 
the  function  from  those  caused  by  the  noise. 

In  literature  there  are  some  papers  about  the  matter.  In  particular  we  have  considered 
[3]  where  Lee  has  tried  to  solve  the  problem  using  the  Gabor  window  together  with  a 
statistical  test  on  the  error  behaviour  while,  in  more  recent  papers  [1]  and  [2],  the  authors 
have  tried  to  detect  that  points,  studying  the  evolution  of  the  wavelet  transform  across 
scales. 

We  present  a  method  that  intends  to  detect  such  irregularity  points,  using  the  wavelet 
transform  together  with  a  statistical  study  of  the  error  spectrum. 
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An  introduction  to  Fade  approximations 

C.  Brezinski 


Laboratoire  d’ Analyse  Numeriqne  et  d’ Optimisation 


UFR  lEEA  -  M3 

Universite  des  Sciences  et  Technologies  de  Lille 
59655  Villeneuve  d’Ascq  Cedex,  France 


This  talk  is  an  introduction  to  Fade  approximations. 

Let  /  be  a  formal  power  series.  A  Fade  approximant  is  a  rational  func¬ 
tion  whose  power  series  expansion  agrees  with  that  of  /  as  far  as  possible. 
After  reviewing  some  of  their  properties,  the  problem  of  their  recursive 
computation  will  be  treated  by  introducing  formal  orthogonal  polynomials. 
Then,  their  convergence  will  be  studied. 

Fade  approximants  have  very  many  applications  ranging  from  the  approx¬ 
imation  of  functions  to  the  inversion  of  the  Laplace  transform  and  to  the 
construction  of  A-acceptable  approximations  to  the  exponential  function. 
Some  of  these  applications  will  be  discussed  and  a  practical  method  for 
estimating  the  error  will  be  given. 

Fade  approximants  have  been  generalized  in  several  directions  including 
the  multivariate  case,  series  in  a  non-commutative  algebra,  Fade-type  and 
partial  Fade  approximants,  vector  Fade  approximation  and  so  on.  Some  of 
these  generalizations  will  be  discussed. 


References 
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3? 


Scattered  Data  Interpolation  over  the  Sphere 


by  Jeff  Brown 


In  1980  Robin  Sibson  developed  a  scattered  data  inteipolant  for  euclidean  spaces 
which  has  local  support,  i.e.  the  interpolated  function  value  at  a  point  p  depends  only  on 
the  function  values  at  data  sites  near  p.  Let  {pj  be  the  data  sites.  Sibson’s  inteipolant  is 
based  on  a  set  of  coordinate  functions  Sf  which  satisfy  the  identities 

=  l  P  =  X^/(P)Pi' 

The  inteipolant  is  defined  throughout  the  convex  hull  of  the  data  sites  and  is  C  ^  except  at 
the  data  sites,  where  it  is  C°.  (In  the  planar  case,  the  speaker  has  extrapolated  Sibson’s 
—-rriinates  to  the  entire  plane  in  such  a  way  that  they  are  still  C  ^  almost  everywhere  and 
Sibson’s  identities  are  preserved.) 

The  basic  Sibson  coordinates  have  been  used .  as  building  blocks  for  other 
inteipolation  schemes  by  Sibson  and  by  Farin. 

The  speaker  will  define  coordinate  systems  for  points  scattered  on  the  sphere  which 
are  analogous  to  Sibson’s  coordinates  in  the  plane.  These  coordinates  provide  an 
interpolant  over  the  sphere  which  has  local  support  and  is  C  \  except  at  the  data  sites, 
where  it  is  C°.  As  in  the  planar  case,  the  local  support  property  makes  these  coordinate 
functions  good  candidates  for  use  in  other  interpolation  schemes. 
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Box  Spline  Interpolation  and 
Solution  of  the  Poisson  Problem 

Gisela  Brumme 

University  of  Rostock(Germany) 

The  Fourier  technique  is  very  convenient  to  describe  cardinal  interpolation 
by  spline  functions  generated  by  integer  translates  of  a  box  spline  (see  K.  Jet- 
ter,  in:  Approximation  Theory  VII,  Academic  Press  1992,  pp.  131-161). 
Concerning  numerical  applications  the  case  of  periodic  interpolation  based 
on  periodization  of  the  cardinal  spline  function  is  more  important  than  the 
cardinal  case.  Even  the  point  of  Fourier  analytical  view  shows  the  close 
connection  between  cardinal  and  periodic  interpolation. 

In  this  talk  we  construct  solutions  of  the  Poisson  equation  on  the  unit  square 
by  box  spline  translates.  Our  method  is  based  on  new  applications  of  bi¬ 
variate  periodic  spline  interpolation  and  attenuation  factors  (see  H.  G.  ter 
Morsche,  in:  Topics  in  Multivariate  Approximation,  Academic  Press  1987, 
pp.  165-174;  M.  Gutknecht,  Numer.  Math.  51  (1987),  pp.  615-629). 

We  obtain  solutions  of  the  Poisson  problem  in  form  of  Fourier  series  by  two 
different  methods.  Both  procedures  can  be  considered  as  a  generalization  of 
fast  Poisson  solvers.  By  our  first  method  we  derive  a  kind  of  eigenfunction 
expansion.  The  second  approach  is  a  collocation  method. 

The  Fourier  representation  of  our  solutions  provides  very  efficient  numerical 
algorithms  based  on  FFT-techniques.  Numerical  test  results  will  be  presen¬ 
ted. 
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An  energetic  paradigm  for  multimodal  images  fusion 
Lionel  Bmnie,  Philippe  Cinquin 


Laboratoire  d'informatique  medicale 
Faculte  de  medecine,  38700  La  Tronche  (France) 
tel:  (33)  76  63  71  00,  ext.  84  60 
fax:  (33)  76  51  86  67 
e-mail:  lionel@timb.imag.ff 


In  many  applications,  and  especially  in  biomedicine,  the  analysis  and  the  interpretation  of  a  three-dimensional 
scene  can  rarely  be  performed  using  a  single  source  of  information.  In  such  applications,  merging  the 
informations  given  by  various  sensors  can  provide  a  powerful  alternative.  A  good  example  of  this  procedure  is 
the  fusion  of  scanner  pre-operative  3D  images  and  X-ray  intra-operative  2D  images  which  allows  to  infer, 
during  the  intervention,  a  3D  representation  of  the  surgical  scene. 

The  matching  of  two  3D  models  of  a  same  structure  or  the  registration  of  a  3D  model  with  a  set  of  3D  points 
has  been  widely  studied  in  the  literature  and  is  today  pretty  well  solved.  On  the  contrary,  the  estimation  of  the 
pose  of  a  smooth  object  by  matching  2D  images  with  a  3D  a  priori  model  of  this  object  remains  a  widely  open 
problem. 

We  show  that  this  problem  can  be  solved  by  an  iterative  procedure  based  on  a  mechanical  model.  We  will  first 
suppose  that  the  contour  of  the  object  has  been  extracted  in  all  the  2D  images  (at  least  two),  called  "dynamic" 
images.  In  the  "dynamic"  pose  of  the  object,  all  the  projection  lines  (i.e.  the  lines  joining  a  pixel  to  the  focal 
point  of  its  image)  are  tangent  to  the  surface  of  the  3D  model.  For  complex  objects,  we  can  reasonably  suppose 
that  in  a  large  neighborhood,  this  dynamic  attitude  is  the  only  one  leading  to  such  a  geometric  configuration. 
Then,  in  a  given  pose,  we  consider,  for  any  projection  line,  the  spring  of  minimal  length  making  this  line 
tangent  to  the  surface  S.  S  being  globally  attracted  to  the  springs,  the  distance  firom  S  to  the  set  of  lines  will 
decrease.  Based  on  this  intuitive  idea,  we  have  formalized  an  iterative  procedure  in  which  the  estimation  of  the 
dynamic  pose  of  the  object  is  described  as  the  evolution  of  a  3D  model  in  a  spring-type  forces  field.  A  dual 
approach  is  to  consider  the  potential  energy  of  the  set  of  lines  with  respect  to  the  surface.  Indeed,  this  energy 
defines  a  criterion  of  goodness-of-fit  of  a  pose  estimation  which  can  be  minimized  to  infer  the  dynamic  attitude. 

Both  approaches  rely  on  the  computation  of  the  3D  Euclidean  distance  in  the  neighborhood  of  the  surface.  For 
this  purpose,  we  have  defined  several  procedures  for  computing  3D  distance  maps,  either  based  on  dynamic 
programing  or  adaptative  hierarchical  data  structures  (octrees). 

Simulations  have  showed  an  accuracy  of  1  degree  in  rotation  and  1  mm  in  translation.  Biomedical  experiments 
have  been  successfully  run  both  with  synthetic  vertebras  and  volunteers  (merging  of  MRI  data  and  X-rays). 

For  non-segmented  2D  images,  we  define  the  notion  of  matching  credibility.  Based  on  the  merging  of  the 
informations  provided  by  various  procedures  of  image  analysis  and  by  the  knowledge  of  the  context,  a 
credibility  field  is  introduced  which  defines  some  kind  of  probability  for  a  projection  line  to  be  a  contour  line 
(i.e.  based  on  a  contour  pixel).  In  a  first  su-ategy,  we  define  the  notion  of  matching  credibility  (or  attitude 
credibility)  as  a  function  of  the  credibilities  of  the  projection  lines  tangent  to  the  surface. 

Then  the  attitude  of  the  object  is  evaluated  by  maximizing  this  credibility  function  with  respect  to  the  6  attitude 
parameters.  In  a  second  strategy,  this  credibility  field  is  merged  with  the  energetic  model  to  define  a  new 
"energy"  function  which  can  be  minimized  either  by  classical  minimization  procedures  or  by  Kalman  filtering. 


Elastic  Curves  on  the  Sphere 

Guido  Brunnett 
Department  of  Mathematics 
Naval  Postgraduate  School 

Peter  E.  Crouch 

Center  for  Systems  Science  and  Engineering 
Arizona  State  University 


This  talk  deals  with  the  derivation  of  equations  suitable  for  the  compu¬ 
tation  of  elastic  curves  on  the  sphere.  To  this  end  equations  for  the  main 
invariants  of  spherical  elastic  curves  are  given.  A  new  method  for  solving  ge¬ 
ometrically  constraint  differential  equations  is  used  to  compute  the  curves  for 
given  initial  values.  A  classification  of  the  fundamental  forms  of  the  curves 
is  presented. 


Recent  Advances  in  the  Theory  of  Radial  Basis  Functions: 

Interpolation  and  Least-squares  Approximations 
Martin  Buhmann 

Department  of  Applied  Mathematics  and  Theoretical  Physics, 

University  of  Cambridge,  Silver  Street,  Cambridge  CB3  9EW . 

In  this  talk  we  shall  review  recent  results  about  approximations  to  multivariable  functions  from  radial  basis 
function  spaces.  These  are  spaces  spanned  by  radially  symmetric  functions  ^(|1  •  —xjWi)  :TV'  ^Ti,  where 
the  “centres”  Xj  are  given.  The  radial  basis  function  (f>  is  prescribed  too.  Particularly  well-known  and  useful 
examples  include  4>{r)  =  r^logr,  the  thin  plate  spline,  and  (f>{r)  =  \Jr^  -b  c^,  the  multiquadric.  We  will 
focus  on  two  approaches  to  approximating  functions  from  those  spaces:  interpolation  at  the  centres  and 
least-squares  approximations,  and  we  will  give  a  review  of  several  results  that  have  been  found  during  the 
last  two  years,  say. 
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H.  G.  Burchard  Dept,  of  Mathematics  Oklahoma  State  University  Stillwater, 
OK  74078-0613  U.S.A. 

ABSTRACT 

APPROXIMATION  WITH  AESTHETIC  CONSTRAINTS 

In  this  paper  we  show  how  to  compute  a  best  approximating  C^-  and  piecewise 
sp line-curve,  e.g.,  a  circular  spline  (with  sufficiently  small  segments),  subject  to  constraints 
for  the  curvature  to  be  positive,  montone,  and  convex  in  a  certain  discretized  sense,  with 
respect  to  certain  intrinsic  parametrizations.  A  family  of  such  conditions  is  proposed  which 
includes  ordinary  convexity  and  logarithmic  convexity.  Logarithmically  convex  curvature 
is  particularly  well  suited  for  automobile  design  lines.  With  these  constraints,  a  procedure 
is  developed  to  approximate  data  by  a  discrete  list  of  points  for  which  the  existence  of  an 
interpolating  circular  spline  with  positive,  monotone  and  approximately  logarithmically 
convex  curvature  can  be  proved.  Some  of  this  is  joint  work,  going  back  to  the  sixties 
at  GM  Research  Labs,  cf.  H.  G.  Burchard,  J.  A.  Ayers,  W.  H.  Frey,  N.  S.  Sapidis, 
Approximation  with  Aesthetic  Constraints,  Report  GMR-7814,  General  Motors  Res.  Labs. 
(1992),  to  appear  in:  Designing  Fair  Curves  and  Surfaces,  N.  S.  Sapidis,  editor, 
S.I.A.M,  Philadelphia,  PA,  (1993). 


CONSTRUCTION  OF  TRIANGULAR  AND  RATIONAL  PIECEWISE 


SURFACES 


J.C.  CANONNE 


Universite  de  Valenciennes 


First  we  will  give  a  necessary  and  sufficient  condition  so  that  the  projections  on  a 
three-dimension  space, i  =  7rc<;(£)  and  R  =  of  two  surfaces  L  and  R  in  a  four- 

dimension  space,  may  join  with  C^^-continuity. 

Then  we  will  show  how  these  conditions  can  be  explained  in  terms  of  massic  vectors 
in  the  case  of  triangular  rational  surfaces. 

And  eventually  we  will  give  an  efficient  method  to  construct  triangular  and  rational 
piecewise  surfaces. 


References: 

[1] .Canonne,J.C.,J.C.  Fiorot  ,  Une  condition  necessaire  et  suffisante  de  raccordement 
de  surfaces  rationnelles,C.R.  Acad. Sci. Paris,  312,serie  I,1991,p.  897-901. 

[2] .Canonne,J.C.,J.C.  Fiorot,Une  condition  suffisante  de  raccordement  de  surfaces  ra- 
tionnelles  triangulaires,C.R.  Acad.Sci. Paris,  313,serie  I,1991,p.  541-546. 

[3] .Farin,G.  ,  Triangular  Bernstein-Bezier  patches,  Computer-Aided  Geom.  Design  , 
3  ,1986,83-127. 

[4] . Fiorot, J.C.  and  P.Jeannin,  Courhes  et  Surfaces  Rationnelles  .  Applications  a  la 
C'AO,Coll.RMA  12  ,Masson,Paris,1989. 

[5] . Fiorot, J.C.  and  P.Jeannin,  Courbes  Splines  Rationnelles  .  ,Coll.RMA  24  ,Mas- 
son,Paxis,1992. 


Monotonicity  preserving  representations. 

J.  M.  Carnicer  and  J.M.  Pena 

Departamento  de  Matematica  Aplicada.  Universidad  de  Zaragoza. 

Edificio  de  Matematicas,  planta  la.  50009  Zaragoza.  Spain 
JMPENA@CC.UNIZAR.ES 

Abstract:  In  Computer  Aided  Geometric  Design,  shape  preserving  properties  are 
crucial  for  curve  and  surface  design.  In  this  paper  we  deal  with  shape' preserving  proper¬ 
ties  related  with  monotonicity.  In  this  sense,  a  basis  of  blending  functions  for  constructing 
a  parametric  curve  from  its  control  polygon  preserves  monotonicity  if  the  following  prop¬ 
erty  holds:  If  the  projections  of  the  control  points  into  a  line  is  a  secpience  of  ordered 
points,  then  the  projected  parametrical  curve  is  given  by  a  monotone  function.  Further¬ 
more  we  shall  impose  that  the  orientation  is  preser\'ed.  Our  goal  can  be  achieved  using 
algebraic  tools:  A  basis  of  blending  functions  is  monotonicity  preserving  if  and  only  if  the 
corresponding  collocation  matrices  transform  monotone  increasing  vectors  into  monotone 
increasing  vectors.  We  characterize  these  matrices.  In  [Karlin],  the  stochastic  transforma¬ 
tions  preserving  monotonicity  were  characterized  in  terms  of  all  the  minors  of  order  1  and 
2  of  the  matrix.  Here  we  obtain  a  simpler  test  to  prove  that  a  matrix  transforms  a  mono¬ 
tone  vector  into  a  monotone  vector.  We  also  characterize  other  related  classes  of  matrices. 
Finally,  we  study  the  particular  case  when  the  matrix  is  totally  positive  and  stochastic, 
which  corresponds  in  Computer  Aided  Geometric  Design  to  transfomiations  with  good 
shape  preserving  properties  (see  [Goodman]  and  [Carnicer  and  Pena]).  We  prove  that 
these  representations  always  preserve  monotonicity. 

[Carnicer  and  Pena]  Carnicer,  J.  M.:  Pena,  J.  M.:  Shape  preserving  representations  and 
optimality  of  the  Bernstein  basis,  to  appear  in  Advances  in  Computational  Mathematics. 

[Goodman]  Goodman,  T.N.T.:  Shape  preserving  representations  in  Mathematical  Methods 
in  CAD,  T.  Lyche  and  L.  L.  Schumaker  eds..  Academic  Press,  19S9,  pp.  3-3.3-357. 

[Karlin]  Karlin,  S.:  Total  Positivity  Stanford  University  Press,  196S. 
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3D  curve  reconstruction  from  noisy  projections 

B.  Chalmond  and  L.  Lequang 
Universite  de  Cergy-Pontoise,  Etis 

The  data  is  composed  of  a  small  number  of  projections  (three  in  our  application)  of 
a  3D  curve.  Indeed,  these  projections  are  2D  images  obtained  from  aligned  sources:  for 
each  source,  the  curve  projection  is  not  a  curve  but  is  degraded  into  a  noisy  grey  level 
image  spread  over  it,  [Ij. 

Our  approach  is  of  ’’active  curve”  type,  [2,3, 4,5].  Let  u  be  the  generic  curve  and 
g  the  data.  This  approach  consists  in  designing  an  energy  U{cj,g)  that  measures  the 
compatibihty  between  u  and  g,  then  to  compute  iteratively:  min„  U{ui,g)  starting  from 
an  initial  curve  a;(0).  For  most  of  the  applications,  the  active  curve  u{t)  obtained  at  each 
iteration  converges  to  a  local  minimum  of  Cf,  highly  dependent  of  ^<^(0).  Here,  we  have  to 
avoid  this  drawback  since  we  suppose  there  is  a  single  possible  curve  hidden  in  g. 

Our  reconstruction  task  may  be  summarized  as  follows,  w  designates  a  smooth  curved 
and  cj  +  d  is  a  deformation  of  a;  by  a;  is  a  parametrized  spline  function:  oj  =  B6  where 
''  :li0  5-spline  matrix  based  on  a  discretization  of  [0, 1].  ^  is  a  3D  random  vector  with 
a  Gibbs  distribution.  Let  be  the  £th  projection.  If  w*  =  B6*  is  the  thrue  unknown 
smooth  curve,  then  we  define: 

Ue.{u;,6,g)  =  J2UUV%u;  +  S),S,g) 

i 

P'.(.)  =  ||P'(a.  +  «)  -  +  <)  - 

where  are  ’’attracting”  points  computed  as:  sg^d{s—Ui-Si),  d  being  a  distance. 

Our  goal  is  to  estimate  6*.  We  proposed  an  iterative  algorithm  in  the  spirit  of  a 
stochastic  EM  algorithm,  [6].  If  6'  is  the  estimate  input  in  the  current  iteration,  then  the 
new  estimate  is  given  by: 

min  E[Ue{u,S,g)] 

9 

E[.]  =  J2Ueiu,6,g)P{6  \  uj,g,e')P{oj  \  gj') 

w,6 

P{u  I  g,6')  =  lu,=BB-  and  P{6  \  u},g,B')  is  a  Gibbs  distribution  whose  energy  is  a  weighted 
sum  of  a  smoothness  term  and  a  local  distance  between  u  -\-  6  and  g.  The  implemented 
algorithm  alternates  a  deterministic  phase  D  and  a  stochastic  phase  S.  During  the  phase 
D,  we  compute  iteratively  mine  Uo{u!',  0,g).  During  the  phase  S,  given  os',  we  compute  by  a 
Monte  Carlo  method  minj  Ye  ,S,g)P(6  \  u',g,6'),  using  simulations  of  P(^  |  u',g,d'). 

[1]  R.  Azencott,  B.  Chalmond,  F.  Coldefy,  Poc.  ICPR,  The  Hague,  1992. 

[2]  D.  Munford,  J.  Shah,  Proc.  CVPR,  San  francisco  1985. 

[3]  U.  Grenander,  D.  Keenan,  J.  Appl.  Statis.,  1989. 

[4]  F.  Leitner,  1.  Marque,  S.  Lavallee,  P.  Cinquin,  Curves  and  Surfaces,  1991. 

[5]  D.  Terzopoulos,  D.  Metaxas,  IEEE  PAMI,  1991. 

[6]  B.  Chalmond,  Pattern  Recogn.,  1989. 


Dyadic  Affine  Decompositions 
and 

Functional  Wavelet  Transforms 


Charles  K.  Chui 

Texas  A&M  University 
College  Station,  TX  USA 


Decompostion  of  continuous  functions  can  be  accomplished  by  considering  the 
difference  of  consecutive  interpolation  operators. 

In  this  presentation,  we  will  discuss  the  structure  of  such  decomposition  spaces, 
including  the  formulation  of  bases  and  their  duals,  which  lead  to  the  notion  of  the 
functional  wavelet  transform  (FWT).  This  transform  retains  some  of  the  most 
imponant  properties  of  the  integral  wavelet  transform  of  Grossmann  and  Morlet,  such 
as  the  property  of  vanishing  moments  which  has  significant  applications  to  engineering 
problems. 
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Reducing  the  degree  of  Bezier  curves 


M.  A.  COETZEE 

POTCHEFSTROOM  UNIVERSITY  FOR  C.H.E. 


Abstract: 

Various  methods  for  reducing  the  degree  of  Bezier  curves  have  been  inti'oduced,  for  example 
by  Watkins  and  Worsey  in  19SS.  We  found  that  all  of  these  methods  lead  to  discontinuity 
if  you  want  to  reduce  the  degree  of  curves  consisting  of  two  or  more  segments  of  coplanar 
cubic  Bezier  curves.  Watkins  and  Worsey  suggested  that  this  problem  can  be  eliminated 
be  keeping  the  endpoints  of  the  segments  fixed.  This  guarantees  C'°  continuity,  but  it  is 
more  difficult  to  retain  (or  more  spesifically  G^)  continuity.  We  suggest  two  methods 
for  solving  this  problem,  and  discuss  the  generalization  of  these  methods  to  curves  of  higher 
degree. 

References: 

1.  Farin,  G.  Curves  and  surfaces  for  computer-aided  geometric  design,  Academic  Press, 
USA  (1988) 

2.  Gregory,  J.A.  ‘Geometric  continuity’  in:  Mathematical  methods  in  CAGD,  Academic 
Press,  New  York  (1989),  pp  353  -  371 

3.  Watkins,  M.A.  &  Worsey,  A..J.  ‘Degree  reduction  of  Bezier  curves’.  Computer-aided 
design,  Vol  20  No  7  (1988),  pp  398  -  405 
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NEW  METHODS  FOR  2-D  SHAPE-PRESERVING 
INTERPOLATION 

Paolo  Costantini 
Universita'  di  Catania 

It  is  well  known  to  any  researcher  working  in  Numerical  Analysis  or  in  Computer  Aided 
Geometric  Design  that  in  many  practical  applications  it  is  often  necessary  to  produce  curves  or  surfaces 
subject  to  constraints  of  various  kind.  As  a  typical  example,  we  can  consider  the  modeling  of  a  phisical 
phenomena,  described  by  an  (unknown)  increasing  function,  using  a  set  of  experimental  data.  In  this 
case  we  must  use  shape-preserving  interpolation  techniques,  in  which  the  curve  assumes  given  values  at 
the  interpolation  points  and  possesses  the  same  shape  (in  this  case  monotonicity)  of  the  data  . 

The  problems  concerning  shape-preserving  interpolation  (existence,  characterization  and 
algorithmical  construction)  have  been  widely  studied  since  the  seventies  and  a  great  amount  of  papers 
has  been  pubhshed.  However  a  closer  inspection  reveals  that  only  few  papers  are  concerned  with  2-D 
shape-preserving  interpolation  and  many  important  problems  have  not  yet  been  solved  both  for  the 
■’simple”  case  of  data  given  on  a  rectangular  grid  (it  is  not  clear,  for  example,  how  to  obtain  viable 
schemes  for  reproducing  the  convexity  and  concavity  of  the  data)  and  for  the  interesting  one  of  non 
regular  -scattered-  data  (for  which  specifical  methods  have  not  yet  been  devised). 

The  aim  of  this  talk  is  thus  to  present  new  methods  which  can  be  successfully  used  in  shape- 
preserving  2-D  interpolation. 
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Tentative  title: 


Smoothness  of  refinable  functions  obtained  by  convolution 

products 


Stephan  Dahlke 
Wolfgang  Dahmen 
Vera  Latour 

Institut  fur  Geometric  und 
Praktische  Mathematik 
RWTH  Aachen 
Templergraben  55 
W-5100  Aachen 

Abstract 

Based  on  a  given  self-similar  lattice  tiling,  we  construct  a  multi¬ 
resolution  analysis  and  the  associated  wavelet  basis.  A  smooth  gene¬ 
rator  is  obtained  by  the  convolution  product  of  the  characteristic  func¬ 
tions  of  some  specific  self-similar  tUes.  It  turns  out  that  there  is  an 
intimate  connection  between  stationary  subdivision  and  the  concept 
of  self-similar  lattice  tilings.  Employing  some  results  from  stationary 
subdivision  combined  with  a  judicious  choice  of  representers  of  the 
equivalence  classes  induced  by  the  underlying  scaling  matrices  we  are 
able  to  derive  sharp  estimates  for  the  smoothness  of  the  generator. 


SO 


Local  decompositions  of  refinable  spaces  and  some  applications 

Wolfgang  Dahmen 


Abstract:  Multiscale  techniques  have  become  more  and  more  important  in  many 
areas  of  mathematical  theory  and  applications.  They  are  usually  based  on  appropri¬ 
ate  sequences  of  nested  approximation  spaces.  Decomposing  each  such  space  into  a 
direct  sum  of  its  predecessor  and  some  complement  facilitates  successively  updating 
current  approximations  by  including  more  and  more  details  on  finer  scales.  Specifi¬ 
cally,  choosing  orthogonal  complements,  leads  to  wavelet  decompositions.  However, 
in  particular  in  the  multivariate  case,  orthogonality  interferes  with  locality  of  corre¬ 
sponding  bases.  On  the  other  hand,  in  many  applications  it  is  of  crucial  importance 
to  work  with  compactly  supported  basis  functions  while  orthogonality  is  not  so 
essential  and  in  some  cases  not  even  desirable.  This  talk  is  concerned  with  deter¬ 
mining  all  possible  decompositions  of  refinable  spaces  which  are  in  a  certain  sense 
local  and  stable.  Some  of  the  results  to  be  mentioned  were  obtained  in  collaboration 
with  J.M.  Carnicer  and  J.M.  Pena.  It  turns  out  that  these  decompositions  can  be 
conveniently  characterized  in  terms  of  bi-infinite  matrix  relations  in  that  any  such 
decomposition  can  be  generated  from  an  arbitrary  fixed  one  through  certain  matrix 
operations.  In  particular,  this  allows  one  to  identify  decompositions  induced  by 
possibly  nonorthogonal  projectors.  This  has  important  applications,  for  instance, 
in  multilevel  preconditioning  linear  systems  arising  from  Galerkin  discretizations  of 
elliptic  boundary  value  problems  as  well  as  for  matrix  compression  techniques  in 
connection  with  the  numerical  solution  of  integral  or  pseudo  differential  equations. 


Approximation  of  stochastic  processes  with  wavelets 

J.R  D’Ales 

CEREMADE,  Universite  Paris  Dauphine 


Abstract 

The  purpose  of  this  work  consists  in  building  approximations  of  non-stationnary 
random  processes  using  wavelet  bases.  The  case  of  stationnary  processes  has  already 
been  studied  by  A. Cohen,  J.Froment,  and  J.Istas.  We  are  dealing  here  with  transient 
processes,  and  particularly  with  piecewise  stationnary  processes  on  a  time  interval 

[o,r]. 

To  build  the  approximation,  one  performs  a  wavelet  decomposition  of  the  sig¬ 
nal  using  the  well  known  Fast  Wavelet  Transform  algorithm  of  S.Mallat.  One  then 
keeps  the  best  wavelet  coefficients,  and  discards  the  others.  Here  ’’best”  means 
larger  in  amplitude.  Finally  the  approximation  is  obtained  by  applying  the  recon¬ 
struction  algorithm  to  the  shrinked  wavelet  trajisform.  This  confers  the  property  of 
spatial  adaptation  to  our  approximation:  the  asymptotic  performance  is  the  same 
for  stationnary  and  piecewise  stationnary  processes.  That  means  that  the  presence 
of  isolated  discontinuities  in  a  smooth  signal  does  not  alter  the  order  of  the  approx¬ 
imation.  Moreover  we  use  the  orthonormal  wavelet  basis  adapted  on  the  interval 
of  A. Cohen,  I.Daubechies  and  P.Vial.  Thus  the  ’’edge  effects”  are  totally  supressed 
and  the  quality  of  our  approximation  does  not  decrease  near  the  extremities  of  the 
signal. 

The  comparison  has  been  made  with  the  same  type  of  approximations  using 
other  orthonormal  bases  (Fourier,  DCT).  One  shows  that  the  wavelet  approximation 
outperforms  the  Fourier  approximation  in  the  case  of  transient  signals,  but  not  in 
the  case  of  stationnar}'  processes.  Some  numerical  experiments  are  being  processed 
and  should  be  available  at  the  time  of  conference. 
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A  SURFACE-SURFACE  INTERSECTION  ALGORITHM  WITH 
A  FAST  CLIPPING  TECHNIQUE 


Marc  DANIEL,  Alain  NICOLAS 

IRIN  -  Ecole  Centrale  de  Nantes.  1  rue  de  la  Noe  44072  Nantes  cedex  03  -  France 
email :  (daniel  nicolas)(a)cc01.ensm-nantes.fr 

We  present  a  new  approach  for  surface-surface  intersection.  An  hybrid  method  suitable  for  B- 
splines  or  Nurbs  surfaces  has  been  developped.  It  is  decomposed  in  two  main  steps  :  localization 
and  intersections  computations. 

The  localization  is  achieved  through  a  fast  clipping  technique  based  on  the  local  scheme  and  the 
convex  hull  modeling  properties.  The  position  of  a  surface  patch  defined  over  one  parametric 
rectangle  depends  on  the  position  of  a  set  of  neighbouring  control  points.  So,  the  study  of  the 
respective  positions  of  one  surface  control  points  and  the  second  surface  minmax  box  allows  for 
defining  elementary  parametric  rectangles  on  the  first  surface  where  there  are  potential 
intersections.  The  algorithm  is  organized  to  produce  this  result  very  quickly.  Connecting  the 
neighbouring  rectangles  provides,  for  the  first  surface,  free  form  areas  with  potential 
intersections.  Each  area  can  be  bounded  with  a  minmax  box  defined  with  its  associated  control 
points.  With  these  minmax  boxes,  we  can  similarly  compute  the  corresponding  free  form  areas 
for  the  second  surface.  This  clipping  process  is  an  iterative  one  :  it  is  repeated  until  the  different 
areas  on  both  surfaces  remain  unchanged.  The  convergence  is  quickly  obtained. 

Each  set  of  corresponding  areas  is  independently  proceeded.  It  is  first  checked  if  there  is  a 
degeneracy  (tangency  or  partially  identical  surfaces)  or  a  closed  curve.  This  testi^  studies  and 
compares  normal  surface  bounding  cones.  It  has  been  extended  for  degeneracies  detection  ;  the 
difference  being  the  existence  of  an  intersection  point  with  colinear  normal  vectors  in  the  case  of 
degeneracies.  If  a  closed  curve  is  detected  (degeneracies  are  separatly  computed  using  a  net  of 
isoparametric  curves  on  both  surface  and  a  robust  curve-curve  intersection  algorithm),  the  area  of 
each  surface  is  split  into  four  new  areas  with  the  aim  to  only  obtain  opened  curve  segments  on 
each  part.  This  splitting  involves  surface  subdivisions  only  if  the  area  is  defined  over  one 
elementary  parametric  rectangle.  Otherwise,  the  new  areas  are  immediatly  deduced  from  the 
separation  of  the  parametric  domain.  The  process  for  closed  curves  elimination  is  a  recursive 
one.  Finally,  only  opened  intersection  curves  can  exist  on  each  couple  of  areas.  They  necessary 
begin  on  area  boundary  curves  (isoparametric  curves)  and  are  computed  with  a  marching 
method.  The  opened  curve  segments  are  eventually  connected  together  to  produce  the  existing 
closed  curves. 

(1)  T.W  Sederberg,  H.  N.  Christiansen,  S.  Katz,  “Improved  test  for  closed  loops  in  surface 
intersections”.  Computer  Aided  Design,  vol.  21,  n°  8,  pp.  505-508,  1989 
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Curves  and  Surfaces,  Chamonix  1 993. 


Zonoidal  surfaces^. 


O.  Daoudi,  B.  Lacolle,  N.  Szafran, 

Universite  Joseph  Fourier,  Laboratoire  LMC-IMAG,  BP  53  X,  F-38041  Grenoble  Cedex,  F ranee 

and 

P.  Valentin 

ELF-FRANCE,  Centre  de  Recherche  ELF  SOLAIZE,  BP  22,  69360  S‘  Symphorien  D'Ozon,  France. 


In  the  euclidian  space  £  =  IR”,  a  zonotope  is  a  convex  centrally  symmetric  set,  defined  as 
the  Minkowski  sum  of  segments.  The  closure,  in  the  Hausdorff  metric,  of  the  set  of 
zonotopes  in  £,  is  a  set  of  convex  bodies  so  called  zonoids.  Some  papers  deal  with  the 
general  theory  of  zonoids  and  approximation  by  sequences  of  zonotopes.f^ 

The  zonoids  arise  naturally  in  many  contexts.  The  motivation  of  our  study  is  the 
geometric  management  of  mixtures  which  issues  from  actual  applications  in  oil  industry 

In  the  geometric  model  we  use,  any  product  can  be  viewed  as  a  vector  in  IR^^  one  component 
for  each  species,  and  the  mixing  procedure  corresponds  to  vector  addition  in  the  vector  space 

[R”.  Given  a  set  of  basic  products  n^,n2,...,np,  a  product  m  is  feasible  if  and  only  if  m 

lies  in  the  set  Zinun2 . «„)  =  (Y  .  0<A^<1  il<j<p)]  which  is  a  zonotope. 

In  this  paper,  we  consider  the  three  dimensional  euclidian  space  and  we  deal  with 
particular  zonoids  defined  by  : 

Z(/)  =  j/n  =  1(0 /(r)  dt,  1  e  L^(/),  0<1(0<1  (re/)], 

with  /  a  compact  interval  in  !R  and/ a  function  in  °[/,IR^].  For  example,  a  zonoid  in  this 
previous  form  can  be  viewed  as  the  set  of  feasible  mixtures  by  a  distillation  process. 

We  name  the  boundary  of  a  zonoid  as  a  "zonoidal  surface"  and  we  study  some  structural 
properties  of  this  surface.  Using  a  generalization  of  the  projective  diagram  of  a  zonotope,  we 
obtain  the  decomposition  of  the  zonoidal  surface  into  a  set  of  convex  regular  patches.  In  the 
particular  case  of  a  "convex"  projective  diagram,  the  boundary  of  the  zonoid  is 
decomposable  into  one  or  two  patches  and  we  obtain  simple  parametric  equations  of  these 
patches.  This  basic  case  is  used  in  a  general  framework  leading  to  the  parametric  equations 
of  the  patches  for  a  large  class  of  zonoidal  surfaces.  Then,  we  can  show  how  the  regularity 
of  these  patches  is  related  to  the  regularity  properties  of  the  function/.  At  last,  we  present 
numerical  methods  for  visualizing  a  zonoid,  computing  plane  sections  of  a  zonoid,  and  other 
related  problems. 
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An  exclusion-bissection  algorithm 
for  computing  the  intersection  of  curves  and  surfaces 

Jean-Pierre  DEDIEU 

Universite  Paul  Sabatier 
Toulouse,  France 
Email:  dedieu@cict.fr 


The  exclusion-bissection  algorithm  describded  in  this  talk  computes  the  intersection  of 
n  algebraic  hypersurfaces  in  R"  with  an  accuracy  r  in  0(Log(l/r))  steps. 


Fair  Surface  Design  using 
Recursive  Subdivision  Surfaces 


Tony  DeRose 
University  of  Washington 


There  has  been  considerable  recent  activity  in  the  area  of  fair  surface  design. 
Although  several  new  schemes  offer  impressive  improvements,  they  still  suffer  from 
extreme  computation  time,  or  are  limited  in  the  topological  types  of  the  surfaces  they 
can  represent. 

In  this  talk  I'll  argue  that  recursive  subdivision  surfaces,  such  as  the  ones 
developed  by  Doo-Sabin  and  Catmull-Clark,  offer  potential  advantages  when 
attempting  to  overcome  the  above  deficiencies.  In  particular.  I'll  describe  techniques  for 
efficiently  computing  fair  subdivision  surfaces  that  interpolate  the  vertices  of  a  given 
polyhedron. 

Joint  work  with:  Michael  Kass  and  Mark  Halstead  of  Apple  Computer,  Joe 
Warren  of  Rice  University,  and  Michael  Lounsbery  of  the  University  of 
Washington. 
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Abstract 

Image  Processing  and  Hyperbolic  Wavelets 
Ronald  A.  DeVore,  University  of  South  Carolina 


Wavelets  have  had  successful  application  to  several  areas  including  image  processing 
and  surface  compression.  We  shall  review  some  of  these  applications  and  point  out  some 
deficiencies  of  the  usual  wavelets.  We  will  then  introduce  hyperbolic  wavelets,  discuss  some 
of  their  properties,  and  apply  them  to  generic  problems  in  image  processing. 
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Z.  Ditzian 


Title  of  talk 

"Multivariate  Bernstein  polynomial  approximation" 


Abstract 

Relation  between  Bernstein  polynomial  approximation  and  a 
K— functional  defined  with  the  aid  of  a  weakly  elliptic  operator,  is 
established. 
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Spline  Conversion,  Existing  Solutions,  and  Open  Problems 

Tor  Dokken  and  Tom  Lyche 


SINTEF  SI  and  The  University  of  Oslo 
January  29,  1993 


Abstract 

There  is  a  great  need  in  industry  for  methods  to  convert  between  different  representation 
formats.  Such  conversion  involves  much  more  than  transformation  between  different  spline 
formats.  Consider  for  example  the  geometry  received  from  PC  based  CAD-systems.  This 
is  often  represented  by  chains  of  straight  lines  representing  a  more  complex  shape.  Another 
problem  is  conversion  of  trimmed  surfaces.  The  transformed  surface  is  often  not  usable  for 
further  design.  Examples  like  these  show  that  the  standard  conversion  tools  in  use  today 
are  not  sufficient  to  fulfill  industrial  needs.  In  our  talk  we  will  give  an  overview  of  some  of 
the  conversion  problems  facing  industry,  describe  alternative  solutions,  and  list  some  open 
problems  still  existing. 
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SUBDIVISION  SCHEMES 


Nira  Dyn 

School  of  Mathematical  Sciences 
Tel  Aviv  University 
Tel  Aviv,  Israel 


In  this  talk  various  classes  of  subdivision  schemes  for  the  generation  of  curves  and  surfaces 
will  be  reviewed:  the  uniform  stationary  schemes  and  their  extensions  to  non-stationary  and  non- 
uniform  schemes  and  to  integral  subdivision  schemes.  Typical  examples  of  these  classes  will  be 
presented.  Conditions  for  strong  and  weak  convergence  of  subdivision  scheraes-will  be  discussed. 

The  relation  to  multi-resolution  analysis  and  to  wavelets  theory  with  the  analysis  of  a  specific 


example  of  C”  wavelets  will  be  reviewed. 


A  Stepwise  Algorithm  for  Converting  B-Splines 


Matthias  Eck 

T.H.  Darmstadt,  Germany 


In  this  talk  a  method  is  proposed  for  converting  B-spline  curves  of  order  k  (resp. 
B-spline  surfaces  of  order  {k,  1))  approximatively  into  B-spline  curves  of  any  lower 
order  k  (resp.  B-spline  surfaces  of  any  lower  order  (fc,  /)).  Further,  it  is  required  for 
practical  purposes  mainly  that  the  resulting  B-splines  do  not  consist  of  too  many 
spline  segments. 

In  contrast  to  most  other  conversion  methods  described  in  the  literature  and  treating 
the  above  polynomial  problem,  omr  algorithm  carries  out  the  conversion  stepwise. 
Here,  in  each  step  either  the  order  of  one  particular  segment  is  lowered  by  one  or 
two  neighbouring  segments  eire  merged  together  to  only  one  segment.  Then  these 
u  vvo  major  steps  are  repeated  in  a  specific  order  until  the  target  order  of  the  B-spline 
is  reached.  Further,  one  can  hope  that  the  final  number  of  segments  is  low  as  the 
method  works  carefuUy  in  this  regard. 

Both  major  steps  are  based  on  recent  restilts  of  the  author.  So  the  first  step  is  nearly 
identical  with  a  degree  reduction  scheme  for  polynomial  Bezier  curves  and  surfaces 
whereby  certain  continuity  conditions  at  the  boundaries  are  preserved.  Further,  the 
second  step  is  very  close  to  a  knot  removal  scheme  for  polynomial  B-splines. 
Summarized,  the  main  advantages  of  these  two  methods  are: 

1.  The  new  control  points  are  obtained  by  simple  and  explicitly  known  construc¬ 
tions  from  the  given  control  points  only.  Here  it  turns  out  that  the  bivariate 
case  is  just  a  two-fold  application  of  the  curve  construction,  firstly,  to  all  rows 
and,  secondly,  to  all  columns  of  the  control  net.  Altogether,  these  calculations 
are  very  fast. 

2.  Upper  bounds  of  the  maximal  approximation  error  for  all  t  are  apriori  known 
in  each  step.  These  bounds  can  be  used  to  carry  out  the  conversion  within  a 
given  pointwise  tolerance.  Here,  of  course,  we  have  to  subdivide  spline  segments 
where  the  error  tolerance  is  exceeded. 


Knot  Removal  for  B- Spline  Curves 


Matthias  Eck  &;  Jan  Hadenfeld 


T.H.  Darmstadt,  Germany 


The  approximative  problem  of  knot  removal  is  well  known  and  has  been  discussed 
by  several  authors  before  (e.g.  Lyche,  M0rken  ’90).  Here,  it  is  seen  as  a  procedure 
for  representing  a  given  B-spline  curve  X  of  order  k,  defined  over  the  knot  vector 
T  =  (to,  ti, .  - . ,  tn+Jt),  as  good  as  possible  by  another  B-spline  curve  X  of  same  order 
k  but  defined  over  a  shorter  knot  vector  T  =  (to,  ti,  • .  •  ,tTn+fc),  being  a  subset  of  the 
original  knot  vector  T. 

In  this  talk  we  mainly  consider  the  special  case  m  =  n  —  1  which  means  the  removal 
of  only  one  knot.  This  restriction  is  necessary  here  since  we  are  interested  in  explicit 
solutions  for  the  new  curve  X. 

Doing  so,  we  interpret  knot  removal  as  an  inverse  procedure  of  knot  insertion.  This 
idea  leads  us  to  a  very  simple  and  local  construction  of  the  new  control  points  from 
the  ones  of  the  given  curve.  Now,  this  construction  contains  some  additional  degrees 
of  freedom  which  are  determined  explicitly  in  several  ways  for  a  general  knot  vector 
T  as  well  as  for  general  order  k  by  minimizing  different  error  functions  between  X 
and  X. 

The  error  functions,  we  used,  are  either  ’’continuous”  norms,  considering  the  differ¬ 
ence  for  all  parameter  values,  or  ” discrete”  norms,  considering  only  the  difference 
between  the  old  and  new  control  polygon. 

All  methods  are  illustrated  with  several  examples  which  also  show  that  good  re¬ 
sults  can  be  obtained  if  we  remove  more  than  one  knot  by  applying  our  constuctive 
algorithm  for  one  knot  after  another. 
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Shape  from  Motion 


Michael  Eichhard  Knut  Moller 

University  of  Bonn 
Department  of  Computer  Science 
Romerstr.  164 
D-5300  Bonn,  FRG 

January  27,  1993 


Abstract 


The  need  to  naodel  surfaces  is  of  great  importance  in  Computer  Aided  Design  (CAD) 
and  in  a  different  setting  in  control  applications.  While  CAD  research  is  dominated  by 
searching  for  efficient  representations  and  manipulations  of  surfaces,  control  is  concerned 
with  modelling  dynamic  effects.  Often  in  optimal  control  a  prediction  problem  needs  to 
be  solved  (e.g.  grasping  a  rolling  ball  with  a  robotmanipulator  relies  on  a  prediction  of 
the  ball  trajectory,  knowledge  about  the  manipulator  dynamics  etc.).  In  control  surface 
representations  are  mainly  needed  for 


1.  environmental  modelling  (prediction  purposes)  and 

2.  as  so  called  decision  planes  supporting  fast  decision  making  in  real-time  control 
systems. 


Usually  the  construction  of  surface  modells  in  CAD  is  based  on  Bezier-representations 
which  provide  a  powerful  tool  for  approximating  curves  and  surfaces.  These  represen¬ 
tations  have  a  number  of  attractive  features  e.g.  quasi-localit)^,  differentiability,  com.- 
pactness  etc.  But  (from  the  control  point  of  view)  unfortunately  the  CAD-designer  is 
interactively  involved  in  manipulating  representations  until  they  fit  his  needs.  In  learn¬ 
ing  control  a  different  problem  has  to  be  faced.  From  possibly  noisy  observations  (in  our 
case  of  a  moving  object)  an  efficient  representation  of  a  surface  needs  to  be  derived.  We 
were  particularly  interested  how  representations  can  be  utilized  for  control  applications. 

In  the  talk  we  will  give  a  specification  of  the  problem  followed  by  a  description  of  how 
to  include  dynamics  into  the  process  of  generating  a  surface  representation.  Finally  the 
approach  is  compared  with  a  number  of  different  other  approximation  schemes  including 
fashionable  “neural”  learning  methods. 


Best  constrained  approximations  of 
planar  curves  by  Bezier  curves 


Eberhard  F.  Eisele 


In  CAD  systems  methods  for  good  approximation  of  a  given  curve  by  a  certain  class  of 
spline  curves  are  frequently  needed.  For  instance,  a  method  for  approximate  conversion 
of  spline  representations  is  required  for  exchanging  data  between  different  CAD  systems 
or  the  problem  arises  to  approximate  offset  curves  by  spline  curves  of  a  certain  kind. 

Given  a  parametric  planar  curve  p  and  any  Bezier  curve  x  of  degree  n  such  that 
p  and  X  have  contact  of  fixed  order  k  at  the  common  end  points,  we  use  the  normal 
vector  field  of  p  to  measure  the  distance  of  corresponding  points  of  p  and  x.  Then  the 
Bezier  curve  x  is  called  a  best  approximation  to  p  if  the  maximum  norm  ||px  W^c 
distance  (or  deviation)  function  px  is  minimal  for  x  and  the  best  approximation  x  has 
been  characterized  by  an  alternation  property  of  px  in  [Degen  ’92]  and  [Eisele  ’92]. 

This  paper  addresses  to  the  problem  of  characterizing  a  best  approximation  x  sa¬ 
tisfying  certain  constraints.  We  consider  the  cases  that  either  each  x  interpolates  to  p 
at  certain  prescribed  points,  or  that  each  deviation  function  px  has  a  restricted  range;  in 
particular  the  second  case  includes  one-sided  Chebyshev  approximation.  In  both  cases 
we  characterize  a  best  constrained  approximation  x  by  an  alternation  property  of  pxi 
applying  results  of  Chebyshev  approximation  with  side  conditions.  Finally,  Remes  type 
algorithms  can  be  used  for  the  computation  of  best  constrained  approximations. 
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A  Geometry  Based  Approach  to  the  Generation 
of  Unstructured  Surface  Grids 


Stefan  Farestam 

CE?Gm 

European  Centre  for  Research  and  Advanced  Training  in  Scientific  Computing 
42  avenue  Gustave  Coriolis  •  31057  Toulouse  •  France 

We  introduce  a  method  to  generate  unstructured  triangular  surface  grids,  that  uses  the  intrinsic  geometry 
of  the  surface  to  define  the  size  and  the  shape  of  the  triangular  elements.  The  method  is  based  on  the 
Advancing  Front  Technique,  but  instead  of  using  a  background  grid,  a  transformation  tensor,  T,  is  defined 
using  the  principal  curvatures,  and  the  principal  directions  of  curvature.  These  properties  are  intrinsic 
to  the  surface,  and  not  dependent  on  the  particular  parametrisation  used,  the  resulting  grid  is  also 
independent  of  the  parametrisation. 

Given  a  parametrically  defined  surface,  the  two  principal  curvatures  and  the  two  principal  di¬ 

rections  of  curvature  are  defined  by  two  generalised  eigenvalues,  /ci  2,  and  unit  eigenvectors,  €i_2,  of  the 
system; 

De  =  KGe, 

where  G  and  D  are  the  first  and  second  fundamental  matrices  of  the  surface,  respectively.  The  target 
triangle  to  be  generated  at  x  may  now  be  defined  as  a  triangle  inscribed  in  an  ellipse  centered  at  x,  where 
the  two  axes  of  the  ellipse  are  aligned  to  ei,2,  and  their  length  is  proportional  to  k;i,2-  We  refer  to  the 
constant  of  proportionality,  8,  as  the  geometry  resolution,  and  it  is  globally  defined  for  the  surface. 

The  transformation  tensor,  T,  defines  a  mapping  from  the  parameter  space  of  the  surface  to  a  control 
space,  7^2^  in  which  the  triangulation  algorithm  (AFT)  strives  to  generate  equilateral  triangles  of  unit 
size.  Given  the  above,  we  can  now  define  the  tensor  T  as; 

~^2j  Cl  Ml  €ij  Cj  _  ’ 


where 

C  =  ||ei  X  €2!!,  Mi^2  =  5i,2/  Si, 2  =  ^1,2- 

An  equilateral  triangle  of  unit  size  in  will  thus  have  the  desired  size  and  shape  on  the  surface. 


A  sample  grid  on  the  surface  x(^,  77)  =  rj,  sin^  sin  77). 


Projective  cubic  B-splines 

Gerald  Farin 
Computer  Science 
Arizona  State  University 
Tempe,  AZ85287 
USA 

farin(2)  asuvax.asu.edu 


Abstract:  We  will  consider  piecewise  cubics  in  a  projective 
environment.  The  concept  of  differentiability  will  be  based  on  that  of 
osculants,  which  are  certain  blossoms. 

We  will  show  how  to  construct  B -spline-like  curves  with  these  tools, 
and  show  how  they  relate  to  standard  rational  B-splines. 

The  talk  will  include  a  projective  treatment  of  conics  in  the 
Bemstein-Bezier  form. 


ALGORITHMS  INTENDED  TO  FIND  A  SIMPLICIAL  SURFACE 
FROM  A  SET  OF  3D  DATA  POINTS 

Chantal  FAVARDIN 

This  study  consists  in  elaborating  algorithms  based  on  geometric  criteria,  that  allow  to  reconstruct 
a  connected  surface  S,  from  a  discrete  unorganized  representation.  More  precisely,  we  want  to  determine, 
from  the  only  knowledge  of  the  coordinates  of  n  IR^-points  iV//’'  =  i  —  l...n  belonging  to  <S,  a 

polyhedral  surface  with  triangular  faces  (named  simplicial  surface),  whose  vertices  are  the  points  Mp,  and 
that  uniformly  converges  to  the  initial  surface,  when  the  numbers  of  data  points  n,  tends  to  infinity.  Many 
practical  applications,  such  as  pattern  recognition,  computer  vision,  medical  imagery,  where  the  objects  are 
represented  by  a  set  of  data  points  only  known  by  their  coordinates,  aroused  interest  in  this  type  of  study, 
that  also  can  constitute  a  preliminary  stage  for  data  interpolation  problems.  Two  types  of  approach  were 
considered  up  to  now.  The  first  one  assumes  that  the  data  are  divided  up  into  planar  crossed  sections  and 
so  have  an  initial  structure  ([2],  [4]).  These  methods  consist  in  determining  a  triangulation  between  two 
adjacent  crossed  sections.  The  second  one  assumes  that  the  points  have  no  initial  structure  ([l].[3],'5])  and 
that  only  the  coordinates  of  the  points  are  known.  This  is  the  approach  that  we  consider  here.  We  present 
two  new  algorithms  to  solve  this  problem. 

The  first  one  consists  in  maximizing  the  dihedral  angle  between  two  adjacent  faces,  and  concerns 
closed  convex  surfaces.  More  precisely,  let  us  assume  that  the  algorithm  has  already  built  a  simplicial 
surface,  whose  vertices,  taken  among  the  points  are  noted  We  denote  by  P^^,....P^^ 

the  vertices  belonging  to  the  boundary  of  this  simplicial  surface,  (two  points  with  consecutive  suffixes  are 
connected).  To  extend  this  polyhedral  surface,  we  determine  for  each  j  G  {1, (}  a  point  P^  defined  by. 

.  6^,  =  ({Mr,  1  <  t  <  n}\{Pf  ii  <  i  <  ^})  u 

•  P’^  ^  k  maximizes  the  convex  dihedral  angle  between  the  triangles  {^P^.P^.  .  a-^id 

iPh  Pr  M-"^)  for  all  M^  G  Uk  (the  point  designates  the  point  with  which  the  points  P^  and  P^ 
form  a  triangle  of  the  simplicial  surface  built  so  far. 

The  algorithm  stops  when  the  set  of  points  belonging  to  the  boundary  of  the  polyhedral  surface  is  empty. 

We  prove  that  the  criterion  enables  to  obtain  a  triangulation  of  the  convex  hull  of  the  points  M^ .  a^nd 
that  this  triangulation  is  the  only  one  if  four  points  (or  more)  coplanar,  are  never  cocyclic.  So  this  criterion 
fits  quite  the  case  of  points  belonging  to  a  closed  convex  surface. 

The  second  algorithm  enables  to  adapt  this  criterion  to  the  case  of  points  belonging  to  a  closed  non- 
convex  surface.  The  point  P^  is  chosen  in  a  neighborhood  of  P^  and  P^  ^ .  An  insertion  stage  is  also 
added.  The  case  of  surface  with  boundary  curves  is  solved  by  an  heuristic.  The  expected  time  complexity  of 
these  different  algorithms,  that  are  illustrated  by  various  examples,  is  in  O(Tllogn). 

1.  J.D.  Boissonnat, Geometric  structures  for  three-dimensional  shape  representation,  ACM  Trans,  on 

Graph..  Vol.3,  No.4,  (Oct.  1984),  pp.266-286. 
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Processing,  Vol  44,  No  1.  (Oct  1988),  pp.1-29. 
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July  1992.  SIGGRAPH’92  Conference  Proceedings. 

4.  E.  Keppel,  Approximating  complex  surfaces  by  triangulation  of  contour  lines,  IBM  J.  Res.  Develop. 

(Jan.  1975),  pp.2-10. 

5.  J. O’Rourke,  Polyedra  of  minimal  Area  as  3D  object  Models,  in  Proceeding  of  the  international  joint 

conference  on  artificial  intelligence  81  (Vancouver,  Canada)  1981,  pp. 664-666. 
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Interpolating  by  conic  splines 
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*  Universite  de  Valenciennes, 
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BP  311,  F-  59304  Valenciennes  Cedex 
**  Universite  des  Sciences  et  Technologies  de  Lille, 
UFR  de  Calais,  Laboratoire  ANO,  Batiment  M3, 
F  -  59655  Villeneuve  d’Ascq  -  Cedex 


Polynomial  spline  curves  do  not  allow  to  interpolate  given  points  of  the 
plane  and  their  derivatives.  Conic  splines  i.e.  piecewise  conic  curves  solve 
this  Hermite  problem.  This  interpolating  problem  can  consider  points  at 
infinity.  Conditions  are  given  so  that  each  conic  segment  does  not  include 
points  at  infinity  except  possibly  at  its  interpolating  end-points. 
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PLANAR  RATIONAL  CUBICS  AS  BR-CUBICS 
J.C.  Fiorot,  P.  Jeannin  and  S.Taleb 


The  concept  of  massic  vectors  enables  to  represent  all  rational  parametric  curves 
under  B-rational  form,  BR-form  in  abrev:  [Fiorot  k  Jeannin ’89, ’91, ’92].  A  BR-curve 
is  determined  by  a  massic  polygon,  likewise  a  Bezier  curve  is  determined  by  a  control 
polygon.  As  particular  cases,  the  BR-form  includes  polynomial  Bezier  curves,  and 
rational  Bezier  curves  [Farin’90]. 

In  BR-form,  a  rational  cubic  is  determined  via  Bernstein  polynomials  by  four 
massic  vectors. 

In  this  work,  we  treat  the  problem  of  detecting  some  geometric  elements  such  as  - 
cusps,  loops,  inflection  points,  asymptotes  and  parabolic  points  on  the  BR-cubic. 

Conditions  for  the  occurence  of  cusps,  loops,  and  inflection  points  are  investi¬ 
gated.  This  leads  to  a  projective  classification  of  planar  rational  cubics.  This  study 
gives  a  new  approach  and  another  numerical  point  of  view  than  the  algebraic  aspect 
given  in  [Patterson’88],  [Walker’50],  in  order  to  use  them  in  computer  science,  CAGD 
and  CAM. 

We  also  determine  conditions  for  asymptotes  and  parabolic  points  in  terms  of 
massic  vectors. 
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A  GENERALIZATION  OF  BEZIER  CURVES  AND  SURFACES 


I.  GANSCA*,  Gh.  COMAN  and  L.  jAMBULEA’* 


Abstract:  As  it  is  known  a  Bezier  curve  (or  surface)  is  defined  by  using  the  Bernstein  basis 
A:=d^,  ?e[0,l]  and  a  set  of  points  Pj^,  k=Q,n, 

(or  P.J,  i=^,  j=0,n,  for  a  surface) 

In  this  paper  we  change  the  Bernstein  basis  by  the  follov^ng  basis  (D.D.  Stancu) 

■(”^^)  0^k<r 

n- 

k=Q^,  fe[0,l], 

where  r  is  a  non-negative  integer  parameter  conditioned  by  the  inegality  2r  <  n,  and  obtain  curves 
and  surfaces  that  generalize  the  Bezier  curves  and  surfaces,  because,  for  r=0  and  r=l  results 
Wn,r,k(t)=bn,k(t)»  Vte[0,l].  Taking  in  view  that  te[0,l],  k  =  0,n  and  ^ 

new  partition  of  the  interval  [0,1]. 

One  shows  that  the  well  known  de  Casteljau  algorithm  is  applicable,  but,  in  this  case,  for 
the  variable  points  Q,^(t)  =  (l-t)P[(.+tPj.+k  k=0,n-r,  t6[0,l].  For  r=0  results  Qj.(t)=Pij  and  for  r=l 
we  obtain  Qk(t)  =  (l-t)Pk+tPk+n  k=0M,  that  is,  the  points  Q^t)  rezult  from  the  given  points  P^, 
k=dA  after  the  first  step  of  the  Casteljau  algorithm. 

Also,  the  paper  contains  formulas  for  derivatives  and  some  figures  with  curves  and  surfaces 
which  are  compared  with  the  standard  Bezier  curves  and  surfaces  corresponding  to  the  same  set 
of  given  points. 


'')  t\l-ty-^-^*^+  r^k^n-r 

p  n-r<k^n. 
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Comer  cutting  algorithms  and  totally  positive  matrices. 

Mariano  Gasca  and  Juan  Manuel  Pena. 

Depto.  Matematica  Aplicada,  University  of  Zaragoza,  Spain 

Comer  cutting  algorithms  play  an  important  role  in  Computer 
Aided  Geometric  Design.  Many  authors  have  shown  their  conection  with 
Total  Positivity.  In  [3]  Goodman  and  Micchelli  defined  precisely  these 
algorithms  and  proved  that  there  exists  a  corner  cutting  algorithm 
transforming  a  control  polygon  into  another  one  with  the  same  number 
of  vertices  if  and  only  if  the  matrix  which  relates  both  polygons  is 
totally  positive  and  stochastic.  In  fact  they  gave  a  detailed  matricial 
description  of  such  algorithms.  Goodman  and  Said  [4]  proved  that  there 
are  corner  cutting  algorithms  which  transform  the  control  polygon 
associated  to  a  Ball  basis  of  polynomials  of  a  curve  y  into  the  Bezier 
control  polygon  of  y.  They  conjectured  that  the  same  could  happen  for 
the  control  polygon  associated  to  any  totally  positive  basis  of 
polynomials.  The  conjecture  was  recently  proved  by  Camicer  and  Pena 
[1]. 

In  this  paper,  we  use  the  results  to  appear  in  [2]  and  prove  that 
under  an  extra  (but  quite  natural)  condition  to  be  added  to  the  ones  in 

[2]  the  comer  cutting  algorithms  which  are  studied  there  can  be  carried 
out  in  a  unique  way.  More  precisely,  for  a  given  nonsingular  totally 
positive  stochastic  matrix  there  exists  a  unique  factorization  which 
describes  a  comer  cutting  algorithm. 

[1]  Camicer,  J.M.  and  J.M.  Pena :  Shape  preserving  representations 
and  optimahty  of  the  Bernstein  basis.  To  appear  in  Advances  in 
Computational  Mathematics  (1992). 

[2]  Gasca,  M.  and  J.M.  Pena :  A  matricial  description  of  Neville 
elimination,  with  applications  to  total  positivity.  To  appear  in  Linear 
Alg.  and  Apphc. 

[3]  Goodman,  T.N.T.,  Micchelli  C.A.:  Comer  cutting  algorithms  for  the 
Bezier  representation  of  free  form  curves.  Linear  Alg.  AppL  99  (1988). 

[4]  Goodman,  T.N.T.,  Said,  H.B.:  Shape  preserving  properties  of  the 
generalized  BaU  basis.  Computer  Aided  Geometric  Design  8  (1991). 
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It  is  well-known  that  generally  the  offsets  of  a  rational  parametric  curve 
are  not  rational.  We  denote  by  1%  the  set  of  the  rational  curves  with  rational 
offsets. We  give  a  geometrical  characterization  of  the  curves  of  /T^and  the 
explicit  form  of  their  parametrization.The  arc  length  of  a  curve  of  In  is  an 
integral  of  a  rational  function. We  specify  .the  subset  An  of  In  consisting  of 
the  rational  curves  with  rational  arc  length:We  show  that  it  is  the  set  of  the 
caustics  of  the  rational  parametric  curves. 


Composite  Curve  Estimation  for  Human  Body  Modelling 


Alexander  M.  Geurtz 


Riccardo  Leonard!^ 


Murat  Kunt 


Signal  Processing  Laboratory 
Swiss  Federal  Institute  of  Technology 
EPFL-Ecublens,  CH-1015  Lausanne,  Switzerland 

We  will  report  on  the  state  of  our  research  on  the  analysis  of  human  movement  from  digital  image 
sequences.  In  this  work,  the  human  body  is  modelled  as  a  set  of  articulated  objects  that  are 
considered  representative  for  the  overall  shape  and  the  task  is  to  estimate  the  parameters  of  these 
objects  from  features  detected  in  the  image  sequence.  For  the  sake  of  mathematical  and  numerical 
tractability  the  model  is  restricted  to  two-dimensional  planar  ellipses  and  applied  to  static  images 
(although  intended  for  a  future  extension  to  sequences).  As  feature  data  we  use  the  (sparse)  contour 
•Doints  detected  in  the  grey  level  images.  Thus,  a  set  of  ellipses  is  fitted  to  the  contour  points.  As 
the  error  measure  for  the  fitting  of  one  curve  we  use  the  sum  of  squares  of  the  deviations  e,-  = 
/(s,;p)  of  the  implicit  relation  for  the  ellipse  at  contour  points  x,-,  with  p  representing  the  position, 
orientation,  and  size  of  the  curve.  The  segmentation  problem  of  identifying  which  contour  point 
belongs  to  which  ellipse,  a  major  problem  when  fitting  multiple  curves  to  a  data  set,  is  integrated 
in  the  optimization  task  by  applying  a  so-called  fuzzy  representation.  In  this  representation,  a 
normalized  membership  vector  is  associated  with  each  feature  point,  with  the  components  of  this 
vector  indicating  the  importance  (weight)  of  the  relation  between  the  data  point  and  each  of  the 
curves.  The  error  measure  for  a  set  of  ellipses  is  the  sum  of  squares  of  deviations,  weighted  by 
the  membership  coefficients.  The  spatial  relations  between  the  individual  ellipses  are  expressed  as 
nonlinear  constraints  on  the  resulting  optimization  problem.  A  hierarchical  decomposition  allows 
to  choose  the  level  of  detail  of  the  representation  and  helps  reducing  the  risk  of  being  trapped  in 
one  of  the  numerous  local  minima  of  the  nonlinearly- constrained  optimization  problem.  We  will 
report  results  on  both  artificial  and  real  data  sets. 
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Invariant  Approximation  of  planar  star-shaped  contours 


F.  Ghorbel  and  V.  Burdin 

Abstract: 

Euclidean  invariance  is  widely  recognized  as  an  important  property  of 
curve  and  surface  shape  description  algorithms  for  computer-aided  design, 
pattern  recognition  and  image  analysis. 

The  Euclidean  transformations  are  defined  by  combining  rotations  and 
translations.  Euclidean  invariance  is  equivalent  to  axis-independence.  This 
simply  means  that  applying  the  curve  algorithm  to  the  transformed  data  set  is 
equivalent  to  applying  the  algorithm  to  the  original  data  set  and  transforming 
the  resultant  curve  or  surface. 

In  a  series  of  papers,  authors  have  investigated  the  general  functional  form 
of  curve  algorithms  that  are  invariant  under  certain  transformation  groups.  In 
this  paper,  we  use  the  radial  curve  representation  with  arc  length  parameter- 
isation  in  order  to  obtain  axis-independence  for  any  planar  curve  algorithms. 
However,  we  have  a  restriction,  we  only  consider  a  particular  class  of  planar 
shapes  called  star-shaped  curves. 

We  present  a  study  of  the  polar  arc  length  parameterization  and  we 
show  that  the  radial  function  which  is  invariant  with  respect  to  translation 
and  rotation,  is  sufficient  to  represent  completely  the  star-shaped  curves.  We 
obtain  an  invariant  curve  algorithm  by  approximating  this  radial  function 
rather  than  cartesian  coordinates. 

In  many  applications  in  pattern  recognition  or  image  analysis,  this  pa- 
rameterisation  can  be  justified  by  its  adaptation  to  a  wde  range  of  shapes  (for 
example  in  medical  and  biological  imaging).  Finally,  we  present  the  advan¬ 
tages  and  the  limits  of  this  method  by  giving  results  concerning  a  particular 
class  of  planar  contours. 


Polynomial  Approximation  of  Spheres 

Sigurd  Glaerum 
University  of  Oslo 

Abstract.  A  quadrilateral  approximation  scheme  for  spherical  surfaces  is  pre¬ 
sented.  The  scheme  is  based  on  a  simple  generalization  of  circle  approximations. 
It  will  be  shown  that  the  spherical  aproximant  retain  many  of  the  properties 
of  the  circle  approximations.  For  instance,  if  the  circle  approximation  is  of  or¬ 
der  k,  the  spherical  approximation  will  be  of  order  k  in  each  parameter.  The 
geometric  continuity  between  two  neighboring  sphericed  approximations  is  also 
closely  linked  to  the  geometric  continuity  between  two  neighboring  circle  ap¬ 
proximations.  The  problem  of  dealing  with  the  singularities  that  occur  as  a 
result  of  using  quadrilateral  patches  is  briefly  considered. 


Dual  Polynomial  Bases 


Ronald  N.  Goldman 
Department  of  Computer  Science 
Rice  University 
RO.  Box  1892 
Houston,  Texas  77251 


Abstract 

Blossoming,  the  de  Boor-Fix  dual  functionals,  Marsden’s  identity,  the  Oslo 
algorithm,  recursive  procedures  for  evaluation  and  differentiation:  these  are 
topics  usually  associated  only  with  B-splines.  But  all  these  formulas  have  local 
polynomial  interpretations.  The  goal  of  this  talk  is  to  extend  these  formulas  and 
techniques  to  totally  arbitrary  polynomial,  and  locally  linearly  independent 
piecewise  polynomial,  bases.  These  extensions  help  to  unify  the  theory  of 
univariate  polynomials  and  splines,  and  they  also  provide  some  additional 
perspective  on  the  special  status  of  the  B-spline  and  Bernstein  bases. 

We  shall  also  develop  necessary  and  sufficient  criteria  for  extending  the 
blossoming  form  of  the  dual  functionals  to  multivariate  polynomial  bases.  We 
close  with  a  few  open  questions  for  future  research. 


I  interpolation  for  solutions 
of  the  Helmholtz  equation 
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Using  Fourier  analysis  techniques  we  construct  a  Hilbert  space  whose  elements 
are  solutions  of  the  Helmholtz  equation.  A  characterization  of  the  corresponding  re¬ 
producing  kernel  is  given  explicitly.  Thus,  by  minimizing  the  corresponding  norm,  we 
construct  an  interpolant  that  passes  through  the  data  points  on  a  line.  These  values 
sample  the  solution  of  the  two-dimensional  Helmholtz  equation.  The  interpolant,  — 
that  for  any  number  of  prescribed  data  points  is  a  solution  of  the  equation —  is  given 
by  Bessel  functions  of  order  zero  and  one.  Since  these  functions  are  bell-shaped,  the 
corresponding  linear  system  of  equations  has  a  good  condition  number.  This  analysis  is 
relevant  for  the  47r-optical  sampling  problem  by  sensors  on  a  screen  and  for  holographic 
design. 

MSC:  41A05,  41A29,  35J05,  78A40 


Curvature  of  rational  quadratic  splines 
Tim  Goodman 

University  of  Dundee 
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tgoodman@dundee.ac.uk 

Abstract 

Monotonicity  of  curvature  is  investigated  for  a  rational  quadratic  segment  in  Bezier 
form.  In  particular  it  is  shown  that  for  fixed  end-points  and  tangent  directions,  provided 
that  the  tangents  do  not  turn  through  more  than  a  right-angle,  there  is  a  specific  interval 
of  weights  for  which  the  curvature  is  monotone.  An  algorithm  is  then  given  for  minimising 
the  number  of  local  extrema  of  curvature  of  a  piecewise  rational  curve  with  specified  values 
and  tangent  directions  at  knots. 
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B-Spline  Knot-Lines  Elimination  and  Continuity  Conditions 

Raul  GORMAZ 

Departamento  de  Ingenieria  Civil  Matematica 
Universidad  de  Chile 

Following  the  work  of  Dahmen-Micchelli-Seidel  [1],  where  a  new  scheme 
for  multivariate  B-spline  is  proposed,  some  insight  of  its  geometrical  nature  is 
presented.  This  scheme  for  multivariate  B-spline  provide  automatical  C'"~^ 
continuity  with  piecewise  polynomials  of  degree  k  but  the  cost  to  be  payed  is 
a  complex  structure  of  the  piecewise  polynomial  spline.  Moreover,  too  much 
flexibility  is  provided  around  de  edges  of  the  triangulation.  This  excesive 
flexibility  can  produce  irregular  surfaces. 

An  analogy  between  those  B-splines  and  triangular  Bezier  patches  is 
given.  The  affine  conditions  between  the  poles  of  two  adjacents  Bezier  patches 
for  C”'  continuity  are  replaced  by  similar  conditions  between  the  poles  of  the 
B-splines.  This  time  the  conditions  do  not  produce  additional  continuity  but 
some  knot-lines  are  eliminated.  The  “knot-lines  elimination  conditions”  can 
be  use  to  reduce  the  flexibility  of  the  scheme. 

[1]  W.Dahmen,  C.A.  Micchelli  and  H.-P.  Seidel,  Blossoming  begets  B-splines 
built  better  by  B-patches.  Math.  Comp.,  59,  (1992) 
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Applications  of  Constrained  Polynomials  to  Curve  and  Surface  Approximation - 
Author:  R.J.  Goult,  LMR  Systems,  33  Filgrave,  Newport  Pagnell,  MK16  9ET.UK 

Abstract : 

This  talk  will  describe  two  families  of  constrained  polynomials  and  their 
applications  to  problems  in  the  approximation  of  parametric  curves  and 
surfaces  arising  in  the  exchange  of  geometric  product  data. 

A  generalisation  of  the  classical  Chebyshev  economisation  method  provides  a 
method  of  degree  reduction  for  polynomial  functions  with  minimum  error  in  the 
sense  of  the  uniform  norm.  The  new  polynomial  basis  provides  this  degree 
reduction  facility  whilst  at  the  same  time  maintaining  the  value  of  the 
function  and  selected  derivatives  at  the  end  points.  The  application  of  this 
method  has  been  extended  from  simple  polynomial  functions  to  parametric 
polynomial  curves  and  bi— parametric  polynomial  surfaces.  The  method  provides 
upper  bounds  for  the  approximation  errors.  With  this  approach  the  degree 
reduction  is  achieved  in  a  seguential  manner  reducing  the  degree  by  one  at 
■  ~h  step  of  the  computation.  The  ability  to  maintain  the  values  of  end 
points  and  derivatives  makes  it  possible  to  subdivide  the  curve  or  surface 
and  obtain  polynomial  spline  approximations  of  greater  accuracy. 

The  constrained  orthogonal  polynomials  also  provide  an  optimum  approximation 
combined  with  the  ability  to  maintain  continuity  conditions  at  the  end  points 
of  parametric  curves  and  along  the  boundaries  of  bi— parametric  surfaces .  In 
this  case  the  approximation  is  optimal  in  the  least  squares  sense  and  also 
provides  an  easy  method  of  computing  the  mean  square  error  in  the  resulting 
approximation.  At  the  core  of  the  computational  method  is  the  evaluation  of 
integrals  of  products  of  the  curve  or  surface  to  be  approximated  and  the 
basis  functions.  This  means  that  although  the  technique  can  be  applied  to 
degree  reduction  problems  it  can  also  be  applied  to  any  more  general 
approximation  problem  in  which  the  curve  or  surface  to  be  approximated  is 
sufficiently  well  defined  for  the  integrals  to  be  evaluated.  The  method  has 
been  applied  to  obtain  polynomial  approximations  to  rational  curves  and 
surfaces,  to  approximate  offset  curves  and  surfaces,  to  approximate 
procedurally  defined  surfaces  and  to  obtain  a  single  span  high  degree 
parametric  polynomial  approximation  to  a  spline  curve.  A  generalisation  of 
the  method  makes  it  possible  to  a  bi-parametric  polynomial  approximation  to  a 
trimmed  portion  of  a  surface. 


Semiregular  B-spline  surfaces 
Jens  Gravesen 
Mathematical  Institute 
The  Technical  University  of  Denmark 

A  B-spline  surface  can  be  considered  as  the  surface  swept  out  by  a  B-spline  curve  with  a  fixed  knot 
sequence  and  with  the  control  polygon  moving  along  a  B-spline  curve.  If  the  knot  sequence  is  allowed  to 
vary,  we  obtain  a  new  type  of  surface  called  a  semiregular  B-spline  surface.  This  type  of  surface  has  some 
advantages  over  the  ordinary  (regular)  B-spline  surface  with  respect  to  its  interpolating  properties  and  the 
possibility  to  describe  surfaces  with  singularities. 
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Blending  functions  and  scattered  data 
interpolation  using  surfaces 
with  minimal  curvature 

Giinther  Greiner 

Department  of  Computer  Science 
Universitat  Erlangen,  Germany 


Abstract 

Finding  a  blend  (i.e.  a  rather  smooth  transition  between  primary  surfaces)  as 
well  as  scattered  data  interpolation  finally  leads  to  the  following  problem;  One 
has  to  construct  a  surface  which  satisfies  certain  (linear)  constraints  and  is  nice 
looking  (e.g.  does  not  oscillate  too  rapidly).  In  case  of  blending,  the  constaints 
are  boundary  conditions  (determined  by  the  primary  surfaces).  In  the  other  case 
interpolation  gives  the  constraints. 

We  present  a  method  based  on  variational  principles  to  construct  these  surfaces. 
The  basic  idea  is  as  follows.  In  a  (sufficiently  large)  class  of  tensor-spline  surfaces 
one  has  to  determine  the  surface  which  satisfies  the  constraints  and  in  addition 
minimizes  a  certain  functional.  The  functional  will  be  chosen  such  that  a  mixture 
of  mean  curvature  and  surface  area  become  minimal.  We  describe  the  procedure  in 
detail  and  compare  it  with  other  approaches. 
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Design  of  Blending- Surfaces  in  an 
Object-Oriented  Framework 

Giinther  Greiner'^ 

Philipp  Slusallek^ 


Abstract: 

In  computer  graphics  and  geometric  modeling  one  generally  faces  the  problem  of  integrating 
a  variety  of  curve  and  surface  types  into  a  single  program.  Blending  surfaces  are  especially 
difficult,  because  by  their  definition  they  need  to  interface  to  all  other  supported  surface  types. 
Object-oriented  design  offers  the  opportunity  of  using  the  inherent  hierarchical  structure  of 
curves  and  surfaces  to  solve  this  problem. 

This  paper  presents  an  object-oriented  framework  together  with  its  C++  implementation  that 
starts  from  an  abstract  class  of  general  differentiable  curves  and  surfaces  and  in  turn  refines  this 
design  to  curves  and  surfaces  that  are  explicitly  given  in  either  parametric  or  implicit  form.  The 
standard  curve  and  surface  types  are  derived  from  these  abstract  classes. 

An  approach  to  design  blending  surfaces  in  this  framework  is  presented.  Also,  examples  that 
visualize  the  differential  geometry  of  curves  and  surfaces  illustrate  the  approach. 


^Universitat  Erlangen,  IMMD  9  -  GDV, 
{sIusallek,greiner}@informatik.uni-erlangen.de 


Am  Weichselgarten  9,  W-8520  Erlangen,  F.R.G.,  email: 
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Abstract 

‘‘Convex  interpolation  with  minimad  energy” 

Christoph  Henninger  -  Karl  Scherer 

The  following  problem  is  considered:  let  Pq,  . . . ,  be  data  points  in  such  that  the 
polygon  passing  through  the  Pi  is  convex.  Then  determine  a  curve  F{t)  ;  t  £  [a,  6]  — » 
such  that  the  functional 

JiF)  := 

is  minimized  under  the  constraints 


/ 


liF”(t)fdt 


=  Euclidean) 


F{ti)  =  Pi,0<i<N +  1-,  k(F)  :=  F' ^  r' <  0  . 

Here  (g)  denotes  the  cross  product.  The  nodes  satisfy  a  =  Iq  <  ti  <  ...  <  fiv+i  =  b. 

It  is  shown  that  the  above  problem  does  have  always  a  solution.  Every  solution  is  - 
continuous  and  piecewise  cubic  in  each  segment  with  at  most  two  extra  knots. 

Criteria  are  given  which  allow  to  determine  their  number  and  positions  exactly.  Based 
on  this  characterization  an  algorithm  for  determining  a  soltuion  is  constructed. 
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Connection  Matrix  Splines 


Bruce  Hickey 

Department  of  Computer  Science 
University  of  Waterloo,  Canada 
and 

Hans-Peter  Seidel 
Department  of  Computer  Science 
Universitat  Erlangen,  Germany 


Abstract 

Connection  matrix  splines  are  a  generalization  of  geometrically  continuous  splines  that  are 
based  on  arbitrary  connection  matrices.  The  CMSBasis  library  implements  cubic  connec¬ 
tion  matrix  spline  curves  with  arbitrary  knot  sequences  and  arbitrary  connection  matrices 
in  C-1— 1-  as  part  of  the  Waterloo  spline  library.  The  hierarchy  in  this  implementation  is 


Curve  — ■»  LCCurve  BCurve  — »■  CHSCurve  — »  ExplCMSCurve  . 

Part  of  the  Curve  library  Part  of  the  CMSBasis  library 
Underlying  each  ExplCMSCurve  is  a  ExplCMSBasis,  which  has  the  following  ancestry 


FuncBasis  — ^  LCBasis  — >  BBasis  — > 

V - ,, - 

Part  of  the  Curve  library 


CMSBasis  — >■  ExplCMSBasis  . 

- - - / 

Part  of  the  CMSBasis  library 


The  Expl  prefix  to  the  classes  indicates  that  all  the  knot/connection  matrix  information 
is  explicitly  stored,  and  thus  variable. 


A  Local  Convexity  Preserving  Interpolation  Method 


Bas  J.  Hogervorst 


Department  of  Applied  Mathematics 
University  of  Twente 
P.O.Box  217,  7500  AE  Enschede 
The  Netherlands 


Abstract 

We  consider  the  problem  of  interpolating  scattered  3D  data  by  a  geometrically  smooth 
surface.  Assuming  that  data  points,  corresponding  normals  and  a  triangulation  are  given, 
we  propose  a  parametric  local  method  which  constructs  a  smooth  surface,  interpolating  the 
given  data,  with  the  additional  property  that  it  is  convex  if  the  data  (and  the  triangulation 
thereof)  are.  First  we  construct  boundary  curves  and  define  normals  on  them.  The 
next  thing  to  do  is  to  fill  in  each  patch,  which  we  shall  perform  by  a  kind  of  blending. 
This  results  in  a  smooth  convex  surface,  which  is  non-polynomial.  An  analysis  of  the 
parametric  method  is  given  and  some  examples  are  presented. 

E-mail:  hogervor@math.utwente.nl 


Fax:  x-31  53-356695 
phone:  x-31  53-893430  /  3460 


Mesh  Optimization 

Hugues  Hoppe 

Department  of  Computer  Science  and  Engineering 
University  of  Washington,  FR-35 
Seattle,  WA  98195 

hugues(2)cs.washington.edu 


We  present  a  method  for  solving  the  following  problem: 

Given  a  set  of  data  points  scattered  in  three  dimensions  and  an 
initial  triangular  mesh  Mq,  produce  a  mesh  M,  of  the  same 
topological  type  as  Mq,  that  fits  the  data  well  and  has  a  small 
number  of  vertices.  Our  approach  is  to  minimize  an  energy 
function  that  explicitly  models  the  competing  desires  of 
conciseness  of  representation  and  fidelity  to  the  data. 

We  show  that  mesh  optimization  can  be  effectively  used  in  at 
least  two  applications:  surface  reconstruction  from  unorganized 
points,  and  mesh  simplification  (the  reduction  of  the  number  of 
vertices  in  an  initially  dense  mesh  of  triangles). 
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Conversion  of  (trimmed)  integral  and  rational  B-Spline  surfaces 

Josef  Hoschek,  Franz-Josef  Schneider 
TH  Darmstadt 

May  be  given  a  set  of  (trimmed)  integral  or  rational  Bezier  tensor  Product 
surfaces  of  degree  (n,m)  or  a  set  of  integral  or  rational  B-Spline  surfaces  of 
order  (k,l).  These  surfaces  shall  be  converted  to  rational  or  integral  B-Spllne 
surfaces  with  respect  to  a  required  order  (k  ,1  )  and  a  minimal  number  of 
patches  with  respect  to  a  given  (distance)  error  tolerance. 

First  generic  new  boundary  curves  (knots  in  the  knot  vector)  are  introduced 
with  respect  to  the  curvature  of  the  parametric  lines  of  the  given  set  of  surface 
patches,  then  we  choose  a  suitable  set  of  points  on  the  given  surfaces  and 
approximate  these  points  with  the  discrete  norm  and  additionally  parameter 
correction. 

The  trimming  curves  are  proposed  to  be  given  in  the  parametric  domain,  the 
error  tolerance  for  the  given  and  the  approximated  trimming  curves  is  measured 
on  the  surfaces.  This  approximation  is  also  combined  with  parameter  correction. 

If  the  given  set  of  B-Spline  patches  has  C^-  or  C^-continuity.  this  continuity 
class  can  be  hold.  The  method  can  also  be  used  for  merging  of  sets  of  Bezier 
surfaces  and  data  reduction. 
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Conference  on  CURVES  AND  SURFACES 
Chamonix-Mont-Blanc,  June  10  -  16,  1993 

Riesz  Bases  in  Scattered  Data  Approximation 
Kurt  Jetter,  Universitat  Duisburg 

Stable  ils -Approximation  with  translates  ^p{--Xk),  A:  €  /,  of  a  basis  function  :  R'^ 

R,  where  X  =  {xk;k  G  Z^}  C  R'^  is  a  set  of  knots  and  7  C  is  a  finite  index  set,  requires 

that  the  translates  —  Xk),  A:  €  Z'^,  are  a  Riesz  basis  for  =  clos£,25'^,  where 

5°  :=  -Xfc);7  C  Z‘^,card(I)  <  oo}. 

k^J 

Equivalently,  in  terms  of  the  auto- correlation  function  ^  {  with  <f*{x)  =  Tp{-x) 

the  involution  of  93  ),  one  asks  for  uniform  lower  and  upper  bounds  0  <  A  <  A  for  the 
spectrum  of  the  collocation  matrices 

A/  =  -  Xk))k,iei- 

The  talks  aims  at  reporting  on  recent  results  along  these  lines  when  ip  ox  (j)  are  given  by 
preconditioned  versions  of  radial  basis  functions  which  originate  from  univariate  functions 
in  the  way  <t>{x)  =  pdla:!!')  and  g  is  an  m-th  order  completely  monotone  function  on  the 
positive  real  numbers.  In  some  aspects,  results  from  nonharmonic  Fourier  analysis  are 
useful. 
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A  Geometrical  Approach  to 
Interpolation  on  Quadric  Surfaces 


Bert  Juttler 

Technische  Hochschule  Darmstadt 
Fachbereich  Mathematik 
SchloBgartenstraBe  7 
6100  Darmstadt,  Germany 


Quadric  surfaces  like  spheres,  hyperboloids  of  one  or  two  sheets,  elliptic  and  hyperbolic 
paraboloids  are  often  used  in  industrial  applications.  The  paper  presents  a  geometrical 
approach  to  the  interpolation  of  given  points  by  rational  Bezier-  and  B-spline-curves  and 
-surfaces  on  quadrics.  Based  on  a  result  from  number  theory,  a  generalization  of  the  stereo¬ 
graphic  projection  on  the  unit  sphere  has  been  introduced  in  [Dietz  &  Hoschek  &  Juttler’93]. 
This  paper  studies  the  properties  of  this  generalized  projection. 

i.  he  generalized  stereographic  projection  is  shown  to  be  the  composition  of  the  usual  stere¬ 
ographic  projection  with  a  second  map.  The  second  map  is  called  hyperbolic  projection. 
Its  properties  are  presented  and  some  connections  to  kinematics  and  to  advanced  geometry 
cire  outlined.  The  properties  of  the  generalized  stereographic  projection  result  immedi¬ 
ately  from  those  of  the  hyperbolic  projection.  The  extension  of  the  results  to  arbitrary 
nondegenerated  quadric  surfaces  is  given. 

As  a  first  application,  the  generalized  stereographic  projection  is  applied  to  the  construction 
of  rational  Bezier  surface  patches  on  the  unit  sphere  from  given  boundaries.  Furthermore, 
interpolation  with  rational  curves  on  quadrics  will  be  shown  to  be  a  linear  problem.  More 
detained  constructions  of  curves  and  surfaces  on  quadrics  will  be  given  in  a  paper  by 
R.  Dietz  to  be  presented  at  the  3.  SIAM-Conference  on  Geometric  Design  (Tempe,  Novem¬ 
ber  93). 


References 

-  Dietz,  R.,  Hoschek,  J.  and  Juttler,  B.  (1993):  An  Algebraic  Approach  to  Curves  and 
Surfaces  on  the  Sphere  and  on  Other  Quadrics,  accepted  for  publication  in  Computer 
Aided  Geometric  Design 
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Decomposed  Surface  Editing 


Johannes  Kaasa  and  Geir  Westgaard 
SINTEF-SI,  P.O.Box  124  Blindern,  N-0314  Oslo,  Norway. 

February  3,  1993 

Abstract 

Decomposition  is  an  important  tool  for  representation  and  manipulation  of  geometry, 
with  applications  in  a  variety  of  areas.  A  main  idea  is  to  split  a  model  into  a  sequence 
of  components  with  common  properties,  such  that  each  component  can  be  handled  sep¬ 
arately.  The  decomposition  could  for  instance  be  based  on  the  level  of  details  in  the 
representation.  This  gives  a  better  flexibility  for  representation  and  manipulation.  Our 
idea  is  to  use  the  principle  in  reverse,  by  combining  geometrical  components  into  a  com¬ 
plete  form.  An  example  would  be  to  add  local  features  on  a  global  underlying  surface. 

In  this  paper  we  look  at  a  problem  often  met  in  CAD  industry.  A  master  surface 
has  to  be  modified  for  a  specific  use.  For  this  problem  we  have  developed  a  new  surface 
editing  scheme  based  on  decomposition.  We  call  it  the  ‘face  lift’  scheme.  From  the  basic 
surface  and  scattered  points  on  the  modified  surface  this  scheme  generates  the  difference 
between  the  existing  surface  and  the  desired  result.  This  ‘face  lift’  surface  should  also 
satisfy  certain  smoothness  constraints.  The  ‘face  lift’  surface  is  constructed  by  use  of 
scattered  data  techniques  and  tensor  product  B-spline  approximation.  The  illustration 
shows  ‘face  lift’  of  a  turbine  blade. 
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GENERIC  PROPERTIES  OF  THE  SET  OF  CROSS-SECTIONS  AND  THE  SET  OF 
ORTHOGONAL  PROJECTIONS  OF  A  SMOOTH  SURFACE 


Y.L.  KERGOSIEN 
Departement  de  Mathematiques 
Universite  de  Paris-Sud 
91405  Orsay,  France. 

We  give  a  complete  classification  and  description  of  the  generic  local  properties  of 
the  family  of  cross-sections  of  a  surface  and  the  family  of  orthogonal  projections  of  a 
surface.  A  smooth  2-manifold  being  given,  we  look  for  properties  that  are  generic  in 
the  space  of  its  embeddings  in  3-space.  These  properties  will  single  out  qualitative 
types  of  smooth  mappings  mainly  depending  on  types  of  singularities. 

In  the  case  of  cross-'sections,  the  position  of  a  cross-sectional  plane  depends  on  3 
parameters.  Given  an  embedded  surface  and  a  line  in  space,  the  set  of  its  cross-sections 
by  the  planes  which  are  normal  to  the  direction  can  be  studied  through  the  orthogonal 
projection  of  the  surface  on  the  line.  The  generic  properties  of  that  object  are  a  classical 
result  of  Morse  theory.  To  account  for  the  whole  family  of  cross-sections,  one  makes 
this  object  depend  on  the  direction  parameter,  to  be  taken  as  a  point  on  a  sphere  (we 
neglect  antipodal  identification  and  consider  oriented  directions  to  avoid  introducing  2- 
dimensional  projective  space).  The  locus  of  points  in  that  control  space  for  which 
Morse  theory  does  not  apply  (unstability  points,  or  bifurcation  points)  is  generically 
made  of  smooth  curves  with  cusps  and  possibly  crossings.  These  curves  are  he  images 
through  the  Gauss  map  of  the  parabolic  points  of  the  surface  under  study.  Their  cusps 
correspond  to  the  "godrons"  of  the  surface;  observing  them  requires  complete  control 
on  the  cross-sectiomil  plane  (their  codimension  is  2).The  qualitative  patterns  of  the  1- 
family  of  parallel  cross-sections  at  each  of  such  bifurcation  control  points  and  their 
unfolding  when  the  control  direction  is  perturbed  are  described.  The  bifurcation 
diagram  can  be  computed  from  the  equation  of  the  surface.  Applications  to  cross- 
sectional  medical  imaging  and  surface  coding  are  briefly  presented. 

The  case  of  projections  is  more  complex  and  we  only  summarize  the  results.  Given 
a  single  orthogonal  projection,  its  generic  properties  are  a  consequence  of  the  classical 
work" of  Whitney  on  singularities  of  mappings.  As  in  the  case  of  sections,  to  get  the 
whole  family  one  makes  the  projection  depend  on  the  oriented  direction,  i.e.  a  point  on 
a  sphere.  The  observable  qualitative  changes  in  the  projection  as  the  direction  is  moved 
are  now  due  to  two  classes  of  local  surface  configuration  :  parabolic  points  and 
swallow-tail  points.  In  control  space,  one  describes  curves  associated  either  to  beek-to 
beek  or  lip  projection  patterns  (parabolic  class)  or  to  swallow-tail  projection  pattern; 
these  curves  coirespond  to  curves  on  the  surfaceunder  study.  Codimension  2  patterns 
include  the  "godron"  (merging  parabolic  and  swallow-tail  points),  the  "gouttiere" 
(transition  between  beek-to-beek  and  lip  points),  and  the  butterfly  (swallowtail  class). 
The  bifurcation  diagram  can  be  computed  from  the  equations  of  the  surface. 
Applications  to  projection  medical  imaging  (standard  radiology)  are  presented. 
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Mxiltiresolution  Approximation  in  Statistic 
G.  KERKYACHARIAN 

Universite  de  Picardie,  Amiens,  CNRS,  URA  1321,  Paris 

It  is  very  usual  in  Statistics  to  reduce  a  non  parametric  problem  (such  as  estimating  the  common 
density  of  n  independent  random  variables)  to  a  simpler  one  of  finite  dimension  writh  the  help  of  some 
classical  results  in  approximation  theory.  In  this  talk,  we’ll  try  to  show  that  in  some  specific  problems  ,  it  is 
absolutely  necessary  to  use  an  approximation  taking  into  account  different  levels  of  resolution  at  the  same 
time.  And  we  would  like  to  emphasize  that  in  such  cases  the  wavelets  methods  appeared  to  be  much  more 
accurate  than  the  classical  ones  (  Kernel  approximations,  Fourier  series,...  0  We’ll  particularly  illustrate  this 
phenomena  in  two  special  cases  : 

We’ll  first  consider  the  density  estimation  problem  when  there  is  no  artificial  link  between  the  a  priori 
knowledge  on  the  density  and  the  measure  of  the  loss  of  an  arbitrary  procedure.  In  that  case  we’ll  show 
that  the  wavelet  framework  allows  us  to  construct  an  estimation  procedure  which  is  nearly  optimal  almost 
without  any  a  priori  knowledge. 

The  second  example  consists  in  estimating  a  non  linear  function  of  the  density  of  integral  type.  We’ll 
particularly  focuse  on  the  critical  case  occuring  for  densities  of  low  regularity.  In  this  case,  we’ll  see  that 
the  Haar  basis,  with  its  very  specific  properties  (both  of  wavelet  and  martingale  types)  appears  to  be  very 
important. 


Finding  Shortest  Paths  on  Graph  Surfaces 
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Alfred  M.  Bruckstein 
Department  of  Computer  Science 
Technion-I.I.T.,  Haifa  32000,  Israel. 


Abstract 


This  paper  presents  a  new  algorithm  for  determining  minimal  length  paths  between  two 
points  or  regions  on  a  three  dimensional  surface. 

The  numerical  implementation  is  based  on  finding  equal  distance  contours  from  a  given 
point  or  area. 

These  contours  are  calculated  as  zero  sets  of  a  bivariate  function  designed  to  evolve  so 
as  to  track  the  equal  distance  curves  on  the  given  surface. 

The  algorithm  produces  all  minimal  length  paths  between  the  source  and  the  destination 
areas  on  the  surface  given  as  height  values  on  a  rectangular  grid. 

Complexity  and  accuracy  are  governed  by  the  grid  resolution  and  the  distance  step  size 
in  the  iterative  scheme. 

Using  the  distance  maps  from  three  areas  we  also  solve  the  Steiners’problem  (for  three 
points)  on  a  surface. 

Key  Words: 

Curve  evolution,  Equal  distance  contours,  Geodesic  path. 

Partial  differential  equations,  Digital  implementation. 

Numerical  algorithms. 
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The  Construction  of  TPBC(Triplicated 
Piecewise  Bezier  Cubic) -Curve 
as  PC-Graphics  Tool 


Ha-Jine  KIMN 

Dept,  of  Computer  Engineering 
Ajou  University 
Suwon  441-749,  Korea 


ABSTRACT 


The  degree  of  Bezier  curve  constructed  with  a  set  of  discrete  plane 
points  increases  in  proportion  to  the  number  of  data  points  if  we  define  a  whole 
curve  as  one  Bezier  curve.  Thus,  in  computer  graphics  and  computer-aided 
geometric  design,  Bezier  curves  require  higher  computational  complexity. 

To  avoid  this  problem,  it  is  necessary  to  construct  a  modified 

Bezier  curve  which  is  a  piecewise  low  degree  polynomial  cmve.  For  this 

purpose,  we  bind  a  set  of  given  data  points  with  4  data  points  overlapped 
in  turn  and  construct  a  cubic  Bezier  curve  on  each  subinterval,  and  compute 
a  linear  combination  of  these  cubic  Bezier  curves  to  generate  a  smooth*  curve. 
Our  resulting  curve  is  cubic  on  the  whole  interval  and  is  an  approximate  curve 
while  keeping  the  characteristics  of  Bezier  curve.  In  fact,  we  construct  a  new 

cubic  curve  which  fits  very  closely  with  the  Bezier  curve,  which  is  called  as 

a  Triplicated  Piecewise  Bezier  Cubic(TPBC)-Curve. 


Polygonalization  Of  Real  Algebraic  Surfaces 


Reinhard  Klein 
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I  would  like  to  submit  the  following  abstract  to  the  conference 

Curves  and  Surfaces 

Chamonix  Mont-Blanc,  France,  June  10-16,  1993 


An  algorithm  is  presented  for  polygonizing  real  algebraic  surfaces.  The  algorithm  is  based 
on  an  adaptive  marching  cube  technique  combined  with  a  Bernstein-Bezier-representation 
of  the  polynomials.  Using  the  Bernstein-Bezier-representation,  an  estimation  of  the  ma¬ 
ximum  curvature  of  the  surface  in  a  tetrahedral  region  can  be  achieved.  This  estimation 
allows  to  control  the  adaptive  subdivision  in  an  easy  and  efficient  way. 

The  Bernstein-Bezier-representation  is  also  useful  to  detect  those  regions  of  the  space, 
which  do  not  contain  any  parts  of  the  surfaces.  This  property  provides  a  possibility  to 
overcome  one  problem  of  the  marching  cube  technique:  it  fails  in  cases,  where  the  surfaces 
consist  of  several  not  connected  components,  or  where  components  of  the  surfaces  are  only 
connected  at  singular  points. 

We  discuss  this  algorithm  in  detail  and  compare  it  to  existing  methods.  Several  ex¬ 
amples  will  illustrate  the  approach. 


Curves  and  surfaces  -  Chamonix  ,  10-16  Juin  1993 

Enlevement  de  noeuds  et  splines. 

Courbes  parametriques. 

Jean  Claude  KOUA  Brou 
ENST  Bretagne,  dept  MSC 
BP  832,  29285  Brest  Cedex,  France 

On  suppose  donnes  deux  entiers  naturels  n  et  k  avec  k  >  2.  Soit  t  =  une 

suite  croissante  (strictexnent)  de  reels.  On  suppose  que  I’on  dispose  de  n  +  1  vecteurs 
de  IR^*,  ou  d  est  egal  a  2  ou  3  dans  la  pratique — le  cas  des  fonctions  etant  traite 
de  la  meme  maniere  .  On  considere  la  courbe  spline  definie  par  s(x)  = 
ou  les  pour  1  <i  <n,  sont  les  B-splines  d’ordre  k  associees  a  t. 

Si  on  elimine  un  noeud  quelconque  tj  avec  k  +  1  <  j  <n  +  l  det,  un  noeud  interieur 
done,  pour  fomaer  r,  on  cherche  alors  a  determiner  une  courbe  g  approchant  s  avec  comme 
contrainte  ||y  —  s]|  <  £,  ou  e  >  0  est  une  tolerance  donnee  et  ||  .  ||  designe  une  norme 
discrete  a  preciser.  On  cherche  une  courbe  h(x)  =  AiBi^i.,T{^)  definie  a  I’aide  de  r 

et  de  s.  Dans  la  base  formee  par  les  Bi^k,t{^)  elle  se  met  sous  la  forme 


n+l 

t=l 

OU  les  Ai  sont  obtenus  a  I’aide  de  I’algorithme  d’insertion  de  noeuds  de  Boehm  en  fonction 
des  (Ai)iLi  par 

{Ai  i  <j  —  k 

aiAi  +  (1  -  Q;i)Ai_i  j  -  k  +  1  <i  <j  -1 
Aj-i  j  <  i  <  n  +  1 


ou  si  z  =  tj,  on  a  evi  =  (z  —  Ti)/(Ti+k-i  —  n). 

On  determine  les  Ai  en  identifiant  les  Ai  pour  i  <  j  —  k  ou  j  <  i  <  n  +  1  aux  Ci 
correspondant.  Les  autres  Ai  restant  sont  obtenus  pax  substitution  dans  (+). 

On  pent  ainsi  trouver  deux  courbes  hi  et  Aq,  puis  I’on  retiendra  celle  qui  est  telle  que 
I’erreur  \\hi  —  s]|  est  la  plus  petite,  et  que  Ton  note  h.  Alors  tj  sera  un  noeud  qui  peut  etre 
elimine  si  ||A  —  s]|  <  e. 

En  iterant  le  precede  sm  chacun  des  noeuds  interieurs  tj,  on  determine  un  ensemble 
^  =  (^1, .  ■  ■  ,^q)  de  noeuds  susceptibles  d’etre  elimines  et  on  forme  r  =  t  \ ^,  puis  on  trouve 
une  courbe  f{x)  =  E>iBi^k,T{^)  telle  que  ||/  — J]]  <  £■  On  peut  appliquer  le  precede 

sur  la  courbe  /,  et  ainsi  de  suite  ,  en  comparant  a  chaque  fois  avec  s,  tant  que  I’operation 
n’altere  pas  s  de  plus  de  la  tolerance  e. 


MULTIVARIATE  SPLINES  CONSTRUCTED  THROUGH 
PIECES  OF  POLYHARMONIC  FUNCTIONS 
Ognyan  Iv.  Kounchev 

In  the  simplest  case  (the  biharmonic),  the  polysplines  are 

the  solution  to  the  following  extremal  problem 

J  (2lf(x))^  dx  - >  inf, 

D 

where  the  minimization  is  over  a  proper  class  of  functions  f 

N 

such  that  f(x)  =  g(x)  for  x  £  S  =  U  T-  ,  where  g(x)  is 

j=l  J 

a  function  defined  on  N  closed  smooth  surfaces  Tj  in  D,  such 

that  T-  surrounds  T-  j,  and  =  9D.  Also  (9/9nj^)f(x)  =  0, 
J  J 

X  £  9D,  where  n.  denotes  the  inner  unit  normal  vector  to  T-. 

The  solution  s(x)  to  the  above  problem  satisfies: 

(i)  zl^s(x)  =  0  for  X  £  D\S; 

(ii)  If  we  denote  by  Sj(x)  the  restriction  of  s(x)  to 

the  layer  lying  between  Tj  and  ~  Sj(x), 

(a/anj.)Sj.^j(x)  =  (a/anj)Sj(x),  Zls^^^Cx)  =  ^s^Cx),  all  for  x  £ 
Tj,  j  =  1,...,N-1. 

1.  Existence  and  uniqueness  of  polysplines  is  proved  in  proper 
Sobolev-Slobodeckii  spaces. 

2.  Convergence  is  proved  when  N - ^  oo  and  the  data  function 

g  is  defined  in  the  whole  of  the  domain  D. 

3.  Due  to  the  Green’s  formulas  the  radial  basis  generated 
through  the  fundamental  solution  of  the  polyharmonic  equation  is 
dense  in  the  polysplines  defined  above.  Correspondingly,  density 
results  for  such  radial  basis  are  proved  through  point  2. 

Research  partially  sponsored  by  the  Alex,  von  Humboldt  Foundation 
Address:  Dept.  Math.,  University  of  Duisburg,  Lotharstr.65,  4100 
DUISBURG,  GERMANY;  E-mail:  hn314ko@unidui.uni-duisburg.de 
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On  approximation  by  bivariate  incomplete  polynomials 


An  dr  as  Kroo 


Let  AT  C  be  a  star  -  like  domain  and  consider  polynomials  of  the  form 

V{^^y)=  au  , 

dn<k-^lKn 


where  0  <  ^  <  1  is  a  fixed  number.  These  polynomials  are  bivariate  analogues  of 
i^orentz’  ^-inc.  omplete  polynomials.  The  univariate  theory  implies  that  in  order 
that  /  €  C{K)  can  be  approximated  by  these  polynomials  it  is  necessary  that 
/  =  0  on  6~K.  In  this  paper  we  study  the  following  question:  can  any  /  €  C{K) 
such  that  /  =  0  on  6^K  be  approximated  by  polynomials  of  the  form  (=♦=)?  It  turns 
out  that  if  K  is  a  countable  intersection  of  triangles  with  a  vertex  in  the  origin 
then  the  answer  to  the  above  question  is  affirmative.  Thus  the  admissible  domains 
include  triangles  and  parallelograms,  certain  polygons  and  sections  of  discs.  The 
proof  is  based  on  combining  Miintz-type  approximation  with  approximation  by 
bivariate  Bernstein  polynomials  on  triangles. 

We  shall  also  discuss  extension  of  these  results  to  domains  of  higher  dimension 
(countable  intersections  of  r-dimensional  simplices  with  a  vertex  at  the  origin). 
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On  the  Fast  Evaluation  of  Integrals  of  Wavelets  and  Box  Splines 


Angela  Kunoth 
Institut  fur  Mathematik  1 
Freie  Universitat  Berlin 
Arnlmallee  2-6 
D-W-1000  Berlin  33 
Germany 

Abstract 

For  solving  partial  differential  equations  by  a  Galerkin  method,  one  has 
to  compute  the  entries  of  mass  and  stiffness  matrices  which  are  integrals  of 
(derivatives  of)  the  trial  functions.  In  the  concept  of  refinable  shift-invariant 
spaces,  these  expressions  can  be  determined  quickly  and  exactly  (up  to  round¬ 
off)  by  solving  a  simple  eigenvector  problem.  Here,  we  will  focus  on  the 
algorithm  and  the  application  to  biorthogonal  wavelets  and  different  two- 
dimensional  box  splines. 

In  addition,  since  integrals  of  box  splines  are  again  box  splines,  we  will  use 
the  procedure  for  the  fast  evaluation  of  box  splines,  in  particular,  with  many 
direction  vectors. 


MONOTONIZING  PARAMETRIZATION  FOR  INTERPOLATIONAL  SPLINES 
B. I. Kvasov  and  S.A.Yatsenko 

Russian  Academy  ol  Sciences,  Institute  ol  Computational 
Technologies,  Novosibirsk,  630090  , Russia 
Tel.  383-235-3745,  Pax:  383-235-1242,  E-mail;  kvasov<8ict. nsk.su 

Let  P^=(x^,y^),  i=0,1,...,N  be  a  sequence  of  pairwise  different 
points  in  the  plane  XY.  In  order  to  draw  a  curve  passing  through 
these  points  in  general  case  it  is  necessary  to  construct  a  mesh 
A:  a=tQ<t<j<. .  .<tjj=b  and  to  define  a  continuous  vector-function 
G{t)=(C^(t),Cy(t)),  te[a,b3  such  that  Cx^'^i^=^i»  ^y^‘‘^i^=yi* 
i=0,1,...,N,  i.e.  C(tj|^)=P^,  i=0,1,...,N.  The  form  of  the  curve 
depends  both  on  the  knot  selection  and  on  the  interpolation  method. 

In  this  paper,  the  notion  of  a  pararoetrization  correspondence 
with  the  given  data  has  been  introduced.  We  propose  an  algorithm  of 
data  dependent  parametrization  guaranteeing  a  preservation  of  the 
monotonicity  properties  for  sufficiently  arbitrary  initial  data. 
It  is  based  on  the  proved  invariance  properties  in  a  parametric 
space  of  n  degree  interpolating  polynomials  and  splines.  We 
establish  a  close  connection  of  the  suggested  parametrization 
knots  choice  with  sufficient  conditions  of  monotonicity  for  cubic 
splines.  The  algorithm  is  generalized  for  the  case  of  an  multi-form 
surface  approximation  using  NURBS.  The  given  results  of  the 
numerical  calculations  exhibit  our  method  as  more  preferable  than 
generally  accepted  uniform  and  accumulated  chord  length  parametri- 
zations. 

AOt 
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Computations  of  Curvatures  Related  to  Surface-Surface  Filleting/Chamfering. 

Abstract 

Fillet  (and  chamfer)  constnictions  between  two  parametric  surfaces  usually  employ  tracing 
algoritlims.  Accuracy  conditions  require  careful  choice  of  the  step  size  at  every  iteration  of  the 
algorithm.  To  determine  the  step  size,  one  has  to  know  exact  curvature  values  for  contact  and 
.spine  curves  in  tire  object  space  as  well  as  in  parameter  spaces  of  the  surfaces. 

The  paper  shows  how  to  compute  these  values  in  cases  of  constant  radius  fillets  and  general  fillets. 

Genei  al  fillets  can  be  defined  either  locally,  when  every  fiUet  arc  radius  depends  only  on  local 
pro|reities  of  the  siufaces.  or  globally,  when  the  radius  depends  on  the  position  of  the  arc  with 
rosiiect  to  the  entire  fillet  surface.  Constmction  and  multiple  parametrizations  of  the 
corresponding  Voronoi  surface  are  instrumental  in  the  derivations. 

Similar  derivations  cmi  be  applied  to  construction  of  non-isosceles  (elliptical,  etc.)  fillets  and 
chamfers. 
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LEAST-SQUARES  OPTIMIZATION  OF  THREAD  SURFACES 

Michel  LEGER,  Inslitut  Frangais  du  Petrole,  BP  311,  92506  Rueil-Malmaison  Cedex,  FRANCE 

Jesn-Marie  MORVAN,  Universite  Claude  Bernard  LYON  /,  69622  Villeurbanne  Cedex,  FRANCE 
Muriel  THIBAUT,  Insdtut  Frangais  du  Petrole,  BP  311,  92506  Rueil-Malmaison  Cedex,  FRANCE 

In  geology,  faults  may  sometimes  be  considered  as  slipping  surfaces  between  two  (approximately) 
rigid  blocks.  Geometrically,  we  define  these  surfaces  as  follows:  they  are  preserved  by  a  (at  least) 
single-parameter  family  of  moves,  and  we  call  them  threads  by  analogy  with  a  screw  and  a  nut. 
Besides,  we  need  to  optimize  our  subsurface  models  for  oil  exploration  purposes.  Therefore,  we 
propose  a  method  that  optimizes  a  surface  with  respect  to  several  least-squares  criteria  related  to 
curvature,  proximity  to  known  points,  and  the  “thread  property”. 

This  definition  of  threads  is  equivalent  to  following  a  characteristic  property  [suggested  by 
C.  Marie]:  there  exists  at  least  one  nonzero  twistor  IT  such  that  Tl'(P)  £  TpS  at  any  point  P 
on  the  surface  5,  with  TpS  being  the  plane  tangent  to  S  at  P. 

We  represent  the  unknown  surface  S  parametrically  by  a  function  [x.  y.  z)  =  u.  v),  which  we 
discretize  using  B-spline  tensor-products.  The  optimization  problem  consists  in  minimizing  the  overall 
objective  function  -fi  where  twistor  IF  is  an  auxiUiary  unknown. 

is  the  weighted  geometrical  objective  function  =  ITp^-Qp^-  +  WcQc  +  1'1'rQr- 
Qp,r  =  E,  IKP-^Oill'  measures  the  length  of  the  normal  part  of  the  vectors  that  link  known 
points  P^  and  arbitrarily  chosen  points  iv)  of  a  parameterization  $  of  the  surface. 

Q(^-  —  (A'f  -i-  Xr,)dS,  where  Ai  and  Ao  are  the  principal  curvatures  at  each  point  on  the  surface. 
Qj  =  <  N{P).  TF(P)  > -dS,  where  <  -V{P).  IF(P)  >  is  the  scalar  product  at  P  of  unit  normal 
N  and  twistor  I'F. 

=  Wp.:,Qp,,  -f  WdQd  -h  IF.vQ.v  is  an  additional  objective  function.  Objective  functions 
Qp^  =  ||(P,  .  Qd  ^re  necessary  to  constrain  those  of  the  degrees  of  freedom  in  the 

space  of  the  parametric  representations  which  change  function  $  if  they  are  perturbed,  but  preserve 
the  surface  itself,  is  a  smoothing  term  based  on  $’s  second  derivatives.  Twistor  IF  needs  to  be 
nonzero,  because  any  surface  would  be  a  thread  without  this  condition;  this  leads  to  a  normalization 
objective  function  O.v  =  ||I  F]|-(/6’ )/(  /_<;  (/P  )- i]‘.  A  proper  choice  of  weights  IFp,.  .  .  .  enables 

these  criteria  to  be  adequately  balanced. 

The  implementation  of  this  optimization  problem  using  the  Gauss-Newton  procedure  leads  to 


Ao 


satisfactory  numerical  results. 


Approximating  Reachable  Sets  by  Extrapolation  Methods 


Curves,  surfaces  and  higher  dimensional  manifolds,  which  are  implicitly  defined  as 
submanifolds  of  reachable  sets  of  controlled  dynamical  systems,  comprise  a  challenging 
object  of  approximation  methods.  In  this  talk,  our  main  interest  lies  m  difference 
methods,  especially  in  order  of  convergence  results,  for  the  discrete  approximation  of 
reachable  sets  with  respect  to  Hausdorff  distance. 

We  concentrate  on  a  special  approach  for  the  numerical  approximation  of  reachable 
sets  of  linear  differential  inclusions.  This  approach  is  based  on  the  computation 
of  Aumann’s  integral  for  set-valued  mappings.  It  consists  in  exploiting  ordinary 
quadrature  formulae  with  nonnegative  weights  for  the  numerical  approximation  o 
the  dual  representation  of  Aumann’s  integral  via  its  support  functional.  For  set¬ 
valued  integrands  which  are  smooth  in  an  appropriate  sense,  this  approach  yields 
integration  methods  of  arbitrarily  high  order,  e.g.  extrapolation  methods  based  on 
Romberg  integration.  For  problems  which  are  not  smooth  enough  in  that  sense,  this 
approach  yields  further  insight  into  well-known  order  reduction  phenomena. 

Most  important  are  adaptions  of  these  extrapolation  methods  to  linear  differential 
inclusions.  As  a  result,  we  get  higher  order  methods  for  the  discrete  approximation 
of  reachable  sets  of  special  smooth  classes  of  linear  control  problems. 

The  talk  finishes  with  visualizations  of  the  numerical  results  for  several  model 
problems  and  with  an  outline  of  open  questions  and  possible  directions  of  future 

research. 
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Convolution  Kernels  Based  on  Thin-Plate  Splines 


Jeremy  Levesley 

Department  of  Mathematics  and  Computer  Science 
Leicester  University 
Leicester  LEI  7RH,  UK 
jll@uk.ac.le.mcs 

A  radial  basis  function  is  one  of  the  form  ^ojUl,  where  <j)  :  ^  R  is  continuous 

and  nil  is  the  Euclidean  norm.  One  method  of  constructing  approximations  to 
/  €  C(R")  is  quasi-interpolaiion.  Here,  we  approximate 

fi^)  ~  {l~^)  ’ 

z  €  2^ 

for  some  dilation  parameter  h,  where  ^  is  a  radial  basis  function.  In  order  that 
this  approximation  converges  we  require  that  V’  is  integrable  with  integral  equal 
to  1.  In  this  paper  we  seek  functions  rp  of  the  form 

V.(x)  =  f  J—  (  <p{\\x-tu\\)dS^-Hu)dt^{t), 

Jr.  ^n-i 

where  /j.  is  some  appropriate  measure,  5"“^  is  the  unit  n-sphere,  w„_i  its  surface 
area,  and  ^  is  a  function  of  the  form 

CO 

<f>{r)  =  log  r^afcr““*,  r  >  A  >  0,  a  > —n. 

k=0 

We  give  conditions  on  the  measure  n  which  ensure  that  tp  is  integrable  and  show 
that,  under  such  conditions,  ^  has  a  non-zero  integral  only  in  a  very  specific  set 


^  OG 


of  circumstances. 


Vertex  splines  for  smooth  curve/surface  interpolation  with  local  tension 


David  Levin 

School  of  Mathematical  Sciences 
Tel  Aviv  University 
Ramat-Aviv,  Tel  Aviv  69978 
ISRAEL 

Abstract. 

A  new  simple  method  for  curve  interpolation  is  presented.  The  method  is  based  on 
vertex  splines,  it  is  local,  and  it  also  has  local  tension  parameters.  The  method  has  a 
natural  extension  to  surface  interpolation  of  general  (non-gridded)  data.  It  is  a  simple 
k-smooth  local  surface  interpolation  scheme  with  a  second  order  approximation  order,  and 
local  tension  parameters.  A  unique  feature  is  the  possibility  of  forcing  discontinuities  along 
edges,  or  at  vertices. 


Approximation  with  thin  plate  splines  and  exponential  radial  basis 

functions 

Will  Light,  University  of  Leicester 

Over  the  last  ten  years  there  has  been  tremendous  progress  in  the  theoretical  understanding 
of  approximation  by  radial  basis  functions.  These  functions  give  methods  of  approximating 
either  functions  or  data  in  R"  by  linear  combinations  of  functions  with  very  simple  structure. 
Let  xi, . .  be  points  in  IR”.  Then  the  radial  basis  approximating  subspace  associated 
with  this  collection  of  points  is  the  linear  span  of  the  set  of  functions 

WII  • 

Here  ^  is  a  function  from  IR'*'  to  IR.  Suitable  choices  for  4>  usually  depend  on  the  underlying 
dimension  of  the  domain.  For  example,  if  n  is  odd,  then  the  multiquadrics  ^(r)  =  y/r^~+~c, 
c  >  0  are  a  good  choice.  If  n  is  even,  then  the  thin  plate  splines  4>(r)  =  r^ln(r)  perform 
weU.  (It  should  be  observed  that  we  are  in  the  land  where  the  thin  plate  spline  was  born, 
and  that,  as  the  first  of  the  radial  basis  functions  to  be  looked  at  in  detail,  it  received 
curiously  little  attention  at  birth.  However,  it  is  a  happy  reality  that  the  thin  plate  splines 

. ived  some  apalling  post-natal  neglect,  and  have  gone  on  to  become  a  beautiful  part  of 

approximation  theory!) 

Practical  experience  has  shown  that  radial  basis  functions  are  good  for  handling  approxi¬ 
mation  problems  where  the  data  sampling  is  carried  out  at  a  ‘scattered’  or  badly  distributed 
set  of  points.  They  also  perform  well  for  high  dimensional  approximation  problems.  Once 
the  dimension  n  gets  much  above  3,  it  becomes  very  difficult  to  use  many  of  the  more 
sophisticated  triangulation  methods,  but  radial  basis  functions  seem  in  some  ways  to  be 
impervious  to  high  dimensions. 

The  theoretical  understanding,  and  some  of  the  practical  aspects  of  developing  radial 
basis  approximations  are  stiU  lacking  in  several  important  features.  In  this  talk  we  address 
two  of  these  areas.  The  first  concerns  the  approximation  properties  of  radial  basis  functions 
when  the  associated  points  are  evenly  distributed  in  a  compact  domain.  The  main  satis¬ 
factory  treatment  until  now  has  been  for  evenly  distributed  points  throughout  the  whole  of 
IR" ,  where  the  techniques  of  Fourier  analysis  can  be  brought  to  bear.  We  illustrate  a  line 
of  thinking  which  allows  quasi-interpolants  to  be  constructed  for  compact  domains,  once 
the  associated  operator  is  known  for  the  whole  of  ]R".  The  second  concerns  approximation 
using  functions  which  decay  at  infinity.  A  peculiar  dichotomy  seems  to  develop,  brought 
on  by  the  use  of  Fourier  methods.  If  the  radial  basis  function  has  growth  at  infinity,  there 
is  some  chance  of  taking  linear  combinations  which  decay  at  infinity,  and  which  satisfy  the 
Strang-Fix  conditions.  These  conditions  guarantee  that  the  associated  quasi-interpolant 
reproduces  polynomials,  from  whence  come  the  desired  error  estimates  by  the  familiar  Tay¬ 
lor  series  arguments.  If  the  radial  basis  function  has  decay  at  infinity,  then  this  argument 
seems  to  fail.  If  time  permits,  we  will  illustrate  two  ways  in  which  this  apparent  difficulty 
is  overcome. 

I  appreciate  that  the  audience  will  have  a  wide  range  of  interests  and  backgrounds,  and 
will  try  to  make  this  interesting  area  of  research  accessible  to  all! 
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Change  of  Basis  Algorithms  for  Surfaces  in  CAGD 

Suresh  Lodha 

Department  of  Computer  and  Information  Sciences 
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Ron  Goldman 

Department  of  Computer  Science 
Rice  University 


Abstract 

B-patches  were  introduced  by  Seidel  and  shown  to  agree  with  the  multivariate 
B-splines  on  a  certain  region.  Lineal  polynomials,  that  is,  polynomials  formed  by 
products  of  linear  polynomials,  were  introduced  by  Cavaretta  and  Micchelli  and  shown 
to  be  dual  to  B-patches  using  a  multivariate  polynomial  identity.  This  work  establishes 
the  duality  between  B-patches  and  lineal  polynomials  by  generalizing  the  deBoor-Fix 
formula  for  curves  to  surfaces.  This  formula  is  then  used  to  establish  duality  between 
several  important  multivariate  bases  for  representing  surfaces  including  the  duality 
between  Bezier  and  multinomial  bases,  Lagrange  and  power  bases,  and  Newton  and 
Newton  dual  bases. 

Our  purpose  here  is  to  establish  the  change  of  basis  algorithms  between  any  two 
of  these  bases.  This  is  achieved  by  employing  three  basic  techniques:  homogenization, 
blossoming  and  the  duality  principle.  The  underlying  theme  is  unification:  unification 
of  algorithms,  unification  of  derivations,  and  unification  of  representations.  This  ap¬ 
proach  has  the  advantage  of  providing  computational  simplicity  while  at  the  same  time 
revealing  the  beautiful  geometric  and  algebraic  structure  of  these  various  representa¬ 
tions  for  surfaces. 

Our  approach  provides  a  single  consistent  framework  for  many  recursive  algorithms, 
which  are  central  to  surface  representations.  We  shall  demonstrate  that  these  algo¬ 
rithms  including  evaluation,  subdivision,  differentiation,  degree  raising  and  knot  inser¬ 
tion  algorithms  are  different  manifestations  of  the  change  of  basis  algorithms.  Although 
we  present  the  change  of  basis  algorithms  for  surfaces  for  the  sake  of  simplicity,  this 
theory  extends  to  higher  dimensional  hypersurfaces  in  a  straightforward  way. 

Keywords:  algorithms,  surfaces,  CAGD,  deBoor-Fix,  basis  change,  evaluation,  differ¬ 
entiation,  knot  insertion,  degree  elevation,  subdivision,  duality. 


Approximation  of  Circle  Segments 
by  Polynomial  Curves 

Tom  Lyche 
University  of  Oslo 

Abstract.  In  a  joint  work  with  K.  M0rken  we  consider  approximating  a  seg¬ 
ment  of  a  circle  by  parametric  curves  (ar(t),  y(t)),  where  x{t),y{t)  are  polyno¬ 
mials  in  t.  For  a  segment  of  length  h  it  is  known  that  for  degree  n  <  3  there 
exist  approximations  with  error  We  extend  this  result  to  arbitrary  n. 


MATHEMATICAL  MODELLING  OF  FREE-FORM  SURFACES 
FROM  DIGITIZED  DATA  OF  COORDINATE  MEASURING  MACHINES 
USING  NON-UNIFORM  RATIONAL  B-SPLINE  SURFACES 

Weiyin  Ma  and  J.P.  Kruth 

Katholieke  Universiteit  Leuven 

Division  of  Production  Engineering,  Machine  Design  and  Automation 
Celestijnenlaan  300B,  B-3001  Heverlee,  Belgium 

ABSTRACT 

This  paper  introduces  theoretical  and  practical  approaches  for  mathematical  modeling  of 
free-form  surfaces  from  digitized  data  obtained  with  three  dimensional  coordinate  measuring 
machines  (CMMs).  Non-Uniform  Rational  B-SpHnes  (NURBS)  surfaces  are  used  to  represent 
the  free-form  shape.  Least  squares  fitting  techniques  are  applied  to  identify  the  parameters  of 
a  NURBS  surface  from  discrete  points. 

The  procedure  begins  with  the  parameterization  of  digitized  points  (or  allocation  of  u  and 
V  parameters  to  each  of  the  digitized  points)  using  a  first  approximation  of  the  free-form  shape. 
The  first  approximation  is  called  a  base  surface  for  mathematical  modelling  and  is  created  from 
approximate  boundary  information  of  the  underlying  free-form  surface.  The  parameterization  is 
realized  by  projecting  the  digitized  points  to  the  base  surface.  The  parameters  of  the  projected 
points  will  be  used  as  the  parameters  of  the  digitized  points.  After  the  parameterization  of 
digitized  data,  the  degree,  the  number  of  control  points  and  the  knot  parameters  are  selected 
according  to  the  complexity  of  the  underlying  shape.  The  selection  of  the  above-mentioned 
parameters  is  flexible  and  can  be  made  automatic  via  an  iteration  process.  With  the  above 
information  a  set  of  simultaneous  observation  equations  is  constructed.  The  unknown  variables 
in  these  observation  equations  are  the  three  dimensional  control  points  and  the  weights  of 
these  control  points.  Householder  transformations  is  then  applied  to  achieve  the  least  squares 
solutions  of  the  observation  system. 

Algorithms  presented  in  this  paper  have  been  implemented  on  a  VAX/VMS  platform.  A 
software  package  SHAPID  (SHAPe  IDentification)  is  developed  for  this  purpose.  At  the  present 
phase  SHAPID  is  mainly  composed  of  three  functional  modules,  i.e.,  user  and  graphical  inter¬ 
face,  database  management  module  and  an  interactive  least  square  processing  module.  There 
are  also  some  other  supporting  modules  for  input  and  output,  entity  manipulation,  simple 
entity  editing,  etc.  SHAPID  has  been  integrated  inside  a  commercially  available  CAD/CAM 
system  UNIGRAPHICS  (a  product  of  EDS,  formerly  of  McDonnell  Douglas  Corporation)  for 
better  graphical  display. 

This  research  work  is  sponsored  by  the  Commission  of  the  European  Communities  through 
a  Brite/Euram  project  BE-4527,  PROBES:  Surface  Modelling  and  Quality  Control  System 
for  Manual  and  Automatic  Coordinate  Measuring  Machines,  and  the  Katholieke  Universiteit 
Leuven  through  a  doctoral  scholarship. 


A  Method  for  Summability  of  Lagrange  Interpolation 

On  P.L.  Butzer’s  Problem 


Detlef  H.  Mache 


By  using  a  technique  from  numerical  integration  and  the  new  0  -  transformation  (see  [2]) 
we  get  a  discretely  defined  method,  which  is  a  modification  of  the  Lagrange  projection.  We 
will  consider  the  modified  summability  methods  of  Ln  given  by 

n 

j=0 


with  the  generating  polynomial  (see  [2]) 

(l  -  a:  cos  (1  -Tn+2{x))  _  " 


5n(a:)  =  K„- 


(1  -x)[x-  cos  i=o 


where  sin^  and 


^  n-i  +  2  2 

^j,n  -  - — cos  - + 


n  +  2 


- cos  sin 


J-K 


n 


+  2  (n  +  2)  sin  n  +  2  n  +  2 


For  this  discrete  linear  positive  approximation  method  we  have  an  instance  of  alge¬ 
braic  polynomials,  for  which  we  can  prove  estimates  of  Timan  and  Telyakovskii  -  Gopengauz 

type-  .  .  , 

It  is  the  aim  of  this  present  note  to  show  that  the  local  order  of  approximation  by  means 
of  (Z/*)  is  comparable  with  those  furnished  by  the  best  approximation  polynomials. 


At  the  end  of  my  explanations  we  get  a  conclusion  to  prove  an  problem  of  P.L.Butzer 
presented  on  a  conference  in  Budapest  [1]. 
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e-G^  Polynomial  Approximations  of  Non-polynomial 

Surfaces 

Stephen  Mann 

In  CAGD,  there  are  variety  of  surface  representations  that  are  mathematically  complex 
or  computationally  expensive  to  evaluate.  For  example,  the  offset  to  a  polynomial  surface 
is  a  non-polynomial  surface;  S-patches  are  an  n-sided  surface  patch  of  high  rational  degree. 
Additionally,  it  is  sometimes  desirable  to  have  a  parametric  polynomial  representation  of 
an  algebraic  surface.  Unfortunately,  while  parametric  polynomial  surfaces  always  have  an 
algebraic  representation,  the  converse  is  not  always  true.  In  general,  there  is  no  parametric 

polynomial  representation  for  an  algebraic  surface. 

In  this  talk,  I  will  show  how  to  use  the  cubic  interpolant,  a  triangular,  high-order-of- 
approximation  cubic  polynomial  patch,  to  approximate  these  types  of  surfaces.  In  order 
to  use  this  low  degree  surface  patch,  I  relax  the  continuity  conditions  from  to  e-G^  (a 
surface  is  said  to  be  e-G^  if  the  discontinuity  in  surface  normal  is  bounded  by  e).  Despite  the 
discontinuities  in  surface  normals,  the  approximations  are  indistinguishable  from  the  original 


surfaces. 


comonotone  Hermite  interpolation  via  parametric  cubics 


Carla  Manni  -  JJniversita’  di  Pisa 


A  visually  pleasing  interpolating  function  should  generally  retain  the  intrinsic  form 
that  is  implied  by  the  given  data  points.  That  is  if  the  data  are  locally  monotone,  then 
the  interpolating  function  should  be  also.  A  great  deal  of  research  has  been  done  in  recent 
years  in  constructing  comonotone  interpolants. 

In  this  paper  we  assume  that  a  set  of  points  xo  <  Xi  <  . . .  <  Xn  and  the  values 
of  a  function  f{xi)  are  given,  in  addition  we  assume  that  the  derivatives  /'{xi)  at  the 
data  points  are  known  (given  or  computed):  we  seek  to  find  a  function,  s(a:)  €  C'^[a;o,  a:^] 
interpolating  /  and  its  first  derivative  at  the  data  points,  which  is  increasing  (decreasing) 
iir  [xi,Xj-+i]  if  the  data  are. 

It  is  well  known  that  this  problem  can  not  be  solved  by  piecewise  polynomials  with 
preassigned  degree  since  the  derivatives  of  such  a  function  (s'(xi),  s'(xi4.i))  must  belong 
to  a  bounded  domain  of  the  plane  in  order  to  obtain  monotonocity. 

So  the  problem  can  be  faced  via  piecewise  polynomials  with  fixed  degree  adding  extra 
knots,  via  piecewise  polynomials  having  degree  dependig  on  the  data  or  via  piecewise 
rational  functions. 

In  this  paper  we  situate  the  problem  in  the  general  setting  of  parametric  curves:  we 
propose  a  solution  constructed  considering  the  graph  of  the  function  x  s-s  the 

support  of  a  bidimensional  parametric  curve  (x(t),y(t))  where  x(t)  is  strictly  increasing. 

This  point  of  view  is  rather  new;  in  the  papers  following  this  approach  the  shape 
constraints  are  essentially  obtained  in  three  ways:  considering  the  variation  diminishing 
properties  of  the  Bernstein  polynomials,  so  that  only  sufficient  conditions  are  considered, 
minimizing  suitable  functionals,  which  implies  global  schemes,  or,  finally,  varing  heuristi- 
cally  the  parametrization  of  the  curve. 

The  method  proposed  here  constructs  the  required  interpolant  s(x(t)),  x  £  [xi,Xj-+i], 
as  a  single  valued  parametric  curve  having  cubic  components.  The  comonotonicity  of 
the  curve  is  controlled  by  the  amplitude  of  the  tangent  vector  at  the  extrema  of  each 
subinterval;  this  amplitude  is  determined  by  using  the  well  known  necessary  and  sufficient 
conditions  for  the  derivatives  of  monotone  cubics. 

The  method  yelds  a  second  order  Loo  approximation  to  a  monotone  function,  while 
yelds  a  fourth  order  approximation  to  a  function  with  no  vanishing  first  derivative. 

Finally  the  proposed  method  is  local:  a  change  in  the  data  at  the  point  Xi  reflects 
only  in  [xi_i,xi+i]. 


Evaluating  Surface  Intersections  in  Lower  Dimension 

Dinesh  Manocha  Amitahh  Varshney  Hans  Weber 
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Evaluating  the  intersection  of  parametric  and  algebraic  surfaces  is  a  recurring  operation 
in  geometric  and  solid  modeling.  This  problem  has  been  studied  extensively  in  the  literature 
and  the  main  approaches  are  based  on  subdivision,  analytic  and  tracing  methods.  For  rational 
and  algebraic  surfaces,  the  intersection  problem  can  be  formulated  using  algebraic  sets  and 
in  the  last  few  years  a  great  deal  of  attention  has  been  paid  to  evaluating  these  curves  using 
marching  methods.  Robust  algorithms  for  marching  methods  need  to  know  start  points  on 
each  component,  singular  points  on  the  intersection  curve  and  optimum  step  sizes  as  we  are 
tracing  along  the  curve.  When  it  comes  to  tracing,  we  can  choose  to  trace  the  curve  in  the 
higher  dimensional  space  (spanned  by  the  parameters)  or  its  projection  in  the  lower  dimensional 
space,  obtained  by  eliminating  variables  in  the  system  of  algebraic  equations  [1].  Earlier  the 
lower  dimensional  approach  was  considered  infeasible  due  to  efficiency  and  numerical  problems. 
In  our  earlier  work,  we  had  shown  the  lower  dimensional  formulation  of  the  intersection  can 
be  represented  as  the  determinant  of  a  matrix  [2].  A  simple  tracing  algorithm  proceeds  by 
substituting  the  numerical  values  in  the  matrix  and  computing  the  resulting  determinant  using 
Gaussian  elimination  and  using  a  modification  of  Gaussian  elimination  to  compute  the  partials 
of  the  function  (for  tangent  direction  computation)  [2]. 

In  this  talk,  we  highlight  a  new  algorithm  for  tracing  the  intersection  curve  using  the  matrix 
formulation.  The  projection  of  the  intersection  curve  is  represented  as  the  determinant  of  a 
matrix,  M{u,v).  Given  this  formulation,  we  use  the  property  that  only  those  points,  (ui,Vi) 
correspond  to  the  intersection  curve  such  that  M(ui,vi)  is  singular.  The  resulting  algorithm 
for  evaluating  the  intersection  curve  is  based  on  eigendecompositions  and  singular  value  decom¬ 
positions,  as  opposed  to  determinant  computation,  which  improves  its  accuracy.  Furthermore, 
at  each  stage  of  the  tracing  step,  we  reduce  the  problem  to  computing  the  nearest  eigenvalue 
and  eigenvector  using  inverse  power  iteration.  We  also  make  use  of  the  structure  of  the  matrices 
in  performing  these  power  iterations.  A  major  advantage  of  inverse  power  iterations  lies  in  the 
fact  that  its  convergence  is  well  understood  as  opposed  to  applying  local  methods  (like  New¬ 
ton’s  method)  to  solve  a  system  of  polynomial  equations.  As  a  result,  we  are  able  to  develop 
a  robust  tracing  method  using  matrix  computations  and  geometric  reasoning  in  a  lower  dimen¬ 
sional  space.  We  describe  the  implementation  of  this  algorithm,  its  performance  on  a  number 
of  examples  and  our  experiences  in  tracing  curves  with  singularities. 

[1]  C.  Hoffmann,  “A  dimensionality  paradigm  for  surface  interrogations” ,  Computer  Aided 
Geometric  Design,  vol.  7,  pp.  517-532,  1990. 

[2]  D.  Manocha  and  J.  Canny,  “A  new  approach  for  surface  intersection”.  International 
Journal  of  Computational  Geometry  and  Applications,  vol.  1,  number  4,  pp.  491-519,  1991 
(special  issue  of  solid  modeling) . 


CONVERGENCE  OF  NON-STATIONARY  NON-NEGATIVE  SUBDIVISION 

Avraham  A.  Melkman^ 

Abstract:  The  first  part  of  this  talk  is  based  on  joint  work  with  Alfred  Cavaretta. 

Following  Buhmann&:  Micchelli  we  view  subdivision  as  the  iteration  =  Au”,  =  v,  with 
V  £  Here  A  =  is  a  bi-infinite  matrix  which,  in  the  terminology  of  Dahmen& 

Micchelli,  is  two-slanted,  i.e.  for  some  £  >  1 

Gij  =  0  unless  0  <  i  —  2j  <  1. 

The  usual  stationary  subdivision  is  obtained  in  case  =  ai+2,j+i-  More  generally,  call  A  2k- 
stationary  if 

0>i-^2k,j+k  —  ^i,ji  all 

.  .US  is  the  stationarity  associated  with  subdivision  using  &  k  X  k  matrix  mask. 

When  the  subdivision  coefficients  are  non-negative  the  convergence  properties  of  the  scheme 
do  not  depend  on  the  actual  values  of  its  coefficients  but  rather  on  its  stencil,  cr{A),  defined  to  have 
entries  (cr(A))ij  =  1  if  aij  >  0,  and  0  otherwise. 

We  prove  that  a  non-negative  subdivision  scheme  whose  stencil  is  2A;-stationary  converges  if 
and  only  if  there  exist  an  N  and  c>  0  such  that  every  £  consecutive  rows  of  A^  contain  a  column 
all  of  whose  elements  are  bounded  from  below  by  c. 

Specializing  to  the  case  of  stationary  subdivision,  this  theorem  leads  to  a  simple  necessary  and 
sufficient  condition  for  convergence  couched  only  in  terms  of  the  support  of  the  coefficients  of  the 
mask.  To  be  precise,  let  Cij  =  ai-2j  with  aoae  ^  0.  The  condition  reads  then  as  follows:  if  £  is 
even  then  Cf  >  0  for  at  least  one  odd  integer  i,  while  for  £  odd  a,-  >  0  for  at  least  one  even  i  and 
one  odd  i,  0  <  i  <  £. 


^  Department  of  Mathematics  and  Computer  Science,  Ben  Gurion  University  of  the  Negev,  Beer 
Sheva,  Israel,  84105. 
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Knot  selection  for  quadratic  splines 

Jean  Louis  MERRIEN,  Paul  SABLONNIERE,  Rennes,  France 


Resume  :  Soit  =  {  a=x^<  Xj^<...x^=  b  }  une  subdivision  de  I  =  [a,b]  et  S  I'espace  des  splines 

quadratiques  C^  ayant  comme  suite  de  noeuds  X^  ,  ou  X^  =  {  x-  =  d-a^)  Xj  + 

0  <  i  <  n-1 }.  On  s'interesse  au  probleme  d'Hermite  unisolvent  classique  :  trouver  veS  verifiant 
v(xp  =  et  v'(xj )  =  yj  ,  avec  des  donnees  quelconques  {  y.,  y^  ,  0  <  i  <  n  },  ou  au  probleme  de 

Lagrange  associe  quand  les  yj  sont  evaluees  en  fonction  des  donnees  y-.  En  pratique,  on  choisit 

souvent  a-  =  ^  ,  pour  les  noeuds  intermediaires,  mais  ce  choix  n'est  pas  toujours  optimal.  Nous 

tentons  de  definir  des  choix  optimaux  des  a-  associes  a  certains  criteres  definis  sur  v,  par  example 

minimisation  de  certaines  semi-normes,  de  la  courbure,  de  la  longueur  de  la  courbe,  des  normes 
des  operateurs  d’interpolation,  des  constantes  d'erreur.  Un  probleme  lie  est  le  choix  de  X^  donnant 

de  bonnes  approximations  des  yj  lorsque  les  y-  sont  les  seules  donnees.  On  donne  des  exemples 

numeriques  pour  illustrer  les  difFerents  choix  de  noeuds  X^. 

Abstract :  Let  X^  =  {  a=x^<  Xj^<...Xj^=  b  }  be  some  partition  of  I  =  [a,b]  and  S  be  the  space  of 

quadratic  splines  with  knot  sequence  ,  where  X^  =  {  xj  =  (l-a^  x-  +  ®  <«^<1 , 

0  <  i  <  n-1 }.  We  are  interested  in  the  classical  unisolvent  Hermite  problem  :  find  veS  satisfying 
v(xp  =  y-  and  v'(x. )  =  y- ,  for  arbitrary  data  {  y^,  y-  ,  0  <  i  <  n  }  or  the  associated  Lagrange  problem 

when  the  derivatives  are  not  given,  but  computed  in  some  way  from  the  data  {y-}  .  In  practice,  one 

often  chooses  a-  =  ^  ,  but  this  choice  is  not  always  optimal.  We  try  to  define  "optimal  choices  of 
^  2 

a^s  depending  on  the  criteria  to  be  satisfied  by  v,  e.g.  minimizing  some  semi-norms  on  v,  the 

curvature,  the  length  of  the  curve,  the  norms  of  the  interpolation  operator,  the  error  constants.  A 
related  problem  is  the  choice  of  X^^  giving  good  approximate  values  of  (y-}  when  only  the  data  (y^) 

are  given.  Numerical  examples  are  given  to  illustrate  the  different  knot  selections  X^. 


J.L.  MERRIEN,  P.  SABLONNIERE,  Laboratoire  LANS,  INSA  de  Rennes 
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On  the  spectral  radius  of  the  subdivision  operator 
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and 
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The  spectral  radius  of  a  subdivision  operator  is  fundamental  in 
determining  both  convergence  of  stationary  subdivision  schemes  and 
smoothness  of  the  associated  refinable  function.  In  this  talk  we  present 
some  results  on  the  problem  of  determining  the  spectral  radius  of  a 
subdivision  operator. 


Computation  of  the  Monotone  Cubic  Spline  Interpolant 

by 

Charles  A.  Micchelli  and  Florencio  I.  Utreras 

In  this  work  we  apply  the  theory  developed  in  two  preceding  ■  papers  to 
build  a  computational  scheme  for  computing  the  Best  Interpolating  Monotone 
Cubic  Spline,  i.e..  the  function  u  in  jEZ^"[0,1]  which  is  monotone,  interpolates 
non  decreasing  data,  and  minimizes  the  seminorm  f^[ti"(t)]^dt.  We  present  a 
numerical  example  to  illustrate  the  properties  of  the  method. 


SMOOTHING  NOISY  DATA  BY  LOCAL  KRIGING  WITH  NUGGET  EFFECTS 

by 

Pierre  Montes,  ing.  Ph.D 
Ecole  Polytechnique  de  Montreal 
Departement  de  genie  civil 
Montreal,  Quebec,  Canada 


Abstract  -  This  paper  deals  with  the  problem  of  smoothing  noisy  data. 
Given  N  observed  data  ziu^)  =  P{u^)  +  e{Uj)  ,  i  =  1,  .  .  .,iV  ,  where 
e  Q,  Q  s  M"  ,  e(u_^)  are  independent  normal  random  errors  with  zero 

mean  and  known  or  unknown  variance,  P  is  only  known  to  be  smooth.  We 
want  to  estimate  P  in  a  reproducing  kernel  Hilbert  space  so  that  it 
minimizes  the  positive  functional: 


a>(P,) 


7  =1 


[z{u^)  -  (Uj)  ]  2 


+ 


i=N  j=N 


E  E 

I'sl  T=1 


bj^  bj  K{Uj_,  Uj) 


where  K(.,.)  is  the  reproducing  kernel  (generalized  covariance)  ,  bj^  and 
bj  are  scalars,  X  is  a  positive  number,  and  NX  is  the  nugget  effect 
to  be  chosen  in  such  a  way  that  the  solution  preserves  smoothness 

and  fidelity  to  the  data.  This  problem  was  solved  in  the  70 's  and  there 
exist  softwares  which  are  capable  of  findind  the  optimum  parameter  X 
which  minimizes  the  Wahba's  generalized  cross-validation  function 
GCV(A.)  .  In  this  paper,  I  propose  a  method  of  choosing  a  first  guess  of 
X  .  When  the  variance  of  the  errors  is  unknown,  this  choice  is  based 
upon  the  gross  variogram  of  the  data  at  the  vicinity  of  the  origin. 
Then,  cross-validation  method  is  used  to  find  the  optimum  smoothing 
parameter.  The  method  is  applied  in  the  local  kriging  scheme  developed 
earlier  in  order  to  smooth  scattered  data  on  the  line,  in  the  plane,  in 
the  three-dimensional  space,  or  on  the  unit  sphere. 


Shape  Functions  Construction 
by  Different  Polynomial  Approximation 
in  the  Finite  Element  Method  * 
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Abstract 

In  the  analysis  of  a  Finite  Element  Method  (FEM)  given  an  element, 
we  can  describe  his  shape  by  a  set  of  elementary  functions  known  as  shape 
functions.  The  approches  in  describing  these  functions  are  quite  differ- 
ents.  In  spaces  iPf  or  they  generally  are  a  product  of  lagrangian 
polynomials  if  the  coordinate  system  can  be  chosen  with  the  axes  parallel 
to  the  sides  of  the  element,  otherwise  a  system  of  barycentric  coordinate 
(sometimes  called  area  coordinates)  could  be  introduced. 

The  aim  of  this  paper  is  the  description  and  the  representation  of  shape 
functions  when  the  element  has  triangular  or  tetrahedal  shape  (the  sim¬ 
plest).  The  representation  has  been  done  by  using  two  algorithmic  shemes: 
Neville- Aitken  and  De  Casteljau.  For  these  schemes  we  have  deduced  very 
important  properties. 

A  completely  new  approach  seems  their  description  as  attractor  of  a  suit¬ 
able  Iterated  Functions  System  (IFS).  This  idea  becomes  from  the  fractal 
interpolation.  The  main  benefit  of  this  technique  is  the  compacteness  of 
their  internal  representation  on  a  computer  and  the  possibility  of  a  very 
fast  reconstruction. 

A  comparison  of  different  behaviours  using  known  procedures  is  made  and 
an  analysis  of  the  approximation  errors  is  provided. 


*This  work  has  been  done  with  the  support  of  CNR  within  ” Progetto  Straiegico  Sisiema 
Lagunare  Veneziano”  and  40%  of  MURST  on  the  project  ’’Analisi  Numerica  e  Matematica 
Computazionale” 
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High  Order  Approximation 
by  Parametric  Polynomials 

Knut  Mprken 

Department  of  Informatics,  University  of  Oslo, 
P.O.  Box  1080,  Blindern,  0316  Oslo  3,  Norway 
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Abstract.  In  this  joint  work  with  Tom  Lyche,  we  study  approximation  of  smooth  parametric  curves  by 
polynomial  curves,  i.e.,  each  component  of  the  approximating  curve  is  a  polynomial.  It  is  known  that  in 
special  cases,  the  approximation  order  for  plane  curves  is  2n,  where  n  is  the  polynomial  degree.  Here  we 
attempt  to  generalize  these  results. 


Nonnegativity  Preserving  Interpolation  by 
PoweU-Sabin  Splines 

Bemd  Mulansky 
Technical  University  of  Dresden 


Given  a.  set  D  —  C  of  scattered,  data  sites  and  associated  nonns^ativG  data 
values  F  =  {fi}  C  IR,  the  objective  is  to  find  a  nonnegative  interpolating  function. 

In  order  to  define  the  so-called  Powell-Sabin  splines,  each  triangle  (macrotriangle)  of  a 
suitable,  possibly  data-dependent,  triangulation  of  the  data  sites  is  splitted  into  six  micro- 
triangles.  Every  piecewise  quadratic  function  on  the  refined  triangulation  is  uniquely 
determined  by  its  values  and  gradients  at  the  vertices  of  the  macro-triangulation,  as  can 
easily  shown  using  the  Bernstein-Bezier  representation  of  the  quadratic  pieces  on  the 
microtriangles. 

The  simple  sufficient  nonnegativity  conditions,  i.e.,  the  nonnegativity  of  all  B-ordinates, 
result  in  a  solvable  system  of  linear  inequalities  for  the  gradients,  which  is  separated  with 
respect  to  the  vertices.  Introducing  a  suitable  objective  functional,  the  construction  of 
a  nonnegative  Powell-Sabin  interpolant  reduces  to  the  solution  of  independent  quadratic 
optimization  problems  in  two  variables,  namely  the  partial  derivatives,  for  each  vertex. 
The  proposed  method  is  illustrated  by  a  few  mmierical  examples. 
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Progress  in  Data  Dependent  Triangulations 
Edmond  Nadler 

The  study  of  data  dependent  triangulations,  initiated  by  Dyn,  Levin  and  Rippa  in  a  1990 
paper  in  IMA  J.N.A.,  is  concerned  with  the  choice  of  good  or  best  triangulations  with 
a  given  set  of  data  points  Vi  6  as  vertices,  and  also  dependent  upon  a  vector  /  of 
data  given  at  the  points  Vi-  This  contrasts  with  more  traditional  methods,  where  only  the 
points  Vi  in  the  plane  are  taken  into  account,  and  the  triangulation  of  choice  is  often  the 
well  known  Delaunay  triangulation. 

These  problems  are  looked  at  in  the  setting  of  scattered  data  interpolation  with  piecewise 
polynomials  defined  over  a  triangulation.  Given  a  triangulation  A  :=  of  a  set  of 

distinct  and  noncolinear  data  points  V  :=  let  5J(A)  :=  {s  €  C^(Q)  :  s]?;-  6  Vd}, 

with  Q  :=  U^—jTi  and  Vd  fbe  space  of  bivariate  polynomials  of  degree  d.  The  scattered 
data  interpolant  is  a  function  s  €  S^(A)  which  interpolates  the  data  vector  /.  I  focus  on 
the  the  simplest  case  of  5°,  the  continous  piecewise  linear  functions,  which  has  plenty  of 
interest,  and  for  which  the  appropriate  data  vector  /  is  simply  a  function  value  at  each 
point  Vi,  and  for  which  there  is  obviously  a  unique  interpolant. 

In  previous  work,  the  data  dependent  triangulations  have  been  obtained  by  minimizing 
different  criteria  measuring  the  roughness  of  the  interpolating  function.  In  the  present 
work,  possibilities  are  explored  for  using  other  types  of  criteria,  such  as  shape  preservation, 
and  ones  based  on  the  results  of  the  author’s  thesis  (1985)  for  the  optimal  shape  of  triangles 
for  local  best  approximation. 
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Convexity  Preserving  Interpolation  of  Scattered  Data 


Marian  Neamtu 
Department  of  Mathematics 
Vanderbilt  University 
Nashville,  TN  37240,  USA 


The  purpose  of  the  talk  is  to  discuss  some  problems  associated  with  smooth  convexity 
preserving  interpolation  of  scattered  data  in  and  R^.  An  interpolation  method  is 
proposed  which  utilizes  convex  parametric  Powell-Sabin  patches. 


THE  ITERATIVE  SOLUTION  OF  A  NONLINEAR 
INVERSE  PROBLEM  FROM  INDUSTRY: 
DESIGN  OF  REFLECTORS 


Andreas  Neubauer 
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Abstract: 

We  report  about  a  project  tbat  has  been  done  for  and  in  coorporation  with  an  Austrian 
company  which  designs  and  manufactures  lightings  for  industrial  illumination  purposes. 
The  lightings  consist  of  a  lamp  with  known  light  distribution  and  a  reflector.  Our  task 
was  to  develop  a  mathematical  model  and  an  algorithm  (resulting  in  computer  software) 
for  the  problem  of  designing  the  shape  of  the  reflector  in  such  a  way  that  the  illumination 
distribution  on  the  unit  sphere  centered  at  the  lamp,  obtained  by  the  reflected  outgoing 
light  rays  which  are  aU  thought  to  emanate  from  the  lamp,  behaves  in  a  prescribed  way. 
In  the  algorithm  for  this  so  called  “inverse  far  field  problem”,  the  reflector  should  be 
described  in  a  way  that  could  be  directly  used  for  CAM  purposes  in  that  company. 

The  derivation  of  the  model  for  the  “inverse  far  field  problem” ,  is  on  the  one  hand  based 
on  the  reflection  condition,  and  on  the  other  hand  on  the  so  called  “balance  condition” . 
This  condition  means  that  the  amount  of  light  leaving  the  lamp  in  any  solid  angle 
increment  dfi,-  should  equal  the  amount  of  light  reflected  into  the  corresponding  angle 
increment  dfl^  in  the  illuminated  sphere. 

For  deriving  a  numerical  algorithm,  we  have  to  describe  the  reflector  by  a  finite  number 
of  variables.  Since  for  CAM  purposes,  the  company  interpolates  by  splines,  we  have 
chosen  a  bicubic  B-spline  representation  for  the  reflectors.  The  coefficients  of  this  B- 
spline  representation  are  determined  in  such  a  way  that  the  balance  condition  holds  best 
possible  in  the  least  squares  sense.  The  resulting  constrained  minimization  problem  is 
solved  iteratively  with  a  projected  conjugate  gradient  method. 

Finally,  we  present  a  numerical  example  showing  that  the  cdgorithm  works  in  practice. 


Restricted  Splines  with  Modified  Moments 
Gerhard  Opfer 

University  of  Hamburg,  Germany 

In  the  last  vears  the  so-called  histo-splines  gained  considerable  attention.  Roughly  speak¬ 
ing  these  splines  are  used  to  approximate  histograms,  which  are  in  mathematical  terms 
step  functions.  Normally,  the  splines  are  constructed  in  such  a  fashion,  that  the  area  of 
the  underlying  step  function  is  stepwise  preserved.  There  are  several  papers  by  Jochen  W. 
Schmidt  and  co-authors  (Dresden)  on  this  subject.  In  many  cases  the  histosplines  repre¬ 
sent  counting  results  which  are  non-negative  or  population  statistics  which  are  increasing. 
Therefore  it  is  natural  to  try  to  include  these  conditions  in  the  construction  of  the  ap¬ 
proximating  splines.  One  can  try  to  construct  splines  of  this  type  directly  by  choosing  a 
certain  class  of  splines. 

We  will  investigate  the  mentioned  problem  a  little  differently  and  more  generally.  We 
shall  treat  the  following 

Problem.  Given  a  mesh  with  N  >  I  mesh  points  on  the  real  line; 

A  :  a  =  <  <1  <  •  ■  •  <  tiV  = 

'■•nd  two  fixed  integers  fc,n  >  0.  In  addition  there  are  iV  —  1  real  numbers  aj,  j  — 
i,2, . . . ,  A’  —  1  and  a  fixed  polynomial  p  €  !!„  given.  With  these  data  we  want  to  find 
fuuccions  X  €  W^[a,b]  (.t  €  C^'[a,6]  is  essentially  the  same)  which  minimize 


subject  to 

(2)  :=  /  p(t)x{t)dt  =  aj,  j  =  1,2,. . .  ,N  -  1. 

Jtj 

The  functional  Xj  will  be  called  modified  n-th  moment  in  [tj,tj+i].  If  we  choose  Pni^)  = 
we  obtain  the  ordinary  moments,  including  the  area.  In  a  &st  step  we  shall  derive  necessary 
conditions  for  the  stated  problem  and  in  a  second  step  we  will  include  non-negativity 
restrictions.  The  latter  problem  will  consist  of  finding  some  additional  mesh  points  with 
certain  conditions. 


Local  Surface  Approximations  from 
Unorganized  Points 

Karsten  Opitz 

Forschungsgruppe  fiir  Scientific  Visualization 
Gesellschaft  fiir  Mathematik  und  Datenverarbeitung 

A  method  for  constructing  local  smooth  approximations  to  a  surface 
from  a  triangulated  set  of  unorganized  points  in  3-space  is  described. 
For  a  given  vertex  of  the  triangulation,  a  local  model  surface  is  deter¬ 
mined  which  interpolates  a  set  of  neighbors  of  the  vertex.  The  necessary 
correspondence  between  a  neighborhood  of  the  given  vertex  and  the  pa¬ 
rameter  plane  is  established  through  an  approximation  to  the  inverse 
exponential  map.  The  application  of  bivariate  scattered  data  methods 
then  yields  a  parametric  representation  of  the  local  model  surface.  The 
model  surface  can  be  used  for  the  estimation  of  geometric  measures  such 
as  curvature  elements  as  well  as  input  for  the  construction  of  triangular 
surface  patches. 
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G'2-continuous  cubic  algebraic  splines  and  their  efficient  display 

Algebraic  splines  are  built  with  arcs  of  algebraic  curves  (i.e.  given  by  implicit  poly¬ 
nomial  equations)  and  have  been  introduced  by  Sederberg  in  the  80  s.  The  spline  is  G  - 
continuous  if  at  each  joint  the  adjacent  arcs  have  the  same  curvature.  A  method  is  described 
to  construct  a  G'-continuous  cubic  algebraic  spline  consisting  of  n  convex  segments  that 
interpolates  n  -h  1  joint  points.  The  full  spline  depends  on  the  interpolation  points  and 
five  homogeneous  parameters.  Each  arc  of  the  spline  fits  a  triangle,  two  of  whose  sides 
and  vertices  are  the  lines  tangent  to  the  arc  at  its  endpoints.  Our  spline  is  characterized 
by  the  fact  that  the  implicit  equation  of  each  of  its  segments,  when  expressed  in  terms  of 
the  barycentric  coordinates  with  respect  to  its  corresponding  triangle,  has  always  the  same 
expression; 

F{s,t,  u)  =  as~u  +  bsit~  —  cst~  —  dFu  +  tstu  =  0. 

vhere  s,t,u  are  the  barycentric  coordinates.  These  triangles  compose  the  control  polvgon 
of  'Gio  spline.  The  spline  may  be  modified  globally  to  follow  more  faithfully  or  fit  more 
tightly  the  control  polygon.  The  five  homogeneous  parameters,  which  determine  the  full 
spline  are  expressed  in  terms  of  shape  handles  that  control  this  behavior.  To  track  the 
spline  we  determine  the  barycentric  coordinates  of  the  points  of  each  segment.  Since  all 
segments  have  the  same  equation,  it  is  enough  to  track  only  one.  The  tracking  algorithm 
uses  the  special  knowledge  one  has  on  the  behavior  of  one  cubic  segment:  assuming  we 
start  at  the  endpoint  with  barycentric  coordinates  {s,t,u)  =  (1,0,0)  then  as  we  march 
along  the  curve  the  t-coordinate  increases  up  to  a  certain  value  t  =  to  and  then  decreases 
monotoneously  to  t  =  0,  till  the  second  endpoint  (0,0, 1)  is  reached.  To  display  the  spline 
the  actual  cartesian  coordinates  for  each  segment  have  to  be  computed. 

Marco  Paluszny.  Dept  de  Matematicas.  Universidad  Central  de  Venezuela.  Aptado 
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B-Spline  Approximation  for  CAD/CAM 

and  Visualization 


N.  M.  Patrikalakis  *  L.  Bardis  ^ 
S.  T.  Tuohy  *  F.  -E.  Wolter  * 


To  facilitate  the  exchange  of  geometric  data  between  different  geometric  modeling  systems, 
the  approximate  conversion  of  high  degree  B-spline  or  procedurally  defined  geometries  to 
lower  degree  representations  is  required.  High  degree  representations  arise  from  fitting 
high  degree  polynomial  functions  to  point  data,  from  intersections,  from  lofting  high  de¬ 
gree  curves,  and  from  trimming  low  degree  surfaces.  Procedural  definitions  are  produced, 
for  example,  by  offsetting,  sweeping,  trimming  or  blending.  In  this  lecture,  we  will  review 
the  algorithms  we  have  developed  to  approximate  high  order  and  procedurally  defined  con¬ 
tinuous  univariate  and  bivariate  functions  with  B-splines  in  an  efficient  manner.  Examples 
for  the  approximation  of  surface  curves,  and  generalized  cylinder,  offset  and  blending  sur- 
illustrate  the  method.  These  approximation  algorithms  are  based  on  an  iterative 
scheme  in  which  the  approximation  error  is  estimated  by  adaptive  sampling  and  is  reduced 
through  adaptive  knot  addition.  We  will  also  discuss  local  bounds  for  the  approximation 
error  as  a  motivation  for  these  algorithms.  For  the  approximation  of  high  order  B-splines, 
the  effectiveness  of  a  global  error  bound  is  discussed. 

The  representation  of  functions  describing  a  measured  physical  property  (e.g.  via  sparsely 
scattered  or  densely  defined  point  data)  requires  methods  for  the  creation  of  B-spline  ap¬ 
proximations  that  capture  inherent  features  of  the  data.  Such  representations,  for  example, 
facilitate  archiving,  data  storage  reduction,  visualization  and  other  more  general  interro¬ 
gation.  Since  the  measuring  devices  or  sensors  used  to  produce  such  functions  are  of  finite 
precision,  there  is  a  need  to  represent  the  function  values  together  with  their  uncertainty 
and  the  uncertainty  in  the  independent  variables.  In  this  lecture,  we  will  discuss  the  con¬ 
cept  of  enveloping  or  interval  B-splines.  This  concept  allows  the  representation  of  both 
data  and  associated  uncertainty.  We  will  present  an  algorithm  for  the  creation  of  interval 
B-splines  based  on  minimization  with  linear  constraints.  We  will  illustrate  the  method 
using  geophysical  ocean  data. 


*MIT,  Department  of  Ocean  Engineering,  Design  Laboratory,  Cambridge,  MA  02139-4307,  USA 
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DIFFERENTIABILITY  CONDITIONS 
OF  SELF-AFFINE  NON-POLYNOMIAL  CURVES 


Ronaldo  Marinho  Persiano 
Antonio  A.  F.  Oliveira 


LCG/COPPE 
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Self-affine  property  is  usually  associated  to  "wild"  curves.  However,  all  parametric 
polynomial  curves  own  the  self-affinity  property  that  is  closely  related  to  the 
subdivision  principle  natural  to  the  De  Casteljau’s  algorithm  [1].  By  an  adequate 
generalization  of  the  De  Casteljau’s  subdivision  scheme  a  bigger  class  of  (possibly 
non-polynomial)  parametric  self-affine  curves  may  be  built  satisfying  the  usual 
properties  of  variation  diminishing  and  invariance  to  affine  transforms  shared  by 
polynomial  curves.  Careful  choices  of  the  parameter  of  the  subdivision  procedure  aUow 
the  definition  of  non-polynomial  curves  of  class  C^  from  a  set  of  at  least  2k+l  points. 
If  less  of  2k-i-l  control  points  is  given  then  the  construction  scheme  under  the  C^ 
condition  reduces  to  the  usual  De  Casteljau's  subdivision  method  producing  a 
polynomial  curve  of  degree  2k. 


Reference 

1.  Prautzch,  H.  and  Micchelli,  C.:  Computing  curves  invariant  under  halving.  CAGD, 
4(1-2):133-140,  1987. 
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Spline  meshes  for  smooth  blending  of  free  form  surfaces 
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Abstract 

In  the  manner  of  B-splines,  but  without  the  regularity  restrictions  on  the  control  mesh,  the  algorithm 
presented  here  smoothes  a  mesh  of  points  into  a  parametrized  surface  of  low  degree.  Admissible  input 
meshes  may  outline  an  open  or  closed  free-form  shape  of  arbitrary  genus.  There  are  no  restrictions  on  the 
number  of  neighbors  to  a  point.  Blending  ratios  control  the  position  of  patches  and  the  change  of  the  surface 
normal.  The  user  may  choose  freely  from  piecewise  cubics  and  bicubics  for  the  surface  parametrization. 


General  form  of  generalised  B-spline  functions 
on  a  simple  quadrilateral  mesh 
after  k-th  succesive  subdivision  of  a  cube 

by 

A.  R.  M.  Piah 

Department  of  Mathematics  and  Computer  Science 
The  University 
Dundee  DDl  4HN 
Scotland. 

Abstract 

Goodman  has  suggested  an  approach  to  generate  geometrically  smooth  closed  surfaces  by 
defining  biquadratic  splines  whose  domains  are  on  the  faces  of  a  cube.  Upon  a  regular  subdivi¬ 
sion  of  a  cube,  a  simple  quadrilateral  mesh  (SQM)  of  24  faces  is  obtained.  Further  subdivision 
will  result  in  a  new  SQM  of  96  faces.  Following  the  idea  of  Goodman,  the  author  has  con¬ 
structed  biquadratic  splines  on  24-faces  SQM  and  piecewise  tensor  product  polynomial  func¬ 
tions  of  degree  at  most  three  on  96-faces  SQM,  both  with  minimal  support.  Even  though  the 
mesh  is  successively  subdivided,  the  same  continuity  conditions  across  mesh  lines  are  always 
kept  fixed. 

Upon  k-th  arbitrary  regular  subdivision  of  a  cube,  a  new  SQM  of  6(2^)  faces  will  be  obtained. 
In  this  paper,  we  will  construct  the  general  forms  of  generalised  B-spline  functions  on  this 
SQM  with  small  support  and  sum  to  1.  These  basis  functions  are  useful  for  designing,  approxi¬ 
mating  and  interpolating  closed  surfaces.  Because  the  same  continuity  conditions  across  mesh 
lines  are  always  kept  fixed,  data  reduction  can  then  be  used  on  closed  surfaces  defined  by  these 
basis  functions. 


Orthogonal  periodic  splines  and  spline  wavelets 


Gerlind  Plonka 

University  of  Rostock,  Germany 


Let  Nm  be  the  cardinal  B-spline  of  order  m  and  Pm  the  V-periodization  of  Nm- 
Then  the  space  S^{Z)  of  all  iNT-periodic  spline  functions  of  order  m  is  the  set  of  all 
linear  combinations  of  shifts  Pm{-  —  j)  {j  =  0, N  —  1). 

Using  Euler-Frobenius  polynomials  we  shall  introduce  an  orthogonal  basis  of  S^(Z). 
With  the  help  of  this  basis  it  can  be  shown  that  there  is  a  close  connection  between 
S^(Z)  and  the  space  of  trigonometric  polynomials  Tm  :=  span{exp(2xij  •  /V);  j  — 

In  contrast  to  other  approaches  we  shall  prefer  a  representation  which  is  based  on 
the  consequent  application  of  the  symbol  of  spline  interpolation.  In  this  way  some 
known  results  of  M.  Kamada,  K.  Toraichi  and  R.  Mory  (J.  Approx.  Theory  55 
(1988),  27  -  34)  and  S.L.  Lee  and  W.S.  Tang  (Proc.  Edinburgh  Math.  Soc.  (2)  34 
(1991),  363  -  382)  can  be  simplified  and  improved. 

Furthermore,  our  results  can  be  applied  to  cardinal  spline  wavelets  and  iV-periodic 
spline  wavelets  obtained  by  iV-periodization. 

We  shall  give  best  possible  explicit  lower  and  upper  Riesz  bounds  for  the  spline 
basis  as  well  as  for  the  wavelet  basis.  These  bounds  are  very  important  for  finding 
numerically  stable  algorithms  for  the  computation  of  spline  wavelets.  The  properties 
of  the  Euler-Frobenius  polynomials  will  play  an  important  role  in  our  procedure. 
Our  results  improve  known  results  by  C.K.  Chui  („An  Introduction  to  Wavelets", 
Academic  Press,  1992). 


REGULARITY  ANALYSIS  OF  NON  UNIFORM  DATA 
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ABSTRACT 

A  wavelet  which  is  the  derivative  of  a  smoothing  function  is  well  fitted  to  characterize 
function  singularities  of  order  <k  through  multiscale  wavelet  analysis.  A  particular  class  of 
wavelet,  derivatives  of  B-splines,  leads  to  fast  and  efficient  algorithms  for- contours  detection 
of  1-D  and  2-D  signals  by  a  discrete  dyadic  finite  scale  analysis  (Mallat  and  Zhong). 

For  this  class,  the  scale  function  0  is  a  B-spline  of  order  (a-k)  if  the  wavelet  \ir  is  the  k^ 
derivative  of  a  B-spline  of  order  a.  We  express  the  multiscale  analysis  algorithm  in  terms  of  B- 
splines  derivation  formulas  for  the  wavelet  transform  and  by  using  a  knots  insertion  scheme  for 
the  change  of  scale.  This  formulation  allows  us  to  generalize  the  algorithm  for  non  uniformly 
sampled  data.  Under  some  assumptions  on  spacing,  we  obtain  convergent  subdivision 
schemes.  This  algorithm  is  used  to  detect  the  singularities  of  a  function. 

RESUME 

L’analyse  multiechelle  en  transformee  en  ondelettes  est  bien  adaptee  a  la  caracterisation  des 
singularites  d’ordre  <k  d’une  fonction  si  Ton  prend  comme  ondelette  \|/  la  derivee  d  une 
fonction  de  lissage  9.  Une  classe  particulere  d’ondelette,  derivees  de  B-splines,  conduit  a  un 
algorithme  rapide  et  efficace  pour  la  detection  des  singularites  dans  des  signaux  1-D  ou  2-D 
par  une  analyse  en  ondelette  dyadique  discrete  (Mallat  et  Zhong). 

Pour  cette  classe  particuliere,  la  fonction  d’echelle  <|)  est  une  B-spline  d’ordre  (a-k)  si 
I’ondelette  \|/  est  une  spline  derivee  k^®™®  d’une  B-spline  d’ordre  a.  On  formalise  Falgorithme 
d’analyse  multiechelle  en  exprimant  la  transformee  en  ondelettes  a  I’aide  des  formules  de 
derivation  des  B-splines  et  le  passage  d’une  echelle  a  I’autre  par  un  schema  d’insertion  de 
noeuds.  Cette  formulation  permet  de  generaliser  Falgorithme  au  cas  des  donnees 
echantillonnees  irregulierement.  Pour  une  certaine  regularite  du  pas,  on  obtient  des  schemas 
de  subdivision  convergents.  On  utilise  ces  resultats  pour  detecter  les  singularites  d  une 
fonction  dans  le  cas  de  donnees  non  equidistantes. 

Reference  :  S.  Mallat  and  S.  Zhong,  “Characterization  of  Signals  from  Multiscale  Edges”, 
IEEE  T.  on  PAMI,  July  92,  vol  14  n°7  pp7 10-732. 


Special  Functional  Forms  for  Rational  Curves  and  Surfaces 
Helmut  Pottmann,  Technical  University  of  Vienna 


An  important  part  of  geometric  design  is  the  use  of  the  degrees  of  freedom  provided  by 
NURBS  for  achieving  special  functional  properties  or  advantages  in  manufacturing.  We 
present  solutions  to  several  problems  arising  in  applications. 

For  the  mathematical  description  of  milling  procedures  and  the  representation  of  thick 
plates,  offset  curves  and  surfaces  are  important.  However,  offsets  of  rational  curves  or 
surfaces  are  in  general  not  rational  and  therefore  need  to  be  approximated.  In  an  attempt 
to  avoid  this  approximation  which  is  particularly  subtle  for  surfaces,  we  derive  an  explicit 
representation  of  all  rational  curves  and  surfaces  with  rational  offsets.  It  based  on  the  du^ 
representation  and  on  rational  parameterizations  of  the  circle  and  the  sphere.  There  is 
still  enough  flexibility  in  the  obtained  class  of  curves  and  surfaces  to  make  these  elements 
"andidates  for  use  in  free-form  design. 

As  an  example  for  the  use  of  rational  curves  with  rational  offsets  we  show  an  algorithm 
for  the  construction  of  NURBS  representations  of  curves  of  constant  width.  Such  curves 
axe  desirable  profiles  in  certain  cam  mechanisms.  Special  rational  surfaces  with  rational 
offsets  such  as  tubular  surfaces  are  useful  as  blending  surfaces. 

In  several  applications,  the  manufacturing  process  requires  developable  surfaces.  Unlike 
previous  approaches,  we  consider  developable  rational  surfaces  in  their  dual  representa¬ 
tion.  This  allows  us  to  apply  curve  methods.  Particularly,  an  algorithm  for  connecting 
two  curves  with  a  developable  rational  surface  is  presented.  It  is  well-known  that  even 
if  both  curves  are  rational  the  connecting  surface  is  in  general  not  rational  and  therefore 
this  is  an  approximation  problem. 
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Trigonometric  Interpolation  and  Wavelet  Decompositions 
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The  aim  of  this  talk  is  to  describe  decomposition  and  reconstruction 
algorithms  for  nested  spaces  of  trigonometric  polynomials  in  a  suitable  ma¬ 
trix/ vector  form.  The  scaling  functions  of  these  spaces  are  defined  as  funda¬ 
mental  polynomials  of  Lagrange  interpolation.  The  interpolatory  conditions 
and  the  construction  of  dual  functions  are  crucial  for  the  approach  presented 
here.  It  turns  out  that  the  use  of  de  la  Vallee  Poussin  kernels  leads  to  the 
desired  results.  The  investigations  are  based  on  the  works  of  A.  A.  Privalov 
on  orthogonal  bases  of  the  space  and  a  paper  by  C.K.  Chui  and  H.N. 
Mhaskar  on  trigonometric  wavelets.  This  strategy  is  especially  important  in 
applications  as  signal  analysis  whenever  a  separation  into  frequency  bands 
or  octaves  -  and  not  just  single  sine/cosine  frequencies  -  is  necessary. 
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BIORTHOGONAL  WAVELET  PACKETS 
FOR  PARALLEL  SOLUTION  OF  INTEGRAL  EQUATIONS 
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ABSTRACT 

It  is  well  known  that  the  numerical  solution  of  an  integral  equation  is  a  very  ill-posed 
problem.  Classical  tools  to  obtain  an  approximated  solution  are  given  by  the  Truncated 
Singular  Value  Decomposition  method  and  the  Tikhonov  regularization  method  [5]. 
Parallel  version  of  both  these  methods  have  been  already  considered  [4], [6]. 

In  this  work  we  look  for  an  approximated  solution  of  integral  equations  belonging  to 
the  spline  space  Yj,  generated  by  the  m*^  order  B-spline  N^(x)  whit  knots  at  2~^k, 
k  e  Z.  This  space  can  be  decomposed  as  the  direct  sum  of  wavelet  packet  spaces 
Un  [3].  As  consequence  of  this  decomposition  we  obtain  a  block  partitioning  of  the 
matrix  of  the  linear  system  arising  after  discretization  of  the  integral  operator  [1],[2]. 
Based  on  this  special  structure,  we  have  realized  a  parallel  algorithm,  weU  suited  for 
distributed  memory  multiprocessor,  which  achieves  very  good  efficiency  for  the  low  cost 
of  the  communication  and  syncronization  overhead.  Numerical  results  are  presented. 
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allelo,  (1991). 
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multiprocessor.  Parallel  Computing:  Problems,  Methods  and  Applications,  P.  Messina 
and  A.  Murh  (Eds.)  Elsevier  (1992). 
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Spline  wavelets  on  a  bounded  interval  -  Part  I 


Ewald  Quak*  and  Norman  Weyrich 
Center  for  Approximation  Theory 
Texas  A&M  University 
College  Station,  Texas 


In  two  talks,  some  detailed  results  will  be  presented  concerning  a  mul¬ 
tiresolution  analysis  of  spline  wavelets  which  is  not  defined  on  the  whole 
real  axis  but  only  on  a  bounded  interval. 

In  Part  I,  the  basic  ideas  will  be  addressed,  including  a  comparison 
to  other  approaches  for  a  bounded  interval  such  as  the  ones  by  Meyer 
and  Cohen-Daubechies-Jawerth.  Special  boundary  scaling  functions  and 
wavelets  are  constructed  using  B-splines  with  multiple  knots.  To  this  end, 
the  so-called  Chui-Wang  B- wavelet  approach  is  suitably  modified.  As  spline 
wavelets  are  only  semi-orthogonal,  it  is  also  necessary  to  consider  the  cor¬ 
responding  dual  scaling  functions  and  wavelets. 

The  construction  of  the  decomposition  and  reconstruction  sequences 
and  related  algorithmic  and  implementational  questions  will  be  adressed  in 
the  second  part  of  this  presentation. 
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Spline  wavelets  on  a  bounded  interval  -  Part  II 


Ewald  Quak  and  Norman  Weyrich* 
Center  for  Approximation  Theory 
Texas  A&M  University 
College  Station,  Texas 


Using  the  tools  described  in  “Spline  wavelets  on  a  bounded  interval  - 
Part  I”  we  compute  cubic  spline  wavelets  on  the  interval  [0, 1]  as  well  as  the 
dual  scaling  functions  and  the  dual  wavelets.  We  describe  the  usual  decom¬ 
position  and  reconstruction  algorithms.  Additionally,  we  give  an  improved 
version  of  the  decomposition  algorithm  which  works  in  0{j)  operations  on 
each  level  j  and  without  knowing  the  explicit  formulas  of  the  dual  functions. 
Numerical  examples  show  that  the  boundary  effects  in  the  decomposition 
can  be  significantly  reduced  by  using  special  boundary  wavelets  instead  of 
^  - •"-'ating  the  wavelets  defined  on  the  whole  real  axis. 


Which  criteria  to  evaluate  the  quality  of  a  quasi- interpolant  ? 


Cbristophe  Rabut 


This  talk  is  a  reflexion  on  the  main  properties  we  demand  to  a  quasi-interpolant;  we 
give  and  compare  the  interest  of  the  most  used  criteria,  and  explain  the  importance  of 
the  “interpolating  power  of  ip  ”  to  evaluate  how  fax  from  the  data  the  function 
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We  propose  to  say  that  a  function  is  a  quasi-interpolant  if  and  only  if 

a.  p  is  Fi-reproducing,  i.e.  VpGFi  ,  E  P(i)  ^(*  -  j)  =  P 

b.  ||1  -  p\\2  <  1 

(Ffc  is  the  set  of  all  polynomials  of  degree  at  most  fc,  and  p  is  the  “transfer  function”  (or 
the  “symbol  function”)  of  the  function  p,  i.e.  the  function  defined  by 

p{uj)=  E  ■ 


Furthermore,  we  say  that  the  quality  of  a  quasi-interpolant  p  may  be  evaluated  by 
the  following  four  criteria  (in  the  order  of  importance)  : 

i.  II 1  -  ^11 2  is  low  (the  closer  to  zero  ||1  —  ^||2  is,  the  better  the  quasi-interpolant 
is). 

ii.  ,  p{u))  >  0  (^  is  the  Fourier  transform  of  </?  ). 

iii.  p  is  Fn-reproducing  for  high  value  of  n  (the  bigger  is  n,  the  better  the  quasi- 
interpolant  is). 

iv.  p  is  F„' -interpolating  for  high  value  of  n'  (the  bigger  is  n\  the  better  the 

quasi-interpolant  is)  (a  function  p  is  said  to  be  F„' -interpolating  if  and  only 
if  VpeP„.  .  V^:6ZZ^  E  =  ?(*:))■ 
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This  paper  examines  a  problem  arising  in  Numerical  Control.  The  case  studied  is  specific:  a 
toroidal  cutter,  with  vertical  axis,  moves  in  simultaneous  contact  with  a  horizontal  cylinder 
and  with  a  vertical  plane.  The  cutter  climbs  the  cylinder  until  it  reaches  the  ridge,  at  which 
point  it  must  translate  horizontally  through  a  distance  equal  to  the  major  radius  of  the  torus 
before  descending.  Without  this  translation  gouging  of  the  cylindrical  surface  will  occur.  Two 
geometrical  analysis  of  this  situation  are  given,  using  different  representations  of  the  torus. 
The  first  analysis  provides  physical  contacts,  but  features  a  discontinous  trajectory  at  the  top 
of  the  cylinder.  The  second  analysis  shows  that  the  discontinuity  comes  from  a  bifurcation 
phenomenon.  The  bifurcation  theory  provides  both  the  nature  and  the  shape  of  the  singularity. 
Thanks  to  this  information,  the  complete  and  continous  trajectory  can  be  computed. 


Holder  Regularity  of  Subdivision  Schemes  and  Wavelets 


Olivier  Rioul 


CNET  Palis  B  -  France  Telecom 
CNET/PAB/RPE 
38-40.  rue  du  General  Leclcrc 
92131  Issy-Les-Moulineaux,  France 

Email:  rioulQbalsa . issy . cnet.fr 


Abstract.  We  study  special  functions  obtained  a.s  limits  of  subdivision  schemes  in  one  dimension. 
Convergent  dyadic  subdivision  schemes  are  used  for  designing  regular  or  fractal  curves,  which  yield  compactly 
supported  wavelet  bases  under  additional  conditions  on  the  subdivision  masks.  A  “discrete”  approach  is 
undertaken  to  study  the  existence  and  Holder  regularity  of  limit  functions  of  dyadic  or  “p/ 9-adic”  subdivision 
schemes — the  latter  is  a  flexible  “rationar  generalization  which  has  already  found  application  in  digital  signal 
processing. 

Sharp  regularity  estimates  are  derived,  which  are  proven  to  be  optimal  in  the  dyadic  case;  in  constrast 
with  earlier  ones,  they  can  be  easily  implemented  on  a  computer  and  include,  as  a  special  case,  exact  Sobolev 
regularity  estimates. 

The  main  difficulty  of  the  rational  case,  as  compared  to  the  dyadic  one,  is  the  lack  of  shift  invariance: 
This  results  in  an  infinite  number  of  different  limit  functions  instead  of  shifted  copies  of  a  single  one.  In  this 
case,  a  direct  application  of  former  ideas  becomes  impossible.  However,  our  “discrete”  approach  allows  us 
to  extend  the  results  of  Daubechies  and  Lagarias  (on  Holder  regularity  estimates  based  on  infinite  products 
of  matrices)  to  the  rational  case. 


Riesz  Bases  and  Frames 
for 

Shift-Invariant  and  AfRnely-Invariant  Subspaces  of  L2{'R^) 

Amos  Ron 

University  of  Wisconsin-Madison,  U.S.A 

Frames  for  principal  shift-invariant  and  finitely  generated  shift-invariant  subspaces 
of  L2  are  studied  first.  The  characterization  of  frames  is  done  via  the  decomposition  of 
the  space  into  spectral  “fibers”,  and  imposes  no  smoothness  or  decay  conditions  on  the 
generators.  The  results  draw  interesting  connections  between  frames  and  tight  frames 
on  the  one  hand,  and  the  notions  of  quasi-stable  and  quasi- orthogonal  generating  sets  (as 
introduced  bv  de  Boor,  DeVore  and  Ron  in  [1])  on  the  other  hand.  The  characterization 
is  then  extended  to  infinitely  generated  shift-invariant  frames  (and  Riesz  bases). 

The  results  apply  directly  to  Weyl-Heisenberg  frames,  but,  seems,  however,  to  only 
reproduce  (essentially)  known  results  about  this  relatively  accessible  case. 

The  wavelet  (=:affine)  case  is  approached  then  via  the  idea  of  one-sided  frames.  It  pro¬ 
vides  a  characterization  of  affine  Riesz  bases  (generated  by  one,  finitely  many,  or  countably 
many,  functions)  for  T'2(R^).  We  believe  that  result  to  be  new  even  for  the  case  d  1,  and, 
in  fact,  know  of  no  complete  characterization  of  Riesz  bases  even  for  univariate  wavelet 
Riesz  bases  generated  by  a  single  function.  Necessary  conditions  and  sufficient  conditions 
are  easily  derived  from  the  above  mentioned  characterizations.  The  approach  also  results 
similar  characterization  of  wavelet  frames. 

The  entire  work  is  joint  with  Zuowei  Shen. 

[1]  C.  de  Boor,  R.  DeVore  and  A.  Ron,  The  structure  of  shift  invariant  spaces  and  appli¬ 
cations  to  approximation  theory,  CMS-TSR  92-08,  University  of  Wisconsin-Madison, 
February  1992,  J.  Functional  Anal.,  to  appear. 


Interpolation  with  an  Arc  Length  Constraint 
John  A.  Roulier 

Department  of  Computer  Science  and  Engineering 
University  of  Connecticut 
Storrs,  CT  06269-3155 
E!-mail:  jrou@brc.uconn.edu 

Theory  and  algorithms  will  be  presented  for  the  generation  of  a  parametric  curve 
f  into  or  R?  which  interpolates  a  finite  set  of  data  points  Pq,  Pi,  ■  ■  ■ ,  Pn,  preserves 
the  shape  of  the  data  if  possible,  and  attains  a  specified  arc  length  L.  Generally, 
unit  tangent  vectors  Uq  and  u„  are  specified  at  the  two  end  points  Pq  and  P„.  We 
also  allow  for  specification  of  tangent  directions  at  all  of  the  data  points.  Of  course, 
we  must  have  L  >  I  PkPk+i  \  in  order  for  such  an  f  to  exist  since  the  shortest 

curve  interpolating  the  data  is  the  piecewise  linear  curve  connecting  the  consecutive 
data  points.  Furthermore,  if  L  is  too  large,  the  shape  indicated  by  the  data  may 
not  be  preserved  since  inflections  or  loops  may  be  forced  on  f . 

This  work  extends  and  uses  recent  work  by  Roulier  and  by  Roulier  and  Piper  on 
length  specification  for  the  case  of  two  data  points  and  unit  tangent  vectors.  In  this 
work,  the  arc  length  is  reduced  to  a  function  L(a)  of  a  single  variable  a.  Conditions 
are  given  in  which  L{a)  is  strictly  concave  upward  on  the  real  line  and  tends  to  -|-oo 
as  a  ^  ±oo.  The  goal  then,  is  to  find  a  value  of  a  so  that  L{a)  =  L.  Once  such  an 
a  is  found,  f  is  determined  by  that  a. 

The  form  of  the  function  L{a)  allows  the  use  of  the  secant  method  to  rapidly 
produce  an  acceptable  value  of  a.  An  algorithm  will  be  presented  which  implements 
these  ideas  for  our  case.  No  user  input  of  starting  values  for  the  secant  method  is 
required.  Several  examples  will  be  given. 
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Reconstruction  de  Courbes 

Jean-Christophe  ROUX 

LMC,  Tour  IRMA,  51  rue  des  Mathematiques 
BP  53X,  38041  GRENOBLE  cedex. 
e-mail:  roux@busard.imag.fr 

Etant  donne  un  ensemble  de  N  points  dans  le  plan  que  nous  supposerons  appartenir  a  une  courbe  de  ce 
plan,  nous  proposons  ici  une  methode  de  reconstruction  d’un  ordre  sur  les  points  conforme  a  celui  unduit  par 
un  parametrage  naturel  de  la  courbe,  tel  que  celui  de  I’abscisse  curviligne  par  exemple. 

L’interet  pour  ce  genre  de  probleme  a  ete  suscite  par  de  nombreuses  applications  en  reconnaissance  de 
formes  ou  en  imagerie.  Une  approche  de  reconstruction  basee  sur  des  criteres  de  proximite  a  ete  developpee  par 
Dedieu  et  Favardin  [DF  91]  et  consiste,  pour  un  point  a  une  etape  donnee,  a  rechercher  le  point  suivant  le  plus 
proche  sur  la  courbe. 

Nous  presentons  une  approche  differente  non  basee  sur  une  simple  proximite,  mais  sur  le  fait  que  le  cercle 
osculateur  en  un  point  P  d’une  courbe  permet  d’ordonner  localement  les  points  dans  un  voisinage  de  P  ([B-G]). 

Ainsi,  la  methode  que  nous  employons  dans  le  cas  de  courbes  sans  auto-intersections  procede  en  trois  etapes 
principales,  qui  sont  la  decoupe  de  la  courbe  en  arcs  ouverts  au  moyen  d’un  quadtree,  puis  I’approidmation  par 
un  cercle  de  Fare  contenu  dans  chaque  feuille  de  I’arbre  construit  precedement,  ceci  s’accompagnant  de  la 
projection  des  points  initiaux  sur  le  cercle  caJcule  dans  le  but  de  les  ordonner,  et  enfin  une  derniere  etape 
consistant  a  raccorder  entre  eux  les  differents  arcs  ainsi  calcules  pour  obtenir  le  contour  final  ordonne  de  la 
totalite  de  la  courbe. 

Nous  examinons  done  ici  les  problemes  lies  au  calcul  d’un  cercle  approximant  un  nuage  de  points,  et 
presentons  la  methode  de  minimisation  avec  contraintes  que  nous  avons  retenue,  ainsi  qu’une  interpretation 
geometrique  de  ces  calculs.  Nous  etudions  de  plus  les  proprietes  de  convergence  locale  vers  la  courbe,  ainsi 
que  les  bornes  sur  le  rayon  du  cercle  approximant  en  fonction  des  proprietes  de  courbure  de  la  courbe.  Nous 
examinons  aussi  dans  quelles  limites  I’ordre  des  points  projetes  sur  le  cercle  correspond  a  I’ordre  induit  par  le 
parametrage  sur  la  courbe. 

Nous  exposerons  d ’autre  part  les  choix  algorithmiques  que  nous  avons  faits  pour  resoudre  le  probleme  du 
raccord  des  arcs  de  courbe,  ainsi  que  la  methode  retenue  pour  traiter  le  probleme  des  auto-intersections. 

La  complexite  de  I’algorithme  sera  enfin  etudiee,  ainsi  que  son  comportement  et  sa  robustesse  dans  le  cas 
de  jeux  de  donnees  perturbees. 
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Lagrange  and  Hermite  interpolation  by  bivariate 
quadratic  splines  on  nonuniform  type  2  triangulations 


Paul  SABLONNIERE,  Rennes,  France 
Fatima  JEEAWOCK-ZEDEK,  Lille,  France 


Let  Xj  =  {  X-,  0  <  i  <  2m  }  (resp.  Yj  =  { yj  ,  0  <  j  <  2n  }  )  and  Xq  -  {  X2p  0  <  i  <  m } 

(resp.  Yq  =  { y2j  ,  0  <  j  <  n  }  )  be  two  embedded  partitions  of  I  =  [  ,  X2j^]  (resp.  J  =  [  yQ  ,  y2^] ). 

Let  T  be  the  triangulation  of  fi  =  I  x  J  obtained  by  drawing  diagonals  in  rectangles  whose  set  of 
vertices  is  x  Y^  and  let  .*2  =  X^  x  Yq  .  Let  LS^(n,r)  be  the  subspace  of  quadratic  splines 

whose  partial  derivative  (resp.  D^^S)  is  piecewise  linear  and  continuous  on 

partition  Xq  (resp.  Yq)  of  horizontal  (resp.  vertical)  meshlines  determined  by  Yq  (resp.  Xq).  We 
prove  that  dim  LS2(0,t)  =  3(m+l)  (n+1)  and  that  each  element  S  of  this  space  is  uniquely  defined 
by  {  S(A),  DS(A),  A6j«  }.  With  feC^(O),  we  associate  its  Hermite  interpolant  Hf  on  j«.  By 
replacing  partial  derivatives  by  finite  differences,  we  construct  a  Lagrange  interpolant 
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^f  =  Zf  (A)L^  off,  where  L^  is  a  Lagrange  quadratic  spline  with  compact  support.  For  feC  (fl)  , 
we  compute  constants  in  error  bounds  for  the  norms  of  Hf  -  f  and  ^f  -  f  (and  their  derivatives). 

We  also  compute  the  L°°-  norm  of  £.  All  these  quantities  depending  heavily  on  subpartitions 
Xj  c  Xq  (resp.  Yq),  we  try  to  define  a  reasonable  choice  of  these  subsets  in  terms  of  the 

characteristic  parameters  of  partitions  Xq  and  Yq  . 


P.  SABLONNIERE,  Laboratoire  LANS,  INSA  de  Rennes 
20,  avenue  des  Buttes  de  CoSsmes,  35043  Rennes,  France 
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The  Ubiquitous  Ellipse 
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Abstract 

We  discuss  three  different  affine  invariant  evolution  processes  for  smoothing  planar 
curves.  The  first  one  is  derived  from  a  geometric  heat-type  flow,  both  the  initial  and  the 
smoothed  curves  being  continuous.  The  second  smoothing  process  is  obtained  from  a 
discretization  of  this  affine  heat  equation.  In  this  case  the  curves  are  represented  by  planar 
polygons.  The  third  process  is  based  on  B-spline  approximations.  For  this  process,  the 
initial  curve  is  a  planar  polygon,  and  the  smoothed  curves  are  continuous.  We  show  that, 
in  the  limit,  all  three  affine  invariant  smoothing  processes  collapse  any  initial  curve  into 
an  elliptic  point. 
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Axial  convexity  -  a  well  shaped  shape  property 
Thomas  Sauer 

Mathematical  Institute 
University  of  Erlangen-Nuremberg 
Bismarckstr.  1  1/2,  8520  Erlangen 
sauer@  mi .  uni-erlangen  .de 

Motivated  by  the  non-preservation  of  convexity  by  multivariate  Bernstein 
polynomials,  the  notion  of  axial  convexity  was  introduced  in  order  to  have  a  reasonable 
generalization  of  classical  convexity  to  the  multivariate  setting  that  is  more  compatible 
with  Bernstein  polynomials  and  Bezier  patches,  respectively. 

Indeed,  axial  convexity,  defined  by  plain  geometric  arguments,  is  perceptable  by 
the  human  eye  and  in  this  sense  is  called  a  "shape  property  .  It  provides  various 
properties  that  are  well-known  and  appreciated  for  the  case  of  classical  convexity, 
hence,  it  seems  to  be  an  appropriate  generalization  of  convexity  to  the  multivariate  case. 
More  precisely,  one  can  show  that,  if  resonably  defined,  axial  convexity  of.Bezier  nets 
is  preserved  by  degree  elevation  and  causes  the  resulting  Bezier  patch  to  be  axially 
convex,  too;  moreover,  axially  convex  functions  have  monotonous  convergence  of  the 
associated  Bernstein  polynomials  and  satisfy  a  Jensen  inequality  as  well  as  a  maximum 
principle. 

On  the  other  hand,  this  notion  enables  the  design  of  surfaces  that  e.g.  contain 
saddle-points  (and  hence  are  "visibly"  non-convex)  having  nevertheless  a  describable 
convex  behavior. 
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PARTIAL  DIFFERENTIAL  EQUATION  TECHNIQUES  FOR 
3D  SURFACE  DETERMINATION 

C.  Sbert  et  al. 

Address:  Dept.  Matematiques  i  Informatica.  Univ.  de  les  Hies  Balears, 
07071  Palma  (Balears)  Spain. 

Abstract: 

Current  methods  for  3D-surface  reconstruction  of  medical  images  (MNR, 
CAT  . . . )  are  essentially  bidimensional;  segmentation  is  performed  for  each 
section  and  the  surface  is  generated  and  visualized  by  a  spline-model. 

In  this  paper  we  present  two  direct  methods  to  detect  the  surface: 

1)  We  present  a  model  of  anisotropic  diffussion  equation  proposed  by 
Alvarez-Lions-Morel,  which  detect  the  object  and  by  classical  thersholding  we 
detect  their  contour. 

2)  The  second  one  is  a  new  model  of  active  contour  (snakes)  based  on  a 
partial  differential  equation  proposed  by  Caselles-Catte-Coll-Dibos.  ^From  this 
we  obtain  the  surface  of  the  object  as  limit  of  the  evolution  of  an  initial  surface 
containing  the  object. 

To  visualize  the  detected  surface  we  construct  a  model  wire-frame  from  the 
contours  of  the  object  and  render  it  with  the  classical  techniques  of  computer 
graphics. 

We  present  results  on  MNR  images  of  brain  tumours. 


^50 


Polynomial  Reproduction  Properties 
of  Radial  Basis  Function  Interpolants 


Robert  Schaback,  Gottingen 


For  radial  basis  function  interpolation  of  scattered  data  in  a  bounded  region  of 
d-space,  the  approximative  reproduction  of  high-degree  polynomials  is  studied. 
On  the  theoretical  side,  error  estimates  and  convergence  orders  are  proven,  while 
on  the  practical  side  some  numerical  experiments  are  reported.  The  results 
roughly  imply  that  all  currently  available  radial  basis  functions  will  reproduce 
polynomials  of  arbitrary  d  egree  with  arbitrary  precision  on  arbitrary  compact 
sets,  provided  that  sufficiently  many  well-distributed  centers  are  used.  H  owever, 
the  attainable  uniform  convergence  orders  on  compact  sets  will  strongly  depend 
on  the  given  radial  basis  function  and  the  degree  of  the  given  polynomial. 


Lifting  of  Curves  and  Surfaces  in  CAGD 

Markus  Schichtel 


Center  for  Practical  Mathematics 
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The  lifting  of  curves  and  surfaces  has  become  of  great  interest  for  instance  in  the  au¬ 
tomobile  industry.  Car  bodies  are  designed  using  B-Spline  curves  and  surfaces  today. 
Manipulating  control  points  is  a  widely  used  tool  to  alter  the  shape  of  a  B-Spline  surface. 
However  in  many  applications  the  designer  would  rather  like  to  manipulate  a  point  on  the 
surface  itself.  A  direct  manipulation  could  be  the  displacement  of  a  single  point  on  the 
surface  into  the  ambient  space.  After  the  designer  has  defined  such  a  displacement,  he 
wants  the  original  surface  to  follow  this  displacement.  This  operation  is  called  LIFTING. 

The  first  part  of  the  talk  will  present  solutions  for  the  lifting  of  a  single  curve  resp.  surface, 
where  it  is  assumed  that  continuity  should  be  retained  at  the  boundary.  The  question 
of  how  to  preserve  the  shape  of  a  curve  resp.  surface  during  the  lifting  process  will  also  be 
addressed.  In  this  context  shape  preservation  is  understood  in  the  sense  that  for  instance 
convex  curves  resp.  surfaces  remain  convex  after  lifting  is  carried  out. 

The  second  part  of  the  talk  is  dedicated  to  the  case,  where  the  effect  of  lifting  a  single  sur¬ 
face  stretches  out  to  neighboring  surfaces.  In  this  case  lifting  will  introduce  gaps  between 
adjacent  patches.  The  treatment  of  how  to  close  these  gaps  includes  also  trimmed  surfaces, 
which  are  most  common  in  the  design  of  a  car  body.  The  talk  gives  an  overview  on  the 
different  types  of  boundary  adjustments  and  the  methods  associated  with  them.  Basically 
adjusting  a  parameterline  boundary  to  a  given  curve  in  3-D  space  must  be  treated  entirely 
differently  than  adjusting  a  trimming  curve  boundary.  The  solutions  are  in  most  cases 
approximate  ones.  One  reason  for  this  is  the  fact  that  representing  a  trimming  curve 
exactly  would  require  a  very  high  polynomial  degree. 

The  talk  concludes  with  some  examples  discussing  the  several  types  of  boundary  adjust¬ 
ments. 
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Variation  Diminution  and  Blossoming  for 


Curves  and  Surfaces 

Gerd  Schmeltz 

The  variation  diminishing  property  of  some  curve  representation  schemes 
has  proven  useful  in  many  ways.  This  talk  explores  an  idea  to  extend  this 
notion  to  surfaces.  While  it  is  known  that  direct  translations  of  the  geometric 
definition  result  in  properties  not  shared  by  Bezier  surfaces  (which  seems 
inappropriate),  another  approach  is  used  here:  First,  generalized  blossoms 
are  defined  and  shown  to  exist  iff  the  underlying  curve  representation  scheme 
is  variation  diminishing.  Then  analogous  blossom  structures  for  surfaces  are 
formulated  and  investigated.  They  turn  out  to  be  far  more  more  restrictive 
there  than  in  the  case  of  curves.  It  is  shown  that  only  rational  Bezier  triangles 
with  positive  weights  have  this  property. 

Gerd  Schmeltz 
Vanderbilt  University 


Nashville,  TN  37209 


Interpolation  with  C^-Splines  under  Positivity  Constraints 
Jochen  W.  Schmidt,  Dresden  University  of  Technology 

There  are  several  problems  of  practical  interest  which  lead  to  the  interpolation 
of  data  sets  under  positivity  constraints,  and  the  interpolating  functions  should 
be  of  C^-continuity  at  least.  For  positive  interpolation  polynomials  are  unsuited 
while  the  use  of  spline  functions  is  much  more  promising.  These,  of  course, 
should  be  as  simple  as  possible.  As  it  is  well-known,  positive  interpolation  with 
quadratic  splines  is  not  always  successful  if  the  spline  knots  are  equal  to  the  data 
abscissas.  There  exist  data  sets  with  positive  ordinates  such  that  all  quadratic 
spline  interpolants  axe  not  nonnegative  everywhere.  But  we  can  offer  particular 
rational  quadratic  C^-splines  such  that  positivity  is  always  preserved  provided 
the  occuring  rationality  parameters  are  chosen  accordingly.  This  approach  is 
extended  to  the  positive  interpolation  of  gridded  data  in  2D  as  well  as  in  3D.  The 
rational  bi-  and  tri-quadratic  splines  used  then  are  tensor  products  of  the  splines 
in  ID  mentioned  before. 
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Approximation  with  Helix  Splines 

GeraJd  Seemann 

Technische  Hochschule  Darmstadt 
Fachbereich  Mathematik 
SchlofigartenstraBe  7 
6100  Darmstadt,  Germany 

Helix  splines  are  spline  curves  with  (piecewise)  constant  curvature  and  constant  torsion. 
They  are  of  interest  for  milling  machines.  A  helix  describes  a  simple  motion,  so  helix  splines 
can  be  useful  for  kinematics. 

We  approximate  a  ^ven  set  of  points  P*  in  by  a  helix  spline  with  n  segments  minimizing 
a  least  squares  objective  function.  This  nonlinear  problem  is  solved  with  nonlinear  opti¬ 
mization  programs.  In  every  step  of  an  optimization  program  a  helix  spline  is  constructed 
from  the  unkowns.  We  discuss  diflFerent  approaches  to  define  the  unknowns. 

To  run  the  optimization  program  an  initial  guess  b  necessary.  We  calculate  it  by  using  a 
heuristic  with  differential  geometric  arguments  which  gives  also  an  estimate  for  the  needed 
number  of  segments  n. 

Since  a  helix  is  a  transcendental  curve  it  cannot  be  described  by  polynomial  or  rational 
curves.  To  deal  with  the  problem  in  systems  based  on  rational  polynomials  we  introduce 
an  approximating  rational  Bezier  spline  with  small  error  (including  small  error  in  curvature 
and  torsion)  from  which  the  exact  data  of  the  helix,  i.e.  axis  direction,  opening  angle  and 
pitch,  can  be  reconstructed  without  error. 
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A  Solution  to  the  Subdivision  Problem 
for  Triangular  B- Splines 

Hans-Peter  Seidel 
Department  of  Computer  Science 
Universitat  Erlangen,  Germany 


Abstract 

A  new  triangular  B-spline  scheme  has  recently  been  developed  by  Dahnien,  Mic- 
chelli  and  Seidel.  The  scheme  is  based  on  blending  functions  and  control  points  and 
allows  to  model  C"“^-continuous  piecewise  polynomial  surfaces  of  degree  n  over 
arbitrary  triangulations  of  the  parameter  plane.  The  surface  scheme  exhibits  both 
affine  invariance,  and  the  convex  hull  property,  and  the  control  points  can  be  used 
to  manipulate  the  shape  of  the  surface  locally  [1].  Any  piecewise  polynomial  can  be 
represented  by  the  new  scheme  [2].  A  first  test-implementation  of  the  new  scheme 
has  succeeded  in  demonstrating  the  practical  feasibility  of  its  fundamental  algo¬ 
rithms:  Quadratic  and  cubic  surfaces  over  arbitrary  triangulations  can  be  edited 
and  rendered  in  real-time.  Applications  like  the  filling  of  polygonal  holes,  or  the 
construction  of  smooth  blends,  demonstrate  the  potential  of  the  new  scheme  when 
deahng  with  concrete  design  problems  [3]. 

In  this  talk  we  provide  a  solution  to  the  following  subdivisioa  problem  for  tri¬ 
angular  B-splines:  Let  T  be  a  triangulation  of  IR'^  and  let  T  be  an  arbitrary  subtri¬ 
angulation.  Is  it  possible  to  represent  any  triangular  B-spline  surface  F  over  T  as  a 
B-sphne  surface  over  T? 

VVe  show  that  the  answer  to  this  problem  is  always  positive  (independent  of  the 
specific  choice  of  knots  assigned  to  every  vertex  (knot  clouds)),  and  provide  explicit 
formulas  for  the  new  control  points  cj  of  F  over  T  in  terms  of  the  old  control  points 
cj  of  F  over  T. 

The  resulting  algorithm  generalizes  the  well-known  knot  insertion  algorithm  for 
univariate  B-splines  to  the  multivariate  setting. 


[1]  W.  Dahmen,  C.A.  Micchelli,  and  H.-P.  Seidel.  Blossoming  begets  B-splines  built 

better  by  B-patches.  Math.  Comp.,  .59:97-115,  1992. 

[2]  H.-P.  Seidel.  Representing  piecewise  polynomials  as  linear  combinations  of  mul¬ 

tivariate  B-splines.  In  T.  Lyche  and  L.  L.  Schumaker,  editors.  Curves  and 
Surfaces,  pages  559-566.  Academic  Press,  1992. 

[3]  P.  Fong  and  H.-P.  Seidel.  An  implementation  of  multivariate  B-spline  surfaces 

over  arbitrary  triangulations.  Computer-Aided  Geom.  De.sign,  10(1),  1993. 
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Modeling  of  geological  surfaces  using  finite  elements 
Torn  Springer,  GeoForschungsZentrum  Potsdam 

Many  geological  structures  are  described  by  surfaces  which  represent  the  boundaries  between 
different  geological  units.  These  geological  surfaces  are  known  very  incompletely  and  must  therefore 
be  reconstructed  by  use  of  interpolation  methods. 

In  regions  of  undisturbed  sedimentation  almost  horizontal  surfaces  develop  which  may  be  described 
by  a  continuous  function  z  =  u(x,y),  where  z  is  the  vertical  coordinate  (depth).  However,  tectonic 
stresses  cause  discontinuities  (faults)  and  more  complicated  surfaces  (folds) . 

Modeling  of  general  geological  surfaces  leads  to  the  following  mathematical  problems: 

1.  From  a  set  of  given  points  i  =  l,...,n,  determine  a  function  u(x,y)  with 

u(Xj,yi)=Zj,  i  =  l,...,n.  The  resulting  surface  should  be  "as  smooth  as  possible"  to  show  not  more 
features  than  are  included  in  the  data.  This  is  the  classical  two-dimensional  interpolation  problem. 
For  application  on  geological  surfaces  it  should  be  possible  to  prescribe  additionally  the  slope  of 
the  surface  (i.e.  the  partial  derivatives  u^  and  Uy)  at  some  or  all  of  the  data  points. 

2.  The  same  interpolation  problem  is  to  be  solved  if  discontinuities  are  allowed  along  given  lines  in 
the  xy  plane.  The  allowed  discontinuities  usually  lead  to  a  smoother  surface  in  other  regions. 

3.  In  the  presence  of  folds  a  one- valued  function  u(x,y)  may  no  longer  suffice,  a  parameter 
representation  of  the  surface  must  be  used  instead  and  the  interpolation  algorithm  must  be 
generalized. 

To  solve  these  problems  we  propose  an  algorithm  using  triangular  finite  elements.  In  problem  1  the 
triangle  comers  are  identical  with  the  data  points,  in  case  2  additional  points  along  the  discontinuities 
are  introduced.  Among  all  functions  satisfying  the  given  constraints  (u  and  possibly  u^  and  Uy  at  the 
data  points)  the  sought  function  shall  minimize  an  integral  (potential  energy  of  a  thin  elastic  plate) . 
For  modeling  a  general  surface  in  space  (problem  3)  beside  the  data  points  a  connecting  triangular 
mesh  and  normal  vectors  at  each  data  point  (computable  from  neighbour  points)  are  necessary.  Over 
each  triangle  the  interpolation  is  done  again  by  a  finite  element.  To  ensure  a  continuous  and  smooth 
surface  over  all  triangles  the  algorithm  is  properly  modified. 


Some  Quantitative  Aspects  of  Wavelet  Bases 
Joachim  Stockier 

Universitat  Duisburg^  Fachbereich  Mathematik 
Lotharstrasse  65,  4^00  Duisburg  1,  Germany 

Different  wavelet  bases  have  been  constructed  for  one-dimensional  and  multi-dimensional 
X^-spaces.  One  typical  requirement  for  such  bases  is  the  fast  decay  of  the  basis  functions 
and  of  their  Fourier  transforms.  Usually  this  is  described  in  terms  of  polynomial  or  expo¬ 
nential  decay.  The  decay  properties  are  essential  in  order  to  obtain  local  information  about 
the  X" -function  f  by  its  wavelet  coefficients  at  the  same  time  as  local  information  about 
its  Fourier  transform  /. 

A  finer  measure  for  the  localization  properties  is  the  so-called  time-frequency  window 
associated  with  the  wavelet  function  tp.  In  this  talk  we  give  a  quantitative  analysis  of 
the  time-frequency  window  of  univariate  wavelets.  The  wavelets  under  consideration  are 
defined  by  a  positive  definite  scaling  function  which  generates  a  multiresolution  analysis. 
The  method  employed  by  us  gives  sharp  bounds  for  orthogonal  spline- wavelets  of  even 
order.  The  method  also  applies  to  “non-stationary”  wavelets  which  can  be  obtained  from 
a  multiresolution  of  radial  functions  such  as  Hardy’s  multiquadrics. 


Periodic  Spline— Wavelets 

Manfred  Tasche 

University  of  Rostock  (Germany) 


In  order  to  analyse  periodic  functions  or  functions  defined  on  a  bounded  in¬ 
terval,  it  is  necessary  to  develop  the  theory  of  periodic  wavelets.  In  this  talk 
we  introduce  iV-periodic  spline-wavelets  analogously  to  the  semiorthogonal 
cardinal  spline-wavelets  due  to  C.K.  Chui  (Introduction  to  Spline-Wavelets, 
Academic  Press  1992).  It  is  our  aim  to  present  a  simple  constructive  ap¬ 
proach  to  periodic  wavelets. 

Let  (pj  {j  =  0,1,...)  be  the  iV-periodization  of  Am(2^-),  where  Nm  de¬ 
notes  the  cardinal  B-spline  of  order  m  >  1.  Then  for  j  =  0, 1, ...  we  define 
following  subspaces  of  the  Hilbert  space  L|/- 

Vj  :=  span{(,(?j(  •  -  2"^A;) ;  j  =  0, . . . ,  2^N  -  1}  ,  V-i  :=  {0}  . 

Then  these  subspaces  form  an  iV-periodic  multiresolution  analysis  related  to 
Nm-  In  particular,  using  Euler-Frobenius  polynomials  we  are  able  to  present 
rational  optimal  Riesz-bounds  for  the  spline  basis  of  Vj.  More  details  on  this 
topic  and  on  the  new  technique  of  proof  will  be  given  in  the  talk  of  G.  Plonka 
(Rostock) . 

We  describe  a  new  connection  between  reconstruction  relations  and  decom¬ 
position  relations  based  on  a  block  circulant  matrix  and  Fourier  transform 
technique.  The  condition  of  this  matrix  is  determined  by  the  optimal  lower 
Riesz-bound  of  the  spline-basis  of  Vq.  One  can  see  that  spline-wavelets  of 
order  m  >  6  leads  to  numerical  instabilities.  Finally,  we  get  the  reconstruc¬ 
tion  algorithm  and  the  decomposition  algorithm.  Applying  FFT  we  present 
some  numerical  tests  with  periodic  spline-wavelets. 

If  we  extend  our  results  to  cardinal  spline- wavelets,  then  we  get  some  im¬ 
provements  of  the  known  approach  of  C.K.  Chui. 


Splines  in  a  Topological  Setting 

Geovan  Tavares  and  Sinesio  Pesco 
Departamento  de  Matematica 
Pontificia  Universidade  Catolica 
Rio  de  Janeiro 
Brasil 

Abstract 

We  define  a  topological  data  structure,  called  Handle- Edge,  and  its 
corresponding  structural  operators,  Morse  Operators,  whicli  can  deal,  in 
the  same  environment,  with  surfaces  embeded  in  n-dimensional  spaces.  The 
whole  philosophy  is  object-oriented. 

One  of  the  main  features  of  this  structure  is  that  it  can  handle  surfaces 
represented  in  several  forms  like  splines,  imphcit,  etc.  in  the  same  setting 
allowing  the  manipulation  and  visualization  of  these  surfaces  in  a  unique 
computational  sheme. 

In  this  talk  we  will  be  describing  how  splines  can  profit  from  a  topolog¬ 
ical  environment  by  extending  the  Handle-Edge  [1]  data  structure  to  this 
case. 

Slides  and  video  will  be  shown. 

[1]  G.  Tavares,  A.  Castelo,  H.  Lopes  -  Handlebody  Representation  for  Sur¬ 
faces  and  Morse  Operators.  In  J.  Warren,  Proc.  on  Curves  and  Surfaces 
in  Computer  Vision  and  Graphics  III.  SPIE  (The  International  Society  for 
Optical  Engineering), 


Designing  of  a  ’’progressive  lens” 


M.  TAZEROUALTI 

LMC-IMAG 

Universite  Joseph  Fourier 
BP  53X,  38041  GRENOBLE  cedex. 

e-mail  :  tazer@busard.imag.fr 


Abstract: 

We  introduce  a  new  method  for  designing  a  ’’progressive  lens”.  This  lens  is 
represented  by  a  sufficiently  regular  surface.  On  this  surface  we  prescribe  conditions 
on  its  principal  curvatures  in  some  regions  (far  vision  region  and  near  vision  region) 
and  other  conditions  on  its  principal  directions  of  curvature  in  other  regions  (  nasal 
region  and  temporal  region).  The  surface  is  described  with  tensor  product  B-splines 
of  degre  four.  For  its  computation,  we  have  to  minimize  a  non-q^uadratic  operator. 

This  minimization  is  then  processed  by  an  iterative  method  for  which  quick 
convergence  is  numerically  tested. 


Approximation  of  Parametric  Surfaces  with  Discontinuities 
by  Discret  Smoothing  D^-splines 

Juan  Jose  Torrens 

Departamento  de  Matematica  Aplicada,  C.P.S.,  Universidad  de  Zaragoza, 

Man'a  de  Luna,  3,  E-50015  Zaragoza,  Spain 

In  this  work,  we  consider  the  problem  of  fitting  discontinuous  parametric  surfaces 
when  the  data  points  are  part  of  the  nodes  of  a  curvilinear  grid.  This  problem  may 
appear  in  some  fields,  such  as  Geophysics  or  CAD,  for  modeling  reverse  faults  or 
objects  with  cuts. 

We  state  the  problem  in  the  following  way:  Let  5  be  a  surface  defined  by  a 
parameterization  f  :  D  We  suppose  that  /  is  discontinuous  (giving  rise  to 

faults  or  holes  in  the  surface)  and  that  /  belongs  to  the  Sobolev  space  R^), 

with  m  >  2.  For  any  n  6  N,  let  P„  be  a  given  subset  of  the  set  of  vertices  of  a 
curve  network  lying  on  S  that  divides  this  surface  into  a  finite  number  of  rectangular 
(eventually  triangular)  patches.  The  objective  is  then  to  construct,  for  any  u  €  N, 
an  approximating  surface  from  the  set 

First,  we  assign  parameter  values  to  the  data  points.  For  doing  'this,  we  fix  a 
suitable  planar  domain  fi,  tesselated  into  a  grid  of  equal  squares,  and  then  we  at- 
:ach,  in  order,  the  points  in  Pn  to  part  of  the  vertices  of  these  squares.  We  also 
associate  the  fault  edges  with  a  subset  $  of  fi.  Next,  we  construct  a  parametric 
finite  element  space  V^,  contained  in  the  Sobolev  space  R^),  in  which  search¬ 

ing  for  the  parameterization  of  an  approximating  surface.  This  space  is  defined 
from  the  Bogner-Fox-Schmit  generic  finite  element  of  class  (7*,  making  some  slight 
modifications  in  the  basis  functions  in  order  to  allow  discontinuities  on  $. 

Then  we  consider  the  following  minimization  problem:  find  cr^  such  that 

(V)  i  ^  V^, 

1  Vu6V;^  Ucrr.)<Uv\ 

Jn  being  the  functional  defined  on  V)f  by 

Ib(^)  +Snblm,n,R3, 

beB 

where  B  stands  for  the  set  of  points  in  0,  attached  to  data  points,  pb  is  the  data 
point  corresponding  to  b  E  B,  e„  is  a  given  real  positive  number,  and  |  .  1^,0, is 
the  Sobolev  semi-norm  of  order  m  in  the  space  R^). 

It  is  easily  shown  that,  under  suitable  hypotheses,  problem  (P)  has  a  unique 
solution,  called  discret  smoothing  D^-spline^  which  is  also  the  unique  solution  of 
an  equivalent  variational  problem.  From  this  one  we  finally  obtain  three  linear 
systems  having  the  same  symmetric,  positive  definite  matrix.  These  systems  yield 
the  coefficients  of  cr„  as  a  linear  combination  of  the  basis  functions  of  V^. 

The  convergence  of  the  method  is  studied  and  some  numerical  and  graphical 
examples  are  also  given. 
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NATURAL  SPLINES 


LEONARDO  TRAVERSONI 

DEPARTAMENTO  DE  INGENIERIA  DE  PROCESOS  E  HIDRAULICA 
DIVISION  DE  CIENCIAS  BASIC  AS  E  INGENIERIA 
UNIVERSIDAD  AUTONOMA  METROPOLITANA 
IZTAPALAPA 
MEXICO  D.F. 

MEXICO 


The  concept  of  Natural  Splines  has  been  recently  developped  starting  from 
the  well  known  idea  of  natural  neighbors  created  by  R.  Sibson  in  the  early 
eighties.  Redefining  Sibson’s  interpolant  using  covering  spheres  and  Bernstein 
polynomials  some  interesting  results  are  found.  The  first  one,  a  clear  diferen- 
tiation  of  the  field  of  definition  of  each  interpolating  function.  The  second,  a 
natural  way  of  upgrading  the  aproximating  polynomials  and  a  good  geometrical 
interpretation  of  such  upgrading.  Finally,  the  redefinition  provides  the  link  with 
spline  theory  creating  Natural  Splines  that  is,  splines  based  on  natural  neigh¬ 
bors  whose  field  of  definition  are  the  same  of  Sibson’s  interpolant.  It  is  shown 
that  the  combination  of  all  this  concepts  becomes  a  useful  tool  in  many  fields, 
and  that  a  similar  reasoning  can  be  used  to  find  finite  element  aplications. 


Multiplication  as  a  General  Operation  for  Splines 


Kenji  Ueda 


Abstract 

Degree  raising  and  knot  insertion  are  useful  tools  in  making  different  B-spline 
curves  and  surfaces  compatible.  Multiplication,  which  is  a  basic  B-spline  func¬ 
tion  operation,  can  make  different  rational  B-splines  compatible  to  have  a  com¬ 
mon  weight  function. 

This  talk  presents  a  blossoming  method  for  obtaining  spline  function  prod¬ 
ucts  and  show  its  equivalence  to  Mqrken’s  algorithm  based  on  discrete  B-splines. 
This  spline  multiplication  generalizes  not  only  degree  raising  but  also  knot  in¬ 
sertion  and  rational  spline  compatibility  operation.  The  method  is  applied  to 
surface  design  problems  demonstrating  its  effect  and  added  flexibility. 


Spline  Finite  Element  Software 
A1  Vermeulen 

Rogue  Wave  Software 
PO  Box  2328 
Corvallis,  OR  97330 
USA 

alv(a)  roguewave.com 
fax  503-757-6650 


Abstract: 

The  spline  field  approximation  method  is  used  to  determine  the  load- 
deflection  behaviour  of  a  structure.  For  concreteness,  we  consider  in 
this  talk  the  non-linear  load-deformation  behaviour  of  a  two 
dimensional  bicycle  wheel  model.  The  deflection  of  each  spoke,  and  the 
deflection  of  the  rim,  are  approximated  by  B-spline  functions;  the 
control  points  for  these  splines  are  determined  by  solving  equations 
obtained  using  the  Rayleigh-Ritz  procedure.  In  order  that  the  spokes 
and  the  rim  are  joined  properly,  the  spoke  deflection  splines  are 
constrained  so  that  the  deflection  of  the  spoke  at  the  rim  is  equal  to  the 
rim  deflection. 

The  software  used  to  implement  this  application  falls  into  two 
catagories:  application  independent  software  objects,  and  application 
dependent  objects.  Most  of  the  application  independent  objects 
represent  concepts  from  CAGD,  such  as  splines,  or  concepts  from 
general  mathematics,  such  as  non-linear  equations.  The  application 
dependent  objects  serve  largely  as  glue  which  holds  the  application 
together  and  allows  the  different  pieces  to  communicate  with  each  other. 
Well  show  how  these  various  software  objects  relate  to  each  other,  and 
discuss  the  potential  for  creating  flexible  code  using  object  oriented 
design. 


WAVELET  ANALYSIS  OF 
REFINEMENT  EQUATION  SOLUTIONS 

Lars  F.  Villemoes 

Department  of  Mathematics 
Royal  Institute  of  Technology 
S-10044  Stockholm,  Sweden 
e-mail:  larsv@math.kth.se 

A  basis  of  smooth  wavelets  can  be  used  to  analyze  the  regularity  of  a  compactly 
supported  solution  to  a  refinement  equation.  As  a  result  of  this  procedure,  it  turns 
out  that  there  is  a  simple  relation  between  the  exact  Besov  smoothness  of  the 
solution  and  the  spectral  radius  of  a  subdivision  operator  associated  with  the  the 
coefficients  of  the  equation.  This  operator  acts  on  all  sequence  spaces  £p(Z),  and 
the  spectral  radius  depends  on  p. 

Special  choices  of  p  lead  to  both  well  known  and  some  new  criteria  for  continuity, 
Holder  continuity,  Sobolev  regularity,  differentiability,  integrability  and  absolute 
continuity.  In  some  cases,  estimates  of  the  fractal  (box  counting)  dimension  of  the 
graph  of  the  refinement  equation  solution  can  be  obtained. 

We  illustrate  how  the  theory  can  be  applied  to  give  an  elementary  proof  of  the 
optimal  continuity  range  of  a  one  parameter  family  of  4-point  dyadic  interpolation 
schemes,  considered  both  by  Dubuc  and  Deslauriers  as  well  as  Dyn,  Gregory  and 
Levin.  For  a  special  choice  of  the  parameter,  it  is  known  that  the  solution  is  almost 
everywhere  twice  differentiable  but  not  two  times  continuously  differentiable.  This 
fact  can  also  be  obtained  by  coarse  estimates  of  the  spectral  radii  mentioned  above. 

Finally,  we  discuss  possible  extensions  of  the  method  to  several  dimensions. 


Spline  conversion  for  present-day  CAD  data  exchange  problems 
Annelieke  Vries-Baayens 
TNO  CAD  Centre 
P.O  box  5073 
2600  GB  Delft 
The  Netherlands 


At  present,  data  exchange  between  CAD  systems  still  causes  problems. 
Most  CAD  users  are  not  interested  what  causes  these  problems  but  like  to 
have  them  solved.  To  improve  the  data  exchange  process  between  a 
specific  group  of  CAD  systems  namely  surface  modellers,  spline  con¬ 
version  is  essential. 

The  TNO  CAD  Centre  is  developing  a  so  called  ’Help  Desk’ 
which  must  enable  companies  to  improve  their  data  exchange  process. 

This  ’Help  Desk’  consists,  amongst  others,  of  a  number  of  spline  con¬ 
version  routines  and  spline  comparison  routines.  An  overview  will  be 
given  of  the  structure  of  this  ’Help  Desk’.  The  spline  conversions  which 
are  essential  for  the  improvement  of  the  data  exchange  process  will  be 
presented.  Further,  a  number  of  problems  with  the  spline  comparison 
methods  will  be  discussed. 


Symmetric  Tchebycheffian  B— spline  schemes 


by  M.G.  Wagner  and  H.  Pottmann 
Institute  of  Geometry 
Technical  University  of  Vienna 


Tchebycheffian  B-spline  schemes  are  a  natural  generalization  of  the  polynomial  B- 
spline  scheme  and  possess  essentially  all  properties  which  are  considered  important  for 
CAGD.  Recently,  the  importance  of  these  free  form  schemes  was  shown  by  Lempel  and 
Seroussi  (cf.  [1]).  In  this  paper  we  use  a  more  geometric  approach.  Based  on  the  theory 
of  extended  Tchebycheffian  spaces  and  the  associated  normal  curves,  one  can  define  a 
blossoming  method  and  a  de  Casteljau  algorithm  for  the  construction  of  the  associated 
free  form  curves  (cf.  [2]). 

We  now  determine  all  Tchebycheffian  spline  segments  with  the  following  symmetry 
property:  Over  each  interval,  the  Tchebycheffian  Bezier  curve  to  the  control  points 
bo  •  • .  is  the  same  as  the  curve  to  b,^  ...  bo-  The  question  is  solved  at  hand  of  the 
corresponding  normal  curve  which  turns  out  to  be  an  afhne  W-curve.  The  underlying 
function  space  is  the  null  space  of  a  self-adjoint  linear  differential  operator.  In  case  of 
order  three  one  gets  cubic  spline  curves,  exponential  spline  segements  in  tension,  and  helix 
spline  segments  associated  with  the  function  space  {1,  f,  sin(t),  cos(t)}  as  the  only  possible 
solutions. 


[1]  Lempel,  A.  and  Seroussi,  G.  (1992),  Systematic  derivation  of  spline  bases.  Com¬ 
puter  Aided  Geometric  Design  9,  349-363. 

[2]  Pottmann,  H.  (1993),  The  geometry  of  Tchebycheffian  splines.  Computer  Aided 
Geometric  Design,  to  appear. 
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Ondelettes  et  Interpolation  de  courbes  et  surfaces 
par  des  fonctions  radiales 

Jules  Waku 

Laboratoire  TIMC-IMAG 

Universite  J.  Fourier  BP  53x  F-38041  Grenoble  (France) 

Resume 

L’interpolation  des  courbes  et  surfaces  par  les  fonctions  radiales  [2,5]  passe  par  la 
resolution  de  systemes  lineaires  tres  denses  et  de  grandes  tailles  si  le  nombre  de  donnees 
a  interpoler  est  grand,  ce  qui  est  souvent  le  cas.  Le  caractere  dense  et  de  grande  taille 
du  systeme  sont  les  principales  difficultes  pour  la  resolution  de  ces  systemes  lineaires,  ce 
qui  limite  I’utilisation  de  cette  methode  d’interpolation  qui  est  reconnue  conune  Tune  des 
meiUeures  [2]. 

Pour  faire  face  a  ces  difficultes,  nous  proposons  de  rendre  creuse  la  matrice  du  systeme 
cn  uciiisant  Talgorithme  du  a  Beylkin  et  al.  [6]  fonde  sur  une  representation  de  la  matrice 
dans  une  base  d’ondelettes  [4,3,7]. 

Dans  cet  article,  nous  posons  le  probleme  general  d’interpolation  des  surfaces  par  des 
fonctions  radiales  [5].  Un  bref  rappel  sera  fait  d’une  part  sur  les  conditions  d’existence  et 
d’lmicite  de  la  solution  [1]  et  d’autre  part  sur  la  formulation  matricielle  du  systeme  lineaire 
correspondant  [2,5].  La  mise  en  oeuvre  de  ralgorithme  de  Beylkin  et  al.  montre  que  I  on 
pent  rendre  la  matrice  du  systeme  tres  creuse  tout  en  conservant  une  bonne  precision 
des  resultats,  pour  differentes  fonctions  radiales  [5]  d’ou  I’acceleration  des  calculs.  Nous 
montrons  :  i)  que  les  ondelettes  splines  [4]  ne  sont  pas  adaptees  a  ce  type  de  probleme  ; 
ii)  que  les  ondelettes  a  support  compact  [3]  donnent  de  meilleurs  resultats  et  enfin  nous 
presentons  des  calculs  obtenus  avec  des  ondelettes  biorthogonales  [7]. 
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ABSTRACT 

Given  a  mesh  in  R®  constructed  from  curves  with  continuitv’  C° ,  the  problem  of  finding  a  surface 
of  continuity  of  some  specified  continui^  C^  which  contains  the  given  mesh  has  a  long  and 
venerable  history  and  possesses  a  variety  of  different  types  of  solutions.  In  the  case  that  n  and  k 
equal  1  or  2  solutions  ha\'e  applications  in  Computer  Aided  Design  situations  in  which  designers 
first  construct  models  constructed  from  curved  meshes. 

Solutions  which  satis^’  the  additional  the  constraint  of  "local  control"  -  namely  that  modification  of 
the  surface  by  slightly  changing  one  of  the  curves  of  the  mesh  should  effect  the  resulting  surface 
only  in  a  local  region  of  the  change  -  have  distinct  benefit  to  the  CAD  community.  Existing 
solutions  to  the  general  prdjlem  either  do  not  satisfy  the  constraint  of  local  control  or  impose 
constraints  on  the  curves  of  the  mesh  -  typically  th^  are  assumed  to  be  Bezier  splines. 

The  solution  to  be  presented  satisfies  the  local  control  condition,  allows  the  curves  of  the  mesh  to 
be  aibhraiy,  and  is  of  relatively  low  degree  in  the  case  of  a  mesh  consisting  of  piecewise 
polynomial  curves.  The  solution  employs  the  notion  first  intrwlucedby  CatmuU  and  Rom  that 
judiciousty  chosen  "control  functions"  may  be  substituted  for  control  vertices.  In  particular,  the 
solution  constructed  may  be  considered  as  a  B-spline  averaging  of  surfaces  which  locally  interpolate 
the  given  mesh. 


Multi-sided  rational  surface  patches 
with  independent  boundary  control 


Joe  Warren 

Rice  University,  Houston 


Various  schemes  have  been  proposed  for  creating  rational  surface 
patches  with  a  variable  number  of  sides.  Although  these  schemes  appear 
to  be  unrelated,  they  all  produce  surface  parameterizations  that  have 
parameter  values  for  which  the  numerator  and  denominator  of  the 
defining  rational  functions  vanish  simultaneously.  These  parameter 
values  are  base  points  of  the  rational  parameterization. 

This  paper  shows  that  any  rational  surface  scheme  permitting  the 
independent  control  of  non-adjacent  boundary  curves  must  produce 
parameterizations  possessing  base  points.  Locating  these  base  points 
allows  for  a  classification  of  possible  multi-sided  rational  surface 
schemes. 


Constrained  surface  fitting  using  Powell-Sabin  splines 


Karin  Willemans*  and  Paul  Dierckx 
Katholieke  Universiteit  Leuven.  Department  of  Computer  Science, 
Celestijnenlaan  200A,  B-3001  Leuven,  Belgium 

We  will  consider  the  problem  of  smoothing  a  discrete  set  of  data  values  corresponding 
to  points  arbitrarily  scattered  in  a  bounded  area.  Hereby  data  reduction  is  regarded  as 
an  important  issue  and  the  approximation  should  satisfy  some  additional  conditions. 

Although  tensor  product  splines  have  already  proven  their  usefulness  while  dealing 
with  smoothing  problems,  in  some  cases  they  are  not  very  appropriate. 

Therefore  we  will  consider  in  this  talk  smoothing  with  Powell-Sabin  splines  and  we 
wiii  particularly  pay  attention  to  the  imposition  of  additional  conditions  such  as  con¬ 
vexity  conditions,  boundary  conditions,  .  . .  We  wiU  comment  on  the  translatioii  of  these 
conditions  into  a  suitable  set  of  constraints,  on  the  special  structure  of  the  systems  that 
originate  from  these  constraints  and  on  the  solution  of  these  systems  taking  account  of 
their  specific  structure.  Our  talk  will  be  illustrated  with  many  practical  examples. 
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A  NEW  MULTIVARIATE  BOX-SPLINE  FORMULATION 
FOR  IMAGE  COMPRESSION  AND  RECONSTRUCTION 


Daniel  Williams 
Howard  University 
Washington,  D.C.  USA  20059 
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Abstract.  In  joint  work  with  Harvey  Diamond  and  Louise  Raphael,  we  present  an 
alternative  approach  to  data  compression  and  reconstruction.  Our  algorithm  begins 
hy  first  calculating  an  intensity  function  of  the  images  by  using  a  finite  sequence 
of  approximations  by  multivariate  box-splines  at  successively  halved  resolutions. 
The  key  is  that  these  approximations  are  cpiasi-interpolants,  which  are  nearly  or¬ 
thogonal  projections.  This  fact  ensures  the  local  dependence  of  the  image  function 
approximation  .  This  contrasts  with  customai’y  spline  and  wavelet  methods  which 
employ  exact  orthogonal  projections. 


Fiecev.dse  Cubic  ]}terpolation  Surfaces 
Zhen-Xiang  Xiong 

(Beijing  University  of  Aeronautics  &  Astronautics,  China) 

V7.e  have  constructed  a  kind  of  multivariate  piecewise  interpolation 
polynomials  of  degree  2n-l  that  is  a  kind  of  generalization  from  the 
univariate  interpolation  spline  of  degree  2n-l.  It  may  also  be  called 
’interpolation  spline’,  but  is  different  from  that  defined  by  Prof, 
de  Boor.  We  have  discussed  in  great  detail  about  the  bivariate  cubic 
case.  It  can  be  used  to  fit  and  to  design  surfaces.  Some  results  are 
as  follows; 

The  system  of  equations  and  recurrence  formulae  for  determining 
the  unknown  partial  derivatives  in  its  expression  have  been  obtained 
under  certain  continuity  requirments. 

2.  The  property  that  its  convexities  can  be  determined  by  those  at 
the  knots  has  been  proved. 

3.  A  simple  method  for  finding  the  intersection  of  two  such  sur¬ 
faces  has  been  derived. 

4.  Each  of  its  rth  partial  derivatives  appraoches  the  corresponding 
one  of  those  of  the  fitted  surface  with  the  approximation  order  4:-r 
(r=0,l,2,3). 

This  kind  of  cubic  surfaces  has  been  used  to  fit  blue-print  sur¬ 
faces  successfully.  As  a  fxmction  of  this  kind  of  cubic  surface,  it 
has  been  used  to  solve  the  boundary  problems  of  linear  and  quasi- 
linear  P.D.S.  successfully.  ■ 


Approximation  of  curves  and  surfaces  by  the  exclusion  algorithm. 


Jean-Claude  Yakoubsohn 

Laboratoire  Analyse  Nnmerique,  Universite  Paul  Sabatier 
118,  route  de  Narbonne  31062  Toulouse  cedex. 


In  general  way  an  exclusion  algorithm  permits  to  localize  an  variety  of  i?".We  here  give  the  numerical  results 
concerning  curves  and  surfaces. We  describe  the  two  ingredients  which  are  used. 

1-  the  computation  of  the  “  best  “  exclusion  function. 

2-  the  geometrical  organization  of  the  algorithm. 

We  also  explain  how  we  represent  curves  and  surfaces. 


The  Optimal  Error  Bounds  For  The  Cubic  Spline  Interpolation 


Ye,  Maodong 

Dept,  of  Mathematics 
Zhejiang  University 
Hangzhou,  310027,  P,R.CHINA 


Abstract 


Let  S£(x)  be  the  cubic  spline  interixiiating  to  f(x)  on  the 

equidistant  knois  of  [  0,  1]  .  The  optimal  error  bound  of  the 
inteirpoiation  is  defind  by 

...  //t  /v  , 

[‘’•U 

We  give  the  explicit  expression  of  the  kenel  of  the  cubic  spline 
interpolation,  then  exactly  obtain  the  optimal  error  bounds  for  the 
cubic  spline  interpolation. 


B-Spline  Patches  for  Surface  Reconstruction 


ChangSLeng  Zhao  Roger  MOHR 


Lifia-Irimag  46,  Avenue  Felix  Viallet  38031  Grenoble,  France 
E-mail:  zhao@lifia.imag.fr.  Phone:  (33)  76574809,  Fax:  (33)  76574602 

This  paper  considers  the  problem  of  3D  surface  reconstruction  using  B-spline 
patches  in  computer  vision.  This  reconstruction  is  done  by  observing  the  occluding 
contours  where  the  view  lines  gaze  the  surface.  It  has  been  shown  that  reconstruction 
of  such  a  3D  surface  is  possible  when  the  camera  motion  is  known,  except  for  the  fully 
concave  parts  where  no  view  line  can  be  tangent  to  the  surface. 

Our  original  contributions  with  respect  to  these  related  works  are  mainly  twofold. 
Firstly,  we  introduce  a  direct  regularization  on  the  3D  surface  to  be  estimated  instead 
of  smoothing  the  contours  in  the  2D  image.  This  is  based  on  the  regularized  uniform 
bicubic  B-spline  patches.  Secondly,  we  propose  to  globally  recover  the  whole  surface 
from  small  local  patches.  The  others  have  only  been  content  with  the  recovery  of  some 
local  properties  of  the  surface  such  as  curvature  estimation. 

An  image  point  of  an  observed  occluding  contour  provides  exactly  three  inde¬ 
pendent  equations.  Two  of  them  state  that  the  surface  is  reached  by  the  view  line 
associated  with  the  image  coordinates.  The  third  equation  expresses  that  the  view  line 
is  tangent  to  the  contour  at  the  point  which  generates  a  plane  tangent  to  the  surface. 

Therefore,  each  contour  point  gives  a  set  of  three  non-linear  equations,  this  leads 
to  an  over-determined  system  of  non-linear  equations  for  all  observed  points  set  in 
the  image  sequence.  This  system  can  then  be  solved  by  a  non-linear  least  squares 
estimation,  for  example,  by  the  Levenberg-Marquardt’s  algorithm  which  is  currently 
implemented  in  our  system.  Parametrization  of  the  surface  patches  has  to  be  decided 
for  implementing  this  approach.  The  choice  we  made  is  very  close  to  the  epipolar 
parametrization  described  by  Blake  and  Cipolla.  The  starting  values  of  the  system  is 
provided  by  the  simple  method  of  triangulation. 

After  such  a  patch  is  estimated,  the  problem  remains  of  fusing  these  different 
patches  obtainted  with  different  motions.  This  is  done  by  going  back  to  a  discrete 
representation  with  a  net  of  3D  points. 

We  have  experimented  on  some  simulated  smooth  objects.  The  reconstruction 
process  provided  good  results  and  its  numerical  behaviour  was  studied  by  adding 
different  pixel  noise  to  image  points.  We  showed  also  how  the  different  parameters 
should  be  conveniently  selected  according  to  available  data  through  experiments.  The 
experiments  on  real  image  sequence  are  undergoing. 
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